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PREFACE 


T his book is intended as an introduction to some recent 
developments of Maxwells electromagnetic theory which 
are directly connected with the solution of the partial differential 
equation of wave-motion. The higher developments of the 
theory which are based on the dynamical equations of motion are 
not considered at all. Even with this limitation the subject is a 
vast one, and to bring the work of perusing the literature within 
my power I have omitted an account of the modem theory of 
relativity which has been expounded very clearly in several recent 
publications. 

For a thorough understanding of the present subject a very 
extensive knowledge of mathematics is necessary, but there are 
parts of the subject in which a reader with only a limited 
mathematical equipment may soon feel at home and perhaps do 
useful original work. With the idea of enabling such a reader to 
obtain a quick grasp of the nature of the subject and the results 
obtained, I have thought it advisable to state without proof a 
number of relations of which adequate demonstrations can only 
be obtained by means of complicated and difficult analysis. 
I have also endeavoured to keep the analysis as elementary as 
possible, but in some places where the work is perfectly straight- 
forward a few details are omitted. 

The book is far from being a complete treatise on the subject, 
for I have not given any existence theorems to show that the 
solutions of certain problems exist and are unique, and no 
attempt has been made to enter into the details of numerical 
computations. There are many parts of the subject indeed to 
which a pure mathematician might make useful additions; in 
particular, I might direct attention to p. 21, line 2, and p. 101, 
where there are one or two matters which require further 
discussion. 
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PHEFACE 


Chapter viii and paragraph 5 contain some of my own con- 
tributions to the subject. At present there seem to be several 
different directions in which future developments may be made, 
and so it seems unwise to give a hasty judgment concerning the 
physical significance of the results. Ideas which naturally 
present themselves are that the aether can be regarded as built 
up from singular curves of the type considered in § 43, and that 
§§ 41 and 44 may throw some light on the question of the difference 
between positive and negative elementary electric charges. I 
hope to discuss an hypothesis relating to the first idea in a future 
note, but am unable to give any support at present to the 
second idea. 

I gratefully acknowledge my indebtedness to Sir Joseph 
Larmor,who read the manuscript before it was revised and made 
some helpful suggestions, to Prof. Ames who read the greater 
portion of the manuscript, to Prof Morley and Mr Hass^ who 
helped me with their advice and vigilance in reading the proof- 
sheets, and to the officers and staff of the University Press for 
their careful work and constant consideration shown in matters 
connected with the printing. For the correctness of the new 
formulae and examples I alone am nvsponsible ; if any errors are 
discovered I shall be grateful if my readers will inform me. 

HARRY BATEMAN. 


October^ 1914. 
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CHAPTER I 


FUNDAMENTAL IDEAS 


§ L The fundamental equations for free aether. 

In Maxwell’s electromagnetic theory the state of the aether 
in the vicinity of a point y, z) at time t is specified by means 
of two vectors E and H which satisfy the circuital relations* 


rot H = 


IdE 

c dt^ 


rot E — — 


1 dH 
c dt 


( 1 ), 


and thc‘ solenoidal or sourceless conditions 


div /?= 0, (liv = 0. 

If right-handed rectangular axes are used the symbolf rot H 
(l<‘notes the vector whose components arc of type 

dH, _ dHy 

dy dz ’ 

the three components of H being Hy, H, respectively. 
The symbol div H denotes the divergence of if, i.e. the 
(juantity 

dHy dH, 
dx dy dz 

The vector E is called the electric displacement or electric 
force and H the magnetic force. The quantity c represents the 


^ The equations are written in the symmetrical form in which they were 
presented by 0. Heaviside, Electrical Papers, Vol. 1, § 30, and H. Hertz, 
Electric Waves, p. 138. Sir Joseph Larmor points out that a set of equations 
equivalent to these was hrst used by MacCullagh in 1838 as a scheme 
consistently covering the whole ground of Physical Optics, Collected Workt of 
James MacCullagh (1880), p. 145. 

t We use here the units and notation employed in Lorentz^s Tlie Theory of 
Electrons, Gh. i, except that large letters are used to denote vectors and E is 
written in place of D. Many writers use the symbol curl instead of rot. 
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velocity of propagation of homogeneous plane waves and is 
commonly called the velocity of light ; we shall assume it to be 
a constant, although in the most recent speculations it is treated 
as variable*. 

Some of the modern writers on the theory of relativity 
maintain that the introduction of the idea of an aether is 
unnecessary and misleading. Their criticisms are directed 
chiefly against the popular conception of the aether as a kind 
of fluid or elastic solid which can be regarded as practically 
stationary while material and electrified particles move through 
it. This idea has been very helpful as it presents us with 
a vivid picture of the processes which may be supposed to take 
place, it also has the advantage that with its aid we can attach 
a meaning to the term absolute motion, but herein lies its 
weakness. Larmor, Lorentz and Einstein have shown, in fact, 
that the differential equations of the electron theory admit of 
a group of transformations which can be interpreted to mean 
that there is no such thing as absolute motion. 

If this be admitted, the popular idea of the aether must be 
regarded as incorrect, and so if we wish to retain the idea of a 
continuous medium to explain action at a distance we must 
frankly acknowledge that the simplest description we can give 
of the properties of our medium is that embodied in the 
differential equations (1). 

If we abandon the idea of a continuous medium in the 
usual sense only two ways of explaining action at a distance 
readily suggest themselves. We may either think of the 
aether as a collection of tubes or filaments attached to the 
particles of matter as in the form of Faraday s theory which has 
been developed by Sir Joseph Thomson and N. R. Campbell ; 
or we may suppose that some particle or entity which belonged 
to an active body at time t belongs to the body acted upon at a 
later time ^ + t. From one point of view these two theories are 
the same, for if particles are continually emitted from an active 

* A. Einstein, Ann, d. Phys. Vol. 35 (1911), p. 898 ; Vol. 38 (1912), pp. 366 
and 443. M. Abraham, Phys. Zeitsekr. (1912), pp. 1 — 6, 310 — 314, 793—797 ; 
Ann, d, Phys, (1912), pp. 444 and 1066 ; Fifth International Congress of Mathe< 
maticians, Proceedings^ Vol. 2, p. 266. 
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body they will form a kind of thread attached to it. The first 
form of the theory is, however, more general than the second. 

At present we are unable to form a satisfactory picture of 
the processes that give rise to, or are represented by, the vectors 
E and H, We believe, however, that some points may be made 
clear by studying the properties of solutions of our differential 
equations. 

It will be seen from the investigations of Chapter viii 
that the mathematical analysis connected with these equations 
is suitable for the discussion of three distinct theories of the 
universe, which may be described briefly as follows : — 

Aether Matter 

Continuous medium. Aggregates of discrete par- 

ticles. 

Discontinuous medium con- An aggregate of discrete 

sisting of a collection of tubes particles attached to the tubes, 
or filaments. 

Continuous medium. An aggregate of discrete 

particles to which tubes are 
attached. 

The last theory may be supposed to include that form of 
the emission theory of light in which small entities are projected 
from the particles of matter under certain circumstances and 
produce waves in the surrounding medium. This theory might 
be justly ascribed to Newton*. 

For other theories of the aether the reader is referred to 
Prof. E. T. Whittakers recent workf A History of the Theories 
of the Aether, 

In the first part of this book the analysis is adapted almost 
entirely to the first theory, the high development of which we owe 
to the pioneer work of Maxwell, FitzGerald, Hertz, Rayleigh, 
Heaviside, J. J. Thomson, Lorentz and Larmor. The other 
theories have not yet received much attention but it is hoped 

* A form of the theory in which the entities are electric doublets has been 
developed by W. H. Bragg and applied to the X and y rays. British Association 
Beports (1911), p. 840. 

t Dublin Univ. Press ; Longmans, Green and Co. (1910). 
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that the analysis of Chapter viii will lead to further develop- 
ments so that a comparison can be made between the different 
theories. It is quite likely that one theory will be enriched by 
the developments of another. 

§ 2. Electromagnetic fields. 

For many purposes it is convenient to work with a complex 
vector* * M^H ±iE, where i = V — 1 and the ambiguous sign 
+ is independent of the ambiguity which occurs in the determi- 
nation of V — 1. The differential equations (1) may then be 
replaced by the simpler equations 

rot iV = + - ^ , (liv iV/ = 0 (2). 

cdt 

When a solution of these equations has been found a pair 
of vectors E and H satisfying equations (1) may be obtained by 
equating coefficients of the ambiguous sign. In working with 
an ambiguous sign it must be remembered that when two 
ambiguous signs are multiplied together the ambiguity is 
removed. The chief advantage in using the two independent 
ambiguities + and V — 1 is that we can assume that the vectors 
E and H are the real parts of expressions of the form Ae"'*** and 
we are at liberty to equate the coefficients of either i or ± in 
any of our equations. 

Definition, A solution of the differential equations (2) or (1), 
which provides us with single-valued vector functions E and H 
for each space-time point {x, y, z, i) belonging to a certain 
domain D, is said to define an electromagnetic field in the 
domain D, 

Since the differential equations are linear the sum of any 
number of solutions is also a solution. The physical meaning 
of this is that when two electromagnetic fields are superposed, 
they are together equivalent to an electromagnetic field. 

Two superposed electromagnetic fields can of course be 
related to one another in some way. When electromagnetic 

* The use of a complex vector recommended by L. Silberstein, Ann, 

d, Phys, Vols. 22 and 24 (1907) ; Phil, Mag, (6), Vol. 23 (1912), p. 790. He does 

not, however, use the ambiguous sign. 
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waves fall upon an obstacle, a secondary disturbance is produced 
which depends in character upon the nature of both the primary 
waves and the obstacle. 

We shall find that in some cases it is possible to find two 
fields in which the vectors {Ey i/), {E\ H') are connected by the 
two relations embodied in the equation 

{MM') = MJf^' + MyMy' + M,M; = 0 (3) 

for all values of {x, y, z, t) belonging to some domain. 

When this is the case the fields are said to be conjugate 
within this domain. 

If we use the notation 

we may write 

= (/f^) - {E^) ± 2i {EH) = ± 2i 

where /, and I., are two (quantities which we shall call the 
invariants*. It is easy to see that when two conjugate fields 
are superposed the invariant Ii for the total field is the sum of 
the invariants /j for the two component fields. Similarly for 
the invariant 1.^, 

Whc*n the invariants are zero over a given domain the field 
may be called self-conjugate for this region f. 

§ 3. The flow of energy. 

An entity whose volume density :J: p is a function of (.r, y, Zy t) 
will vary in a mannei- which can be described as a simple flow 
with component velocities (w-, Vy w) if the equation of continuity 



is satisfied. This equation implies in fact that there is no 

* They are invariants for the group of linear transformations which leave 
the electromagnetic equations unaltered in form. Gf. H. Minkowski, 
Nachr. (1908) ; E. Cunningham, Proc. London Math. Soc. (2), Vol. 8 (1910), 
p. 89 ; H. Poincard, Rend. Palermo (1906) ; M. Planck, Ann. d. Phys. Vol. 26 
(1908). Other invariants are given by these authors, 
t Silberstein calls it a pure electromagnetic wave. 

X The limitations to which the idea of density is subject and the question of 
the continuity of the function p are discussed by J. G. Leathern, “Volume 
integrals and their use in physics,** Cambridge Mathematical Tracts (1905). 
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creation or annihilation of the entity in the neighbourhood of 
{x, y, z, f). 

Now it is easy to see that the equation of continuity is 
satisfied in virtue of equations (1) if we put 

= + pu^ciEyH^^E^Hyl 

and two similar equations. We shall regard p in this case as 
the volume density of the energy contained in the electn^- 
magnetic field. The vector S whose components are of the 
type c (EyHz — E^Hy) can then be supposed to indicate the rate 
at which energy flows through the field. Since 

p2 (c® — XL‘ — w'^) = ^ (f {E^ — H^Y + (EHYf 
it appears that energy travels through the field with a velocity 
which is less than the velocity of light. The velocity c is 
attained only in the case of a self-conjugate field. 

The vector S was introduced by Prof. Poynting* and is 
usually called Poyn ting’s vector. The idea of describing the 
transfer of energy in this way also occurred to Prof. Lamb 
before the publication of Poynting s work. 


Example, Prove that the equation of continuity may he satished hy 
putting 

ld6 „ dS 
1 a^ „ 3(9 


f^,V 



i OB „ zr , tj 


d6 d0 d0 


where B is an arbitrary function. Obtain a similar solution by replacing 
E by if and if by - E, 


§ 4. First solution of the fundamental equations. 

Let us use the symbol flw to denote the Dalembertianf 
of w, viz. 

^ . Id^u ^ dhi d^it 1 

* Phil, Trans, A, Vol. 175 (1884), p. 343. See also H. A. Lorentz, The 
Theory of Electrons, p. 22. 

t This is the name suggested by Lorentz, loc. cit. p. 17. Many writers use 
Cauchy’s symbol □ to denote the Dalembertian, but 1 think Q is preferable 
because its form suggests a wave. Murphy’s symbol A is also used here in place 
of the usual symbol V®. E. B. Wilson and G. N. Lewis use the symbol 
to denote the Dalembertian of u. Cf. Proc, Amer. Acad, of Arts and Sciences, 
Vol. 48 (1912), p. 389. 
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and the symbol grad U to denote the vector whose components 

m dU dU 

<11 G 

OX dy oz 

respectively. Let us also use flA, where A is a vector with 
components Ay, A^, to denote the vector whose components 
are HAa-, fiA^, HA^. The equation ^?^ = 0 will be called the 
wave-equation and a solution of* this equation a wave-function, 
A vector function A will be said to satisfy the wave-equation 
when each of its components is a wave-function, i.e. if HA = 0. 
We m<ay now satisfy equations (1) and (2) by writing 

M = ±i rot L= i -h grad A (5), 

where the scalar potential A = and the vector potential 

L — BTiA satisfy the ecjuations 

fiA=0, nL = 0, div/, + -^ = 0 (6). 

C Cv 


The Last three equations may be solved in a general way by 
writing 

L = ± t rot G + grad 


c dt 

A = ~ div Cr - - 

c at 


.(7), 


where the vector (? = F + iTI and the scalar K satisfy the wave- 
equation 

0 ^ = 0 ( 8 ). 

The solution of cqu*ations (1) which is embodied in (5), 
(6) and (7) is a simple extension of Hertzs solution* and is 
suggested by Whittaker’s solution f in terms of two scalar 
potentials. It is cle<ar that the function K drops out when we 
differentiate to find M and so the electric and magnetic forces 
depend only on the vector G, The form of this vector indicates 
that the electromagnetic field can be regarded as the sum of 
two partial fields ; one of these is derived from the vector 11 and 


* yinn, d, Phyii. Vol. 36 (18S8), p. 1. The general solution is given by 
Bighi, Bologna Mem. (6), t. 9 (1901), p. 1 ; IlNmvo Ghnento (5), t. 2 (1901), p. 2. 
He finds suitable expressions for the vectors II and F in a number of oases, 
t Proc. London Math. Soc. Ser. 2, Vol. 1 (1903). 
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will be called a field of electric type, the other is derived from 
the function T and will be called a field of magnetic type. 

This resolution of an electromagnetic field into two partial 
fields is analogous to the one used by H. M. Macdonald* in the 
study of the effect of an obstacle on a train of electric waves. 
The component fields are then of such a type that in one case) 
the magnetic force normal to the obstacle vanishes ovei* th<i 
surface of the latter, in the other case it is the electric force 
normal to the obstacle that vanishes. The same idea has been 
used recently by Micf and Debye J in the treatment of the case 
of a spherical obstacle. 

In Hertz’s solution we have F = 0, = 0 and H has 

components (0, 0, aS'). 

The components of E and H are consequently given by the 
formulae 




dxdz ’ 




1 0-N 
c dy?t 


F ~ 


ff 

” c dxbt I 


//. = 0 


.( 9 ). 


Hertz uses Euler’s wave-function § 


»S' = sin K (r — ct\ 
r 


7-2 == .7.2 + J/- 4. ^2^ 


and obtains in this way a theory of his oscillator||. The electric 
and magnetic forces become infinite at the origin which is there- 
fore a singularity of the electromagnetic field. A singularity of 
this type is called a vibrating electric doublet and is regarded as 
the simplest model of a source of light or electromagnetic waves. 

* Electric Waves, Ch. vi. 

t Ann. d. Phys. Vol. 25 (1908), p. 382. + Ibid. Vol. 30 (1909), p. 57. 

§ Periodic solutions representing a disturbance sent out from n-fold poles 
had been used previously by H. A. Rowland and applied to the elucidation of 
optical phenomena. Amer. Journal of Mathematics, Vol. 6, p. 359 ; Phil. Mag. 
Vol. 17 (1884), p. 423. Of. also Stokes, Cambr. Phil. Trans. (1849). 

II To deal with the case in which the vibrations are damped we assume 

sin K (r-ct). Cf. K. Pearson and A. Lee, Phil. Trans. A, Vol. 193 
r 

(1900), p. 159. 
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The solutions of equations (1) which are obtained by superposing 
elementary solutions of this type are of great importance in 
physical optics. 

When r is very great the most important terms in the 


expressions (9) are 



^* = - 

K-XZS 

II 


r- 

E„ = - 

fC-jfZS 

r* ' 



+ y) s 

--.s ^ 

//. = o. 


where a = sin k (r — ct). All the other terms are of order 1// - or 
l/r\ These expressions give 

(£^)-(//0 = 0 . 

Hence at a vtiry great distance from the origin the field is 
practically a self-conjugate field and so the energy travels with 
a velocity very nearly equal to the velocity of light. The 
expressions indicate that Poynting's vector is ultimately along 
tht» radius from the origin ; now the electric and magnetic 
forces are at right angles to Poynting's vector and so the 
vibrations of the light-vector, whether we take it to be the 
electric or magnetic force, are at right angles to the radius. 
The waves sent out from the source have, then, the character 
of monochromatic light at a great distance from the origin*. 
The amplitudes of the vibrations at points on the same radius 
are proportional to the quantities 1/?’ when r is large, and so if 
the intensity of the light be measured by the square of the 
amplitude the inverse square law is fulfilled. 

Since the electric force is ultimately at right angles to the 
radius there is no total charge associated with the singularity, 
for the charge is equal to the surface-integral of the normal 
electric force over a large sphere concentric with the origin 
and this integral is evidently zero. We are consequently 
justified in regarding the singularity as a doublet and in fact 

* For a fuller discussion see Larmor, Phil, Mag, (5), Vol. 44 (1897), p. 503 ; 
Aether and Matter, Chap, xiv, where it is shown that energy is radiated from a 
moving charge only when the velocity of the charge alters in either magnitude 
or direction. 
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as a simple electric doublet of varying moment as is indicated 
by the way in which the electric and magnetic forces become 
infinite*. The axis of the doublet is along the axis of z. 

The electric lines of force due to a vibrating electric doublet 
have been drawn by Hertzf for various stages of the motion. 
The general character of the lines of force is indicated in Fig. 1. 



It will be noticed that the lines are all at right angles to a 
plane perpendicular to the axis of the doublet. M. Abraham § 
has used a Hertzian doublet to obtain a model of the electro- 
magnetic field produced by the oscillations in a vertical antenna, 
the plane just mentioned being supposed to represent the earth 
which is regarded as a perfect conductor. Zenneck|| has, 
however, pointed out that when the imperfect conductivity of 
the earth is taken into account the circumstances of the 

* See § 42. 

t Ann. d. Phy$. Vol. 36 (1888), p. 1. The case of damped vibrations is 
considered by K. Pearson and A. Lee, loc. cit. 

1 am indebted to the Macmillan Company and A. Gray, Esq., for permis- 
sion to reproduce this diagram. 

§ Phys. Zeiischr. Vol. 2 (1901), p. 329 ; Theorie der Elektrizitiit, Vol. 2, 
§ 34; Encyklop. d. Math. Wiss. Band 6, § 18. 

li Ann. d. Phys. Vol. 23 (1907), p. 846 ; Phys. Zeitschr. Vol. 9 (1908), p. 50 ; 
Ibid. p. 553. 
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propagation are somewhat different. The spreading of electro- 
magnetic waves over the earth’s surface has been investigated 
thoroughly by A. Sommerfeld* * * § and his pupil H. v. Hoerschle- 
mannf, and their results seem to indicate that the imperfect 
conductivity of the earth is an important factor in directing 
electric waves and in enabling their effects to be detected at 
great distances. The ionisation of the air by sunlight is also 
an important factor, as has been pointed out by J. J. Thomson, 
W. H. Eccles J and J. A. Fleming§. Marconi’s experiments have 
indicated that the circumstances of propagation are not yet 
thoroughly understood. No good reason has been given to 
explain why communications by means of electric waves can be 
made more easily when the receiving station is in a north or 
south direction than when the direction is east or west. The 
curious contrasts in the results obtained with waves of different 
frequencies in day and night communications are also un- 
i‘xplained||. 

The use of the vector 11 instead of the scalar 8 was 
recommended by Abraham IT. Von Hoerschlemann has obtained 
in this way a model of Marconi’s bent antenna which gives 
a directed effect to the radiation. A number of arrangements 
of Hertzian doublets that can be used to imitate the action 
of antennae have been described by Fleming** * * §§ , Larmorf"|*, 
Sominerfeld and SlacdonaldJJ. 

In the theory of FitzGerald’s magnetic oscillator§§ we have 

11 = 0, r = (0,0,^), * 

iV being Euler’s wave-function. Whittaker’s solution is 
obtained by adding the solutions of Hertz and FitzGerald. 

* Ann. d. Phys. Vol. 28 (1909), p. 665. 

\ Jahrb. d. draht. Teley. Vol. 5 (1912). 

X Proc. Roy. Soc. A, Vol. 87, p. 79. 

§ British Association ReportSj Dundee (1912). See also O. J. Lodge, 
Phn. Mag. Vol. 25 (1913), p. 775. 

II See Marconi’s address to the Royal Institution, June, 1913. 

II Theorie der FAektrizttfit, Vol. 2, Ch. i. See also Righi, loc. cit. 

** Proc. Roy. Soc. A, Vol. 78, p. 1. 

+t Ibid, in a footnote to Fleming’s paper. 

XX Proc. Roy. Soc. A, Vol. 81, p. 394. 

§§ Tram. Roy. Dublin Soc. Vol. 3 (1883); Scientific Writings, p. 122. 
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§ 5. Second solution of the fundamental equations. 

It is easy to see that equations (2) will be satisfied if we 
can find two functions (a, /3) such that 

M 4 - i 

0 {y, z) ~ c d {x, t) 

i d(a, 0) 

^ d (z, x) “ c 0 (y, t) 

^ d(a, 0 )^ i 9 (a, 

" 0(.r, 2 /) ~'c d{z,t) 

An electromagnetic field that is specified in this way is 
necessarily a self-conjugate field, for if we multiply together 
the two expressions for and do the same for My, M^, we find 
that M^ = 0. A particular pair of functions a, 0 is obtained by 
putting 

a — x cos 6 + y sin 6 iz, 0 — xsinO — y cos 6 ^ct (11), 

where 6 is an arbitrary constant. To generalise this field wt‘ 
multiply the expressions for Mx, My, Mz by an arbitrary 
function* of a, y3, ^ and integrate with regard to 6: we thus 
obtain a very general electromagnetic field in which 


Mx — + i j 

r 2rr 

f{a,^,6)cos6d$\ 


1 

My = 4- j 

1 /(a, 0) sin 9 dO 

0 



= - J 

I f(ci,fi,e)de } 
0 


The components of the electric and magnetic forces are 
obtained by equating the ambiguous and unambiguous parts in 
these equations; it is easy to verify that they are all wave- 
functions. 

It should be remarked that these definite integrals may 
give a representation of the electromagnetic field, required for 
the solution of a problem, only in a certain limited domain of 

* When we speak of an arbitrary function it must be understood that the 
function may be subject to certain limitations which render the integration and 
differentiation under the integral sign intelligible operations. 
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the variables x, y, z,t\ the integi*als may in fact represent 
discontinuous functions. 

The limits of integration could have been taken to be any 
other constants instead of 0 and 27r ; they can also be taken to 
be functions of Xy y, Zy t of the type o), where o) is defined by an 
equation of the form 

X sin ft) — y cos w — ct = F{w)y 
F being an arbitrary function. 

A suitable ])air of functions a, /8 is also obtained by putting 

^=:r~-ct (13), 

and in this case Poynting’s vector is along the radius from the 
origin. A more general type of electromagnetic field in which 
this is true is obtained by multiplying the above expressions 
for the ct)mponents of M by an arbitrary function of a and 

Other pairs of functions a, ^ of a very general nature are 
obtained in Chap. viii. It should be remarked that in all cases 
the functions (a, yS) are of such a nature that if F (a, is an 
arbitrary function of a and /3, F satisfies th® partial differential 
(Mjuation 



which is of fundamental importance in geometrical optics* and 
may be called Hamiltons equation. It is found that this 
(‘quation is also satisfied in many cases by the functions of 
,T, y, Zy t which are the limits of a definite integral representing 
a wave-function, when the function under the integral sign is 
a wave-function for all values of the parameter with regard to 
which we are integrating. Thus the function co just defined and 

the function i ~ - (r -f ro) which will be used later are solutions 
c 

of this equation. 


* For another connection between this equation and the electromagnetic 
equations see A. Sommerfeld and J. Bunge, ** Grundlagen der geometrischen 

Optik,** Ann. d. Phys. Vol. 36 (1911), p. 277. 
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§ 6. The fundamental equations for a material medium. 

For a material medium which is stationary relative to the 
axes of coordinates, the equations ( 1 ) must be replaced by the 
more general equations* 


rot ^ , div H = 0 

c dt 


(15), 


where i) is the electric displacement, E the electric force or 
field strength, H the magnetic force and B the magnetic induc- 

dJD 

tion. The quantity J + represents the total current which 

is made up of a conduction-current (7, a displacement-current 
dD ' 

and a convection-current pv, p being the volume density of 


electricity. 

Various notations have been used for the different vectors 
of an electromagnetic field. Most English writers use (a, 6 , c) 
for the components of the magnetic induction, (a, 7 ) for those 

of the magnetic force, (/, h) for the components of the electric 
displacement and (P, Q, R) or (X, F, Z) for those of the 
electromotive intensity or electric force f. This is not to be 
confused with the mechanical force F of electromagnetic origin, 
whose components are sometimes denoted by (X, F, Z). 


* Lorentz (1892 — 1895) and Larmor (1895) have derived these equations and 
a corresponding set of equations for moving bodies by a process of averaging, 
starting from the fundamental equations of the theory of electrons in which we 
have B = H, D—E^ J—pv. Cf. H. A. Lorentz, A Jtad. van Wetenschappen te 
Amsterdam (1902), p. 305; Encykl. d. Math. Wiss. Bd. 5, § 14, pp. 200 — 210. This 
method of averaging has been developed so as to give results in accordance with 
the Theory of Relativity by M. Born, Math. Ann. Bd. 68 (1910) andE. Cunning- 
ham, Proc. London Math. Soc. Ser. 2, Vol. 10 (1911), p. 116. The Born- 
Minkowski equations differ slightly from those of Lorentz and indicate the 
existence of an electrostatic held due to the motion of a magnetised body. 

It has been realised by the foregoing writers and others that the principle of 
relativity alone is not sufficient to determine a complete set of equations for 
moving bodies, a theory of the constitution of matter is needed. Cf. H. R. Hass4, 
Phil. Mag. Jan. (1914). 

+ Clerk Maxwell, Electricity and Magnetism^ 3rd edition (1892), Vol. 2, 
p. 257. 
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In any material medium there are certain constitutive 
relations connecting the vectors i), Ey By H, J, In moving 
media, crystalline media and ferromagnetic bodies the relations 
are rather complicated, but for an isotropic medium in which 
p = 0 the relations can be represented to a good degree of 
approximation by the simple equations 

D = eEy J = crEy B = fjuH ( 16 )» 

where e, o', are scalar quantities which are generally regarded 
as constants; they will be regarded in fact as the optical 
constants of the medium. The quantity a is called the cm- 
d activity y /x the permeability y and e the dielectric-inductive 
capacity. 

The units that are used here are the so-called modified 

•> 

units*, in which Heaviside’s suggestion of eliminating a factor 
47r has been adopted. We can pass to electrostatic units or 
electromagnetic units by replacing our quantities Ey Hy etc. by 
aEy 0Hy etc., where a, /3 are certain factors which are given in 
the following table : 



e, /J, J 

E 

B 

H 

m 

Electrostatic system 

\/47r 

1 

sj 4’*' 

i 

c 

isJAtt * 

1 

c fj An 

c ijAn 

\ 

Electromagnetic system 

C /v/4Tr 

1 

c ij An 

1 

JAtt 

1 

jAn 



We use e here to denote a quantity of electricity, and m a 
quantity of magnetism. 

• 

§ 7. The energy equation for a material medium. 

If we use 2 as before to denote the vector whose components 
are of type c {EyHz — E^H^y we find that 


dx 9y dz 


* Cf. H. A. Lorentz, EneyklopHdie der Math. JFiss. Bd. 5, § 13, pp. 83 — 87. 
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If /A is a constant and + be regarded as the 

energy per unit volume, the change in the distribution of the 
energy can be described by means of a flow a and a loss per 
unit volume of magnitude <tE^ due to the transformation of 
electric energy into heat (Joule's heat)*. If B does not depend 
on the instantaneous value of H so that is not a constant 
there is a loss of energy due to hysteresis. B may depend 
upon H alone but not be a single-valued function of //, con- 
sequently in a cycle of changes ^{H,dB) is not zero and may be 
taken as the heat per unit volume developed during the 
description of the cycle. Notice that 

1{H .dB) = — /(i^ . dH) round a cycle, 
and is always positive since the value of B for a given value of 
H is greater when H is increasing. The experimental analysis 
and the accompanying theory are due to E. Warburgf and 
independently in much greater development to J. A. EwingJ by 
whom the name hysteresis was applied to such phenomena. 

§ 8. Solution of the fundamental equations for a 
material medium. 

Let us assume that cr, /i, e are constants and that E, H are 
the real parts of expressions of the form where A is 

a complex quantity independent of t Then if we write 

+ i/aoxr kc ,, zr . • c- 

k" = - — ' — , = — — , jyl = iM + %v E (17 ), 

C- 

and regard Ey H now as the complex v(ictors of which they were 
formerly the real parts, the differential ecpiations to be satisfied 


by M are 

rot Jl/ = + kMy div 0 (18). 

These may be solved by putting 

M = rot n ± ^ grad div TI ± fell (19), 


where 11 is a solution of the equation Aa -\-k^ii = 0, and may bt^ 
of the form U ±iV, 

* For a fuller discussion see R. Gans, Einfilhrunp in die Theorie der 
Magnetismus; H. A. Lorentz, Encykl. d. Math, Whs. Bd. 6, § 14, p. 240, Heft 1 
(1903) ; Heaviside, Electrical Papers^ Vol. 1, pp. 437 — 450. 

t Ann. Phys. Chem. (3), Vol. 13 (1881), p. 141. 

J Phil. Trans. A, Vol. 176 (1885), pp. 623—640; Proc. Roy. Soc. Vol. 34 
(1883), p. 39; Phil. Mag. Vol. 16 (1883), p. 381; Magnetic Induction in Iron 
and other Metals y London (1892). 
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Problems which deal with the effect of small obstacles upon 
light or electric waves require the taking into account of the 
properties of the material of which the obstacle is composed : it 
is only for long waves and conducting media that the obstacle 
can be treated as a perfect reflector. This is illustrated by the 
work of Maxwell Garnett* and Mief on the optical properties 
of colloidal suspensions of metals and in the work of Sommerfeld 
to which we have already referred. 

Unfortunately, however, the analytical difficulties are very 
great when imperfect conductivity is taken into account The 
simple-looking problem of the reflection of the disturbance 
produced by a moving charge, when the obstacle is an infinite 
plane sheet of metal or other conducting substance, has not yet 
been solved accurately § and there are many similar problems 
that have completely baffled mathematicians. 

Much more progress has been made with problems dealing 
with perfect reflectors. These problems are to some extent 
ideal but some of the characteristics of actual physical problems 
are often preserved ||. Apart from this, such problems are of 
considerable mathematical interest, and have been studied by 
some writers simply on this account. 

§ 9. Boundary-Conditions. 

The conditions to be satisfied at a surface separating two 
different media are obtained by integrating equtCtions (15) across 
a thin layer of transition*^. Taking the axis of z along the 

* mi, Trmut. A, Vol. 203 (1904), p. 385 ; Vol. 205 (1905), p. 237. 

t a/iH. d, riujH. Vol. 25 (1908), p. 377. 

X An importimt solution of the equations has been given by M. Brillouin, 
“Propagation dans les milieux conducteurs,” Comptes Remius, t. 136 (1903), 
pp. 667, 746. (See Ex. 20, Cb. ix.) For other references see Ex. 5, p. 23. 

§ An approximate solution was suggested by Maxwell and has been developed 
by Larmor, Proc. London Math. Soc.^ Ser. 2, Vol. 8, p. 1. An accurate solution 
for the case in which the sheet is treated as infinitely thin and the charge moves 
with uniform velocity parallel to the plane has been given by G. Picciati, Rom. 
Acc. Line. Rend. (5), lla (1902), p. 221. 

II A metal behaves practically as a perfect conductor to electric waves when 
the oscillations are rapid but slow compared with waves of light. Cf. J. Larmor, 
“Electric vibrations in condensing systems,” Proc. London Math. Soc. Ser. 1, 
Vol. 26, p. 119. 

IT Rayleigh, Scientific Papers, Vol. 1; Phil. Mag. Vol. 12 (1881), p. 81; 
H. Hertz, Electric Waves, pp. 207, 238 ; Larmor, Phil. Trans. (1895), 


B. 


2 
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normal to the surface we shall assume that E, H, D, B and 
their derivatives with regard to oo, y, t are finite within the 
layer and that the conductivity <r is also finite. The equations 


dj^ 

dz 


dy 


dH^ 


dz 




then show that 


dy c dt ' 

and Jx 4* are finite and so their 
dz dz dt 

integrals with regard to z across a thin layer of thickness 0 are 

less than ad where a is a finite positive quantity independent of 0. 

This means that the tangential components of the electric force, 

magnetic force and electric current are continuous in crossing 

the surface. 

Again, the equation 


dBx . dBy . dB, 


+ 


dz 


= 0 


dB, 


dx dy 

shows that is finite and so the normal component of the 

magnetic induction is continuous. This result may, however, 
be regarded as a consequence of the previous one. 

The normal component of the electric displacement may be 
discontinuous, for the equation div D ^ p gives 
f^dD 

= j dz = j pdz 4 terms of order 0, 

where dz is the discontinuity in the electric displacement. Hence 
if a is the surface-charge of electricity per unit area we have dz=a. 
When the media are both conducting we have o' = 0 and the 
normal component of the electric displacement is continuous. 

It is known that in the case of a good conducting body 
rapidly alternating currents are confined within a very thin 
layer close to the surface*. In the ideal case of a perfect 
conductor or perfect reflector the field-vectors are zero within 


p. 733; H. M. Macdonald, Electric Waves, p. 14; H. A. Lorentz, Zeitschrift 
f. Math. u. Phys. Bd. 22 (1877). 

* Of. H. Lamb, Phil. Trans. A (1883) ; O. Heaviside, Electrical Papers^ Vol. 2, 
p. 168 ; J. J. Thomson, Recent Researches, p. 281. For some recent work on 
the subject see a paper by E. F. Northrop and J. E. Carson, Journ. of the 
Franklin Institute, Feb. (1914), p. 125. The results of many other researches 
on the skin-effect and alternating current resistance are given in J. A. Fleming’s 
Ths principles of Electric Wave Telegraphy. 
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the body of the conductor. This case is characterised by or = oo 
and the tangential components of the magnetic force are no 
longer continuous as a point moves across the surface, we have 
in fact if h is the discontinuity in the magnetic force 

K = / dz^ — I crExdz + terms of order d, 

Jo dz cJo 

Hence if K is the surface-current and h the discontinuity 
of the magnetic force, we have 

K = ]Ky. 

0 o 

At the surface of a perfect conductor the tangential com- 
ponents of the electric force must vanish and as a consequence 
of this we can say that the normal component of the magnetic 
induction must also vanish. If the medium outside the con- 
ductor is free aether the surface-conditions are simplified on 
account of the relations D = E, B — H, 

In the case of a very thin conducting sheet it is convenient 
to treat the thickness of the sheet as negligible and regard the 
tangential components of the magnetic force as discontinuous 
when a point moves along the normal from one side of the 
sheet to the other. 

The boundary-condition is then* 

— hy = — Exy hx = 

where ^ = J <^dz is the conductivity of the sheet. 

If we wish to extend the idea of Green’s equivalent layer to 
electrodynamics we must consider electromagnetic fields in free 
aether with surfaces at which the tangential components of the 
electric and magnetic forces are discontinuous ; this requires an 
electric current sheet and a magnetic current sheet on the 

* Cf. T. Levi-Civitfl, Rend, Lined (6), llg (1902), p. 75. These conditions 
are used by Picoiati in his solution of the problem of an electric charge moving 
parallel to a conducting sheet. In some previous papers, Rend, Lincei (5), lli 
(1902), pp. 163, 191, 228, Levi-CiviU had used the electromagnetic potentials to 
determine the effect of a conducting sheet on an alternating current flowing 
along a straight wire parallel to the sheet; the boundary-conditions are then 
determined by the discontinuities of the potentials due to the induced current. 

^2 
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surface*. For a complete generalisation we ought to consider 
the cases when the normal components are also discontinuous 
and when the surface is in motion. The last circumstance 
alters matters to some extent and must now be discussed. 

In the case of electric waves in free aether, the vectors 
E and H may be discontinuous at a wave-boundary. If this 
can be regarded as the limit of a thin layer of transition within 
which equations (1) are satisfied and the vectors Ey H are 
finitef , the values of these vectors on the two sides of the 
boundary must satisfy certain conditions which may be found 
as follows. 

Let the equation of the moving boundary be expressed in 
the form 

( 20 ). 

If now we apply Green’s theorem to the integral 
+ icMy ± icM^ ^ dydz 



+ {m, + icM^ ^ ± icMy dxdy^ (21), 

which is supposed to be taken over a closed surface, we find 
that it vanishes on account of the equations (2) provided t is 
supposed to be expressed in terms of Xy y, z according to some 
definite law which we shall take to be that expressed by (20). 

We now apply this theorem to a disc-shaped surface whose 
two faces very nearly coincide. We shall suppose that on one 
side of the disc the vector M represents the field of the advancing 
waves and that on the other side it represents the field obtain- 
ing just before the arrival of the waves. We shall also suppose 

* Cf. J. Larmor, Proc. London Math. Soc. (2), Vol. 1 (1903), p. 11 ; 
H. M. Macdonald, Electric Waves, p. 16 ; Proc. London Math. Soc. (2), Vol. 10 
(1911), p. 91. 

t The idea is practically due to Stokes, Math, and Phys. Papers, Vol. 2, 
p. 275, but was not worked out in detail. The different possible types of dis- 
continuity are discussed with some care by Love. The case in which E and H 
are continuous but some of their derivatives are discontinuous at the moving 
boundary may be discussed more simply by analysis analogous to that given in 
Hadamard’s Legons sur la Propagation des Ondes, Paris (1903), Ch. 2. See also 
Ex. 2, p. 23, and the references to Duhem and Silberstein on the next page. 
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that the derivatives of M are finite or behave in such a way 
that an application of Greens theorem is justifiable. Now 
let M be the discontinuity of M, i.e. the difference in the values 
of M at two neighbouring points on opposite sides of the wave- 
boundary, Then when the two faces of the disc coincide we 
find that a certain surface-integral over one face of the disc is 
zero. The surface-integral is of the same type as (21) except 
that M is written in place of JIf. Since the face of the disc can 
be chosen arbitrarily the integrand must vanish and so we 
obtain three equations of the type* 

M^ + icMy^^±icM,^~y==0 ( 22 ). 

These equations give 

Hence the wave-front advances with the velocity of light 
and the difference between the two electromagnetic fields at 
the wave-boundary behaves as a self-conjugate field in which 
Poynting’s vector is along the normal to the moving 
boundary. 

If the equation of the wave-boundary be expressed in the 
form 

F{x, y,-^,0 = 0, 

we find on calculating the values of 


^ ^ ^ 

dx* dy' dz' 



This is the differential equation of the characteristics, it 
expresses that the moving boundary moves normally to itself 
with the velocity of light. According to the theory of Stokes "I* 

* Equations equivalent to these are obtained in a different manner by 
0. Heaviside, Electrical Papers ^ Vol. 2, p. 405 ; A. E. H. Love, Proc, Londm 
Math. Soc. Ser. 2, Vol. 1 (1903), p. 37 ; L. Silberstein, Ann. d. Physik^ Vol. 26 
(1908), p. 751 ; P. Duhem, Comptes Rendus^ t. 131 (1900), p. 1171. 

• t Proc. Camh. Phil. Soc. (1896) ; Manchester Memoirs (1897) ; Math, and 
Phys^ Papers t Vol. 4. 
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and Wiechert*, Rontgen rays consist of pulses travelling through 
the aether, the energy in a pulse being confined within a thin 
shell. The above theory indicates that the front and rear 
surfaces of the shell move forward with the velocity of light. 

A slight modification of the preceding method can be used 
to find the conditions to be satisfied at a moving surface which 
is the boundary between two different media. If we write 
X = jB ± iX, N— E+ iH and assume that the surface-charge and 
surface-current can be neglected, the six boundary-conditions 
can be expressed by saying that the three quantities of type 

must be continuous as the boundary is crossed. The equation 
of the moving boundary is expressed as before in the form (20). 
When the moving boundary is the surface of a perfect conductor 
or perfect reflector, the boundary-conditions are simply that the 
three quantities of type 

should vanish f. 

In the last two cases the boundary-conditions do not imply 
that the boundary moves normally to itself with the velocity of 
light ; in fact, the motion of the boundary can be quite arbitrary. 


EXAMPLES. 


1. The surface of discontinuity is the sphere r^ct and the electro- 
magnetic field within this surface is expressed by the equations 


where 


-H, o). 




n as A r ^ e “ « <** ~ **) sin p (c^ ~ + «)• 


* Ahh, d, Physydkon, Ges. zu K'&mgsherg^ 1 Pr. (1896), p. 1 ; Ann. Phyi, 
Chem, Bd. 59 (1896), p. 288. See also J. J. Thomson, Phil. Mag. (5), Vol. 45 
(1898), p. 172. 

t These conditions are equivalent to the condition that the mechanical force, 
which would act on a charge moving with the normal velocity of the surface, 
must be along the normal to the surface. Of. Heaviside, Electrical Papergij 
Vol. 2, p. 514; Electromagnetic Theory ^ Vol. 1, p. 273; Hertz, Electric Waves ^ 
p. 257. 
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The constants 'p and 9 ' being known determine the constants A and c 
in order that the field outside the surface r—ct may be the electrostatic 


field for which the potential * is 


(Cambr. Math. Tripos, Part II, 1904.) 


2 . If the vectors E and H are known for all points (^, y, z, t) of a 
moving surface 

P(x,y, z,t) = 0 

the values of all the derivatives of E and and consequently values of E 
and H at points not on the moving surface can generally be found pro- 
vided F does not satisfy the difterential equation 



(Havelock, Proc. London Math, Soc. Ser. 2 , Vol. 2 , p. 297.) 

3. An electromagnetic field is conjugate to an electrostatic field. 
Prove that the flow of energy in the electromagnetic field takes place along 
the lines of electrig force in the electrostatic field. 


4. Let the line OC of length c be drawn in the direction of Poynting's 
vector at each space-time point 0 of a self-conjugate electromagnetic 
field and let OV represent a velocity v associated with the point 0. 
Prove that if this electromagnetic field is conjugate to another field 
{Ef H) in which cH is the vector product of v and Ey the direction of is 
parallel to VC, 


6 . Prove that when o-, €, are constants, the vectors Ey H of § 6 
satisfy the difierential equations 

^ _ 1 

^2 0^2 


„ ca d'^E a-u dE X , 
^ ei - -grad 




d^B 


crp 


dB 


c2 dt^ dt 


= 0 . 


(Maxwell.) 


Solutions of these equations for the case in which p =0 have been given by 

0. Heaviside, Phil. Mag.y Jan. (1889), p. 30 ; Electrical Papersy Vol. 2 , 
p. 478. 

H. Poincar^, Comptes Rendus (1893), p. 1030 j TMorie analytiqxie de 
la propagation de la chaleury Ch. 8 . 

J. Boussinesq, Comptes Rendus (1894), pp. 162 — 223 ; ThAorie 
analytiqne de la chaleury t. 2 (1903), p. 538. 

Kr. Birkeland, Archives des sciences physiquesy Geneva (1896), p. 5. 

0. Tedone, Rend. Linceiy Mar. 31st (1913), Jan. 18th (1914). 

M. I. Pupin, Trans. Amer. Math. Soc.y Vol. 1 (1900), p. 259. 


^ 6 . Prove that the function 


1 
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satisfies the equation 


dht 0% ^ 


7. If u satisfies the partial difierential equation 

fdu\^ /du\^ /0i^Y _ I 

\dx) ^\d^) Vdz) ’ 

and we make the transformation 


Xo-x+f gj, 

,du 


0< ’ 


where /is an arbitrary function, then u also satisfies the partial differential 
equation 

VaV WJ Wo) “ V^o/ * 

8. Plane electromagnetic waves fall on the convex surface of an infinite 
paraboloid of revolution whose surface is a perfect reflector. If 

the incident waves are given by expressions of the type 




0(a,3) 




where a^y’\-iz^ 3=a;+c^, the boundary -conditions at the surface of the 
paraboloid may be satisfied by subtracting from the primary field a 
secondary field represented by expressions of a similar type but with 
y — iz ^ 

as=a^ , 3=a~r-fc^. 

x+r 


If waves represented by the above expressions with a~y-iZy fi=^x-cty 
fall on the concave surface of the paraboloid the boundary-conditions at 
the surface of the paraboloid may be satisfied by supposing that the 
secondary disturbance is of the form - M' - M" where the fields 

M'\ if'" are represented by expressions of the above type, where a, 3 
have the values 

a=y+iz, ^ = 

// 2/ " ^iz ^ M n* 

a P --a-r-cty 


respectively. With these suppositions the forces are finite at the focus, 
when / is independent of a. 



CHAPTER II 


GENERAL SURVEY OF THE DIFFERENT METHODS OF 
SOLVING THE WAVE-EQUATION 

§ 10. The object to be attained. 

It has been shown in Chapter I that the solution of 
MaxwelFs equations can be made to depend upon the solution 
of a single partial differential equation which is either the 
wave-equation ft w = 0 or the equation -f k^u = 0 which 
is satisfied by wave-functions of the form u = y, z). 

The properties of functions satisfying these equations must 
accordingly be studied at some length. It is desirable, 
also, that all types of such functions should be studied 
and not merely those which admit readily of application to 
physical problems. If certain solutions of the fundamental 
ecpiations must be rejected in the treatment of the boundary 
problems of mathematical physics, a knowledge of their 
behaviour is at any rate useful as it gives a ^lear indication of 
the reason why such solutions must be rejected. There is, 
however, another reason why the scope of the inquiry should 
not be restricted. The theory of wave-functions forms a 
natural extension of the theory of functions of a complex 
variable* and may consequently lead to results of great value 
for the general theory of functions. 

* This point of view is adopted, for instance, by Volterra, Rend. Lincei (4), 
Iii2,pp. 225-330,274-287 (1887); iVi, pp. 107-115, 196-202 (1889); Vi,pp. 158-165, 
291-299, 599-611, 630-640 (1889); Rend. Palermo, Vol. 3, pp. 260-272. See 
also Appell, Acta Math. t. 4, p. 313 (1884) ; Painleve, Toulouse Ann., t. 2b 
(1888) ; BOcher, Bull. Amer. Math. Soc. Vol. 9, p. 455 (1903). The theory of 
functions of two complex variables is closely connected with the theory of wave- 
fi^ctions. Cf. H. Poincar6, Acta Math. Vol. 2 (1883) ; Vol. 22 (1898) ; H. F. 
Baker, Camh. Phil, Trans, Vol. 18 (1899), p. 431; Proc. London Math. Soc. 
Ser. 2, Vol. 1 (1903), p. 14. 
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We shall now describe briefly some of the principal methods 
of solving the wave-equation. 

§ 11. Reduction to ordinary differential equations. 

The aim of this method is to determine elementary solutions 
of the form 

^=/i(«)/.(^)/3(7)/4(S) (23) 

where /j,/2,/3,/4 are particular functions of their arguments and 
7, S are particular functions of x, y, z, t This method was 
used by D. Bernoulli in 1732 in the treatment of the vibrations 
of a hanging chain, the partial differential equation being how- 
ever in this case different. 

The general theory of elementary solutions is due to Lam6* 
who transformed Laplace^s equation into curvilinear coordinates. 
For a historical account of the development of the theory we 
may refer to Prof. Bocher s book Die Reihenentmckelungen 
der Potentialtheoriey Leipzig (1894) and to Byerly's Fourier 
Series and Spherical Harmonics. 

A simple elementary solution of the wave-equation is 
obtained by putting a = ^ = y, 7 = -gr, h — we can then 

take 

= gte+m2/ + n2:=t2>^ (24) 

where the constants 1, m, n, p satisfy the relation 

^ (P 4. m2 -f w") (25) 

and can be either real or complex quantities. 

When p is a purely imaginary quantity and Z,w, n unre- 
stricted, the solution is periodic and more general periodic 
solutions may be derived from this one by summation, 1, m, n 
being regarded as variable parameters subject to the relation 
(25). When 1, m, n are purely imaginary the solution (24) is 
appropriate for the representation of plane waves of mono- 
chromatic light, the intensity and phase of which are the same 
at all points of any plane perpendicular to the direction of 
propagation. 

* Liouville^s Journal^ t. 2, pp. 147-183 ; Legons tur let coordonnies 
etirvilignesy Fails (1859). See also E. Mathieu, Court de phytique matMmatique 
(1873). 
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Taking m = n = 0, i = — ^ = i;so that the axis of x is in the 
direction of propagation, we may write 

— — — — Ey — a cos v{x — ct\ 

Hz = a cos v(x — ct) (26) 

where a is a constant. These waves will be said to be linearly 
polarised in a direction parallel to the axis of y and will be 
called homogeneous because E and H do not depend on y and z. 
It will be noticed that the constant c represents the velocity of 
propagation of a phase of the disturbance. 

A wave-function of the type 

ft = sin vx cos vet 


is appropriate for the representation of standing waves. Ex- 
pressions for E and H may be written down by analogy with 
the above. To obtain a representation of plane waves in a 
conducting medium, we must use a solution of 


dhi 


+ = 0 , 


where 1(? has the complex value given in § 8. Putting F = 0, 
V = (0, 0, 

(27) 

where the real parts of the quantities are retained. If 

where rj is positive, the oscillations of the vector E are damped 
owing to the exponential factor 

The elementary electromagnetic fields that have just been 
found are fundamental in the theory of the reflection and 
refraction of light at a plane surface. This theory is given in 
the text-books on Physical Optics* and need not be reproduced 
here. Various attempts have been madef to prove that any 


* See for instance, Wood’s Physical Optics (1911), Ch. 13; Jeans, Elec- 
tricity and Magnetism (1911), Ch. 18. 

V See for instance a series of papers by G. Johnstone Stoney, Phil, Mag, 
o\, 43 (1897), pp. 139, 273, 368 ; Vol. 44, pp. 98, 206 ; BHt, Assn. Reports 
(1902), p. 539 ; Phil, Mag. Feb. 1903. The idea is probably due to Stokes. 
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electromagnetic disturbance in the aether can be represented as 
the sum of a finite or infinite number of elementary disturbances 
of the character of plane-waves travelling in various directions. 
Such a representation is generally only suitable within a 
restricted domain of the variables x, y, z,t\ nevertheless, it may 
sometimes be employed with advantage. 

When waves of all directions and frequencies are considered, 
the method of summation leads to Whittaker’s formula* 


rir r2ir 

I /[^ sin a cos )S 4- y sin asin /8 + -0 cos a — ct, a, )S] 
JoJo 

X sin a. dad/3 (28) 

for a wave-function. The case when 


= f <0< TT 

= 0 

has been used by Debye f in a discussion of the behaviour of 
waves of light in the vicinity of a focus. In order that an 
integral of the type (28) may represent a wave-function it is 
not necessary for the limits of integration to be those chosen. 
The limits for a may, for instance, be 0 and 0 where d is a 
root of an equation of type 

X sin 6 cos ^ + y sin 0 sin ^ -h z cos 0 — ct = F{0). 


In order to obtain other types of elementary solutions it is 
necessary to transform our differential equations to a system of 
orthogonal coordinates {a, v, w) for which the linear element is 
given by 

, du^ dv^ dw^ 

ds = -f;i + -y,+ -^, ( 29 ). 


If Hu, H^, Hyj are the three components of a vector H in 
directions normal to the surfaces u = const., v = const., w = const, 
through a point {x, y, z\ the corresponding components of rot H 
are of the typej 

'''•'[Uwhlif)] (“»• 

* Math, Ann, (1903). See also G. N. V^atson, Mess, of Math, Vol. 36 
(1906), p. 98. 

t Ann, d, Phys, Vol. 30 (1909), p. 735. ^ 

X See for instance, H. M. Macdonald, Electric Wavesj Ch. 6 ; M. Abraham, 
Math, Ann, Bd. 52, p, 81. Some very general transformation-formulae are 
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the new expression for div H is 

and the wave-equation becomes* 

n-irif s ly^ ^iv 

Law V du) ^dvKwu dv) ^ aw uv aw/_ 


It is also sometimes advantageous to transform the wave- 
equation to a system of coordinates for which 

dx^ + dy^ -1- dz'^ — c^dt^ = A^d^ -h B^dq^ -h C^d^ — D^dr^, 
the wave-equation then becomes 
d (BCDd4>\ d (Cl) A d(j>\ d /DABd(f>\ 


/BOBd<l>\ d /Cl) A d(l>\ d fDABd(j> 
■\ A d^j'^dvK B dv)'*"d^\ a dt 


= i f^BCdj^ 

dT\l) dr, 


§ 12. The generalisation of wave-functions. 

When a solution of the wave-equation has been obtained 
other solutions may be derived from it in various ways. For 
instance, the function obtained by ditterentiating the given 
wave-function any number of times with regard to the coordinates 
a?, y, z, ty is also a wave-function. By adding together arbitrary 
constant multiples of all the wave-functions obtained in this 
way we may obtain a very general type of wave-function. 

Another method of generalisation is to make an arbitrary 
change of rectangular axes. The wave-equation is a covariant 
for such a transformation and so is a wave-function of the new 
coordinates. A number of arbitrary constants can be introduced 
into the solution in this way. We can also make a linear 
transformation of coordinates for which the expression 

doiy -t- dy- -t- dz^ — c^dt'^ 

contained in the papers of V. Volterra, llend. lAncei, Ser. 4, Vol. 5, pp. 599, 
630 (1889), and J. Larmor, Cambr. Phil. Traus. Vol. 14 (1885), p. 121. 

Lam6, Journ. tie VEcole Poly technique y Cuh. 23 (1833), p. 216; Lemons 
gur leg coordonnees curvilitjneg, t. 2. A simplified proof was published by Lord 
Kelvin, Cambr. Math. Journ. Vol. 4 (1843). 
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remains unaltered in form and the preceding remarks still hold 
good. 

To illustrate this let us first of all add together two particular 
cases of Euler's wave-function i f{r + c<), viz. 

1 , 1 
r(r-c0 r{r^cty 


we then see that (r^ — is a wave-function. Generalising 
this by writing x — Xq^ y — ^ in place of x^ y, ty 

we obtain the wave-function 


1 

(x - x^y + {y- yof + zof toY 


,(34). 


When we have obtained a wave-function involving one or 
more arbitrary parameters we may obtain others from it by 
differentiating with regard to the parameters, or by integrating 
with regard to them after having multiplied the expression by 
an arbitrary function of the parameters. For instance, from 
the above wave-function we may derive the more general wave- 
function 


/ 


f(T) dT 

{x - x„y + (y- yoT + {z- ZoY -c\t- rf 


(36), 


where the integration is between constant limits. It is not, 
however, really necessary for the limits to be constant, we may 

for instance take them to be — oo and ^ - (r + Tq), where 

c 

^ 0 ^ = x^ 4- y^ + z^. The resulting integral is then a wave-function 
provided /(r) behaves in a suitable manner. If we take/(T) = 1, 
we obtain a wave-function 


u = 


1^ 

R 


log 


i2 — r — 7*0 
JS + r 4- n 


(36), 


where R^=={x- Xof + {y- yof 4- (-e^ - Zof. 


This function is independent of t and so must be a solution of 
Laplace's equation Au = 0. It is closely connected with the 
function used on p. 3 of Basset's Hydrodynamics, Vol. 2. 
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§ 13. Transformations. 

In addition to the linear transformations that have already 
been mentioned there are certain other transformations which 
enable us to pass from one wave-function to another*. 

The first transformation, which is analogous to inversion, is 
defined by the equations 


X = 


X 






t 


(37) 


where ^=zr^ ^ cH\ 

It is easy to verify that if f{x, y, Zy t) is a wave-function, 
then 

?/(|. S' 5. ?) w 

is also a wave-function f. 

The second transformation isj 
' — 0? • ^ y ' _ ~ ^ f' — ~ ^ 

^ z -- ct' ^ ^ z — cV ^ 2{z~-ct) ’ 2c{z—‘Ct) 

(39). 

It is easy to verify that if /(x, y, z, t) is a wave-function, 
then 

z — l^z — ct* z — cV 2{z — cty 2c{z^ct)\ ^ ^ 

is also a wave-function. Since is a wave-function we 

may deduce in this way that 

1 

-e z-ct 

z-^ct 

is a wave-function. 

It should be noticed that if we put 4 t = — c^, a^z + ct, 
a function of the type 

u = F(x, y, T)e-" (41) 


* For a general account of these transformations see a paper by the author, 
Proc, London Math, Soc, Ser. 2, Vol. 7 (1908). 

t This of course is a simple generalisation of Kelvin’s theorem for Laplace’s 
equation. The generalisation to the corresponding equation in n variables is 
mentioned by B6oher, Bull, of the Amer, Math, Soc, Vol. 9 (1903), p. 459. 

X Proc. Londm Math. Soc. Ser. 2, Vol. 7 (1909). This transformation is 
equivalent to a conformal transformation of a space of four dimensions which 
discovered by Cremona. Cf . Darboux, Legons sur les systhnes orthogonaux 
et lea coordonnies curvilignesy Paris (1910). 
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is a wave-function if F satisfies the differential equation 

fc"; 

The wave-function we have just obtained indicates that 
1 


d^F d^F _ 
da^ dy^ 


= F 


.(43) 


is a solution of the above equation. This solution is fundamental 
in the theory of the conduction of heat. It is evident that any 
solution of the equation of the conduction of heat in two 
dimensions can be used to construct a wave-function. 

Our second transformation theorem for the wave-equation 
also tells us that if F{a;, y, r) is a solution of the equation (42) 
the function 


-i) (44) 

T \T T T/ 


is also a solution. This result is due to J. Brill* and Appellf. 

The corresponding theorems for the two-dimensional wave- 
equation 

027 02 7 1 02 7 




C* 0^* 


are first that if /(oc, y, t) is a wave-function and 


the function 


52 — ^ ^ ^2 _ ^2^2^ 


.(45) 


!• ?) <*«) 


is also a wave-function. This is equivalent to Lord Kelvin’s 
theorem for Laplace’s equation if we simply replace ict by z. 

The second theorem is that if f{x, 3 /, t) is a wave-function, 
then 


2 ( 3 / -cO’ 2c{y-ci) 


(47) 


is also a wave-function. Writing 4t = 3 / — c^, o- = 3 / -f- as 
before we find that u = e^^F (x, r) is a wave-function if 


d^^dF 

0T 


,(48). 


* Messenger of Mathematics (1891). 

+ Liouville's Journal^ Ser. 4, t. 8 (1894). 
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The second theorem can now be used to show that 

^ ( 49 ) 

is a solution of this equation and that if F (x, t) is one solution 
the function 

T~^e-i^F(^, -i) ( 50 ) 


is also a solution. This theorem is likewise due to Brill and 
Appell, it can evidently be generalised to the equation in 
n variables. 


EXAMPLES. 


1. Prove that it is possible for a train of plane electsic waves 
to travel along an infinite isolated slab of dielectric material without 
being dissipated by spreading out into the adjacent empty space. Show 
that if 2a is the thickness and K the inductivity of the slab, the velocity 
of propagation of such waves of length X along the slab, when polarised so 
that their magnetic vector is parallel to it, is 




where 6 is the lowest real or the pure imaginary root of the equation 
27ra^ /A’'2 — AT 


(Larmor, Cambr. Math. Tripos, Part II, 1906.) 


2. Prove that with the notation of § 13, the function 

1 '^ + ^0 

{x — ~ yo)^ + (2 - ■“ Jssq -f xxo + yyo + 2-80 ~ ito 

is a wave-function, yo) k being arbitrary constants. 


3. If y, 2 , 0 is a solution of the wave-equation, the function 


Jo \T* r’ r’ t) 

is, under suitable conditions, a solution of the equation 

dt' 


dr 

2 


The quantity has the meaning assigned to it in § 13 and is supposed in 
the present case to be negative. 


B. 


3 
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4. The function 

9 [\/v{x-vt+ iy) -{-s/vix-vt -iy)] 

J _oo 

satisfies the equation 

^ d^_dV 

dx^ dt ’ 

of the conduction of heat in two dimensions : it is zero over a semi-infinite 
line which covers part of the axis of x and moves in the positive direction 
with uniform velocity v. The isothermal lines at any given instant are 
confocal parabolas. 


5. Prove that, under suitable limitations, the function 


is a solution of 


02 F 02 F 
-- 4- — - X2 F 

0;r2^0y2-^ 


Obtain in this way the particular solution 


6. Prove that if a;2+y2>^2^ the integral 


r- 

Jo 


satisfies 
and if :r2>^2^ the integral 


02 F 02F 02F d_V 
dx^ 0y2 “ 0^2 0^ » 


/: 


:^=F 


satisfies 


fjx^ -H 4a2 — 4at 

^_02F ^ 

dx^ ~ dfi dt ■ 


7. Prove that if po^='^o^+yo^ the integral 

Bink{^ + a)dS 

satisfies the equation 


^P+Po 

0»F 02F 
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§ 14. The elementary solutions. 

If we make Laplace’s transformation 

x — r sin 6 cos <^, y = r sin 6 sin <^, z cos 6 . . .(51), 
the equation Aw -f k^u = 0 becomes 


0^ 2 aw _i ^ 

r dr ^ sin 6 dO 


sin 6 


du\ 


a^w 


de) r^^m^ed4>^ 

This is satisfied by a function of the form 
u=^R(r)®(0)^((f>) 

a^^ 


+ A;“w=:0...(52). 


if 


defy 


q. ^ = 0 


.(53), 


ire id ^ fe) + [« ('»+!)- ® =o-(54). 

I dR r. /-rx 

'• dr 


^ 2 

dr"^ r 


The first equation is satisfied by ^ = cos {m<f> 4- a), the 
second by (cos 6) and (cos 6)y where these are the 
associated Legendre functions. The third equation may be 
written in Bessel’s form 


d^W 1 dw /trn\ 

+ + + 

where w — r^Ry and is satisfied by (kr) and (At). 

In these solutions m and n can have any constant values. 
It should be noticed that when n 4* ^ is an integer the Bessel 
functions that have just been written down are not independent 
anS the second solution ^ ^ (kr) of Bessel’s equation must be 
used. 


3—2 
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In dealing with a problem such as the effect of an obstacle 
on a train of electric waves, the secondary waves sent out from 
the obstacle must have the character of diverging waves at a 
great distance from the obstacle. In the case n—0 the differen- 
tial equation for R is satisfied by 

= (57), 

and if the real part of k has the same sign as od when the 
electric and magnetic forces are the real parts of expressions 
of the form we can obtain a solution appropriate for the 

representation of a diverging wave by taking the positive sign, 
for then we have a function of the form 

r 

Neither of the given solutions of Bessel's equation has the 
required form in fact 

Ji(kr) = \/ = ^cos kr. 

We may, however, obtain solutions of the form (57) by taking a 
suitable combination of the preceding solutions. 

In the case of electromagnetic fields in the free aether the 
physical interpretation of the elementary wave-functions when 
n is zero is as follows* : 

- cos ct) 
r 

- cos k(r ct) 

r 

- sin kr, coQ kct 
r 

To obtain the solution of (66) appropriate for divergent 
waves when n has any value we write f 

* A fuller discussion is given by A. Sommerfeld, Jahresberkht der deutsch. 
math. VereiUy Bd. 21 (1913). 

t The theory is due to Stokes, Phil. Tram. Vol. 158 (1868), p. 447; 
Collected Papers, Vol. 4, p. 321. See also Rayleigh’s Sound, Vol. 2, p..304. 
It should be mentioned that different notations are used by different writers. 
This is the notation used by Debye, Ann. d, Phys, Vol. 30 (1909), p. 57. 


Progressive divergent waves. 
Progressive convergent waves. 

Standing forced waves, source at origin. 
Standing free waves. 
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X. » - (-)"a/? t W “ - J'.+j W V' 

Vn («) = ■fn (®) - t’Xn (®) 

?n (*) = irn (x) + tXn («) 

These new functions rf„, are connected with Hankel’s 
cylindrical functions* by the relations 


rra: 

¥ 


-...(58). 


Vn 




When the real part of x is large and positive we have the 
asymptotic expansion 


2x 1 ! 

1 (n — 1) n (n + 1) (n -f 2) 
'(^2 2 l 


+ ... 


...(59). 


The series terminates when n is an integer and then gives 
a true representation of the function. To get ?n(^) we change 
the sign of i. Various other notations have been used for the 
solution of Bessel's equation that is suitable for the representa- 
tion of diverging waves. Lamb uses Dy{x), (i/= nq- J) to denote 
the solution of equation (56) which has the asymptotic value 

D,{x) = re-*’ (*+i) , 

while many other English writers use K^{ix) to denote the 
solution with the asymptotic form 

The following formulae will be found useful : 


X 1- 

• L 


{kry 


(kr)* 


2(2n-h3) ' 2.4.(2n-f-3)(2?i-f 5) 


-■] ' 


(60), 


• See Nielsen’s Handbuch der Cylinderfunktionen^ p. 16. 
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Ukr) = [l - 

L. (w-l)«(« + l)(w + 2 ) . ( 2 »)! ] 

{2krY 2! 

We have written down the last term in the series on the 
supposition that n is an integer. In this case 

i 4- {— | r | large . . .(62), 

(2n + 1) = (n + 1) f (^)- n {x) . . .(63), 


(2n + 1) (a?) = « (a;) + («)] 

Ky(x)=f e~®“®*'“coshi<a.da . 
Jo 


^ 2^r( ? )T ~ ( ' t<+T ) j 0 -K (»')>- i 


j -00 ’ 


-J (w=i)) 


^ ^0,p>0 (67). 


In the last formula n and p are supposed to be integers. 
For further properties of Bessel functions the reader should 
consult Gray and Mathews’ Treatise on Bessel Functions, 
Nielsen’s Handbuch der Cylinderfunktionen, and Whittaker’s 
Analysis. Tables are given in the first work and in Jahnke 
and Emde’s Funktionentafeln, Leipzig (Teubner). A few ad- 
ditions to the tables have been made recently by J. W. Nicholson, 
Proc. London Math. Soc. Ser. 2, Vol. 11, p. 104 ; Dinnik, Archiv 
der Mathematik und Physik (3), Bd. 20, Heft 3, 1912 ; J. R. Airey, 
Phil. Mag. Vol. 22 (1911), p. 85, Brit. Ass. Reports (1911); 
A. Lodge, Brit. Ass. Reports (1909) ; J. G. Isherwood, Manchester 
Memoirs (1904). 

The best definitions of the generalised Legendre functions 
for unrestricted values of m and n are those given by Hobson*. 


♦ Phil. Trans. A, Vol. 187 (1896), pp. 443 — 531. E. W. Barnes has recently 
given new definitions of the functions as integrals involving Gamma Functions 
which make it possible for the principal formulae to be proved very quickly. 
His definition of (x) differs from that of Hobson by a numerical factor which 
becomes rather troublesome when n is an integer and m is not. 
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We are interested here in the case when the variable is cos d 
and 0 is a real angle, such that 0 <0<7r. We may then put, 
with the usual notation of the Gamma and hypergeometric 
functions, 

P„“(cos0)= w + l; 1-m; 

( 68 ). 

Qn'^ (cos 0) = ^r— ; — ^ rTTT- C 

^ ^ ' 2sm(n -f m)7r r(l — m) 

X jcos (n 4 * m) 7 r cot^" | w, n -f 1 ; 1 — m ; sin^ ^ 


0 / 

— tan”' 2 ^ f + 1, 1 — 


COS'* 


.(69), 


PrT^ (cos 0) = 


cosec”'^ cos (n + ^)<f> ^ 

2^r(^)r(m <f io (2 cos <^ — 2 cos 0)^~'^ 

Jt(7n^i)>0 (70). 


When m is a positive integer, we have 

in fim 

~ 

Qn^ (^) = (1 — Qn (x) 

and when n is a positive integer ^ , 


•(71), 


Pn(^)=- 


1.3...(2n-l) 


1.2...n 


~ n(n — 1) 
I" 2(2n-l)* 


n(n-l)(n-2) (n- 3) _ 

2.4(2n-l)(2n-3) 

Pn (cos 0) = cos”^ — cos"~®^^ . sin’0 


■] 02), 


^ p Ceos 0) — ^ 7)” ^ /I 


0 


.(73), 

•(74), 


1 /^ / m (-1)" / I 1 ^ + A /^rx 

g„ (»,»)._ (_log_) (75). 


f»n+i 
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The last two formulae illustrate the method of deriving 
more complicated wave-functions from simple ones by differen- 
tiation. The functions 


1 ill r-\-z 
- and ^ loff 

7* VtV ^ V Z 

are in fact solutions of (52) when A; = 0. 

The formula 

J {x) {x) dx — 0 

_ 2 

2n 4- 1 ‘ (n — m) 1 


in which m, n, v are positive integers or zero, enables the 
coefficients in an expansion of a function in series of functions 
Pn^ (a?) to be determined by a simple integration. 

For further properties of the Legendre functions we must 
refer to Heine^s Kugelfunktionen, Byerly's Fourier Series and 
Spherical Harmonics, Whittaker s Analysis, Nielsen's Handhuch 
der Cylinderfunktionen and TheoHe des fonctions mitaspheriques, 
and to memoirs by E. W. Hobson* and E. W. Barnes f. 

Tables of the Legendre functions have been published by 
J. W. L. GlaisherJ, J. Perry§ and A. Lodge|| ; some tables of 
the functions P^ (p) have been given by H. TallquistlF. 

For the history of the functions of Legendre and Bessel 
the reader should ^consult the article by A. Wangerin in the 
Encyklopadie der Mathematischen Wissenschaften, Bd. ii. 1, 
Heft 5 (1904), p. 695. 


§ 15. Relations between various solutions. 

We have already remarked that when a wave-function or 
a solution of equation (52) involving arbitrary constants has 
been found, other solutions may be derived from it by the 
method of summation or integration. By choosing our sum or 
integral so that it represents certain simple solutions of the 

* Proc. London Math. Soc. Ser. 1, Vol. 22 (1891), p. 431, and op. cit. 

t Quarterly Journal, Vol. 39 (1908), p. 97. 

X Brit. Asa. Repcyrt (1879). 

§ Phil. Mag., Dec. (1891). See also Byerly, loc. cit. 

II Phil. Tram. A, Vol. 203 (1904). '' 

IT Acta Societatis Fennicae, Vols. 32, 33. 
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fundamental equation, a number of important identities may 
be obtained. A few formulae will be written down to illustrate 
this* 

If 4- — 2rri cos 6, 

I (2n + l)/„^j(^r)/„ + 4(ATOP„(cos0) 

R v(^^i) w=o ^ 

(Heine and Hobson) (77), 

4 = —“ 2 (2n + 1) {hr) (M) Pn (cos 6) (n > r) 

R ^’'^iw=o 

(Heine and Macdonald) (78), 


■*/„+j(At)P„(cos( 9) = ^ j 


We may illustrate the peculiar behaviour of certain definite 

integral solutions of our fundamental equation by the following 

example, in which k is supposed to be real and positive. 

rm 

Let I cosxtx(t)dt, 

Jo 

rm 

where ^ (0 ^ function such that I | % (0 1 convergent. 


then it may be proved by means of Fourier’s double integral 
theorem thatf 

“f (k>m) 

TT j «.Qo r — r*! 


Hence 


sin kR 


fm (^i) dr I 


is a solution of (52) which reduces to either fm {r) or fk (r) 
when 0 = 0, Solutions of (52) which are derived from the 
elementary solutions by integrating with regard to n have 
been employed by H. W. March J. He makes use of an inversion- 
formula § 


* Some very general formulae are given by L. Gegenbauer, Monatsh. /. Math, 
Bd. 10 (1899), p. 189. 

t This equation is obtained in a different manner by G. H. Hardy, Proc, 
London Math, Soc, (2), Vol, 7 (1909), p. 445. 

X Ann, d. Physiky Bd. 37 (1912). See also H. Poincar^, Comptes RenduSy 
t. 1^4 (1912), p. 795 ; W. v. Rybcynski, Ann, d, Phys, Bd. 41 (1913). 

§ This is in some respects analogous to the inversion-formula given by 
F. G. Mehler, Math, Ann, Bd. 18 (1881), p. 161. 
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where 


'Sa-4(C0S7) 


/( 0 ) = f P^_^(c 08 d)f(a)ada 

'S^o-4(c0S7)/(7)8in7d7 

V cos ttTT 2 /*" sin aj3. d/3 

') — ■— — — -| I ( ol ). 

-rros^ Wy {2(cos7-cos/3)}i 


It is sometimes instructive to find how a wave-function, 
depending on an arbitrary function, can be expressed in terms 
of elementary wave-functions. Now in the second example 
of § 5 the electric and magnetic forces are all of the form 


(82), 

r‘' 1_ ' ^ ^ 


or are the sums of terms of this form. Consequently, a function 
of this type may be expected to be a wave-function and it is 
easy to verify that this is the case*. We may now deduce that 




X ± iy^ 


is a solution of (52) and consequently it follows that tan"*^ ^ 

0 

and cot”^ 2 solutions of equation (54) when n = 0. We 
have in fact 

cos rriTT . cot’" ^ “* 2 ) ’ 

In the last formula m must not be zero or a negative integer. 


§ 16. The convergence of series of elementary solutions. 

When A? = 0 our fundamental equation (52) reduces to 
Laplace’s equation Au= 0 and we have the familiar elementary 
solutions 


* The theorem also follows immediately from a result given by A. R. 
Forsyth, Messenger of Mathematict (1898), p. 114, and E. W. Hobson, loc. eit. 
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r^Pn”* (cos 6) cos m{<^ — <^o). P*Qr!^ (cos 6) cos m((f)— ^„) 

^ P„™ (cos 0) cos m(<f>- <f>o), ^ Qn™ (cos 6) cos m (<^ - <f>„) 

(84). 

A pair of series of the type 

xQy„(0.<f,), 2(?)”^p„(0,<^) 

which converge when r=a are suitable for representing harmonic 
functions inside and outside the sphere r == a because the first 
converges absolutely when r < a and the second when r > a. 
The case in which A; 0 is very similar. When n is large 
the function {kr) may be replaced by 

1.3...(2n + l) 

and so a series of the form 


S 4 ^) 

n — 0 

converges like a power series. Again, when kr is real, we have 
I ?’»(A;r)| 2 = 1 + 1 «,(« + 1) 

+ |^(n-l)n(n+l)(« + 2).^, + ... 

2.4... (2n)" ^ ^ 

n being a positive integer. It is clear from this equation that 
I ?n (kr) I decreases as r increases, hence if a series of the form 


2 {kr)/n (0, 
n = 0 


converges absolutely for any value of r it converges absolutely 
for all greater values of r. 

For a discussion of the convergence of series of spherical 
harmonics we may refer to C. Neumann’s book Ueber die nach 
Kreis-y Kugel-und Cylinder- Funktionen foi'tschreitenden Entwicke- 
lungeuy Leipzig (1881); to Heine’s KiigelfunktioneUy Bd. 1, p. 435, 
Bd. 2, p. 361, and to papers by U. Dini, Ann, di Mat (2), t. 6 
(1874); H. Poincar^, Comptes RenduSy t. 118 (1894), p. 497 ; 
S. Chapman, Quarterly Joumaly Vol. 43 (1912), p. 1 ; T. H. 
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Gronwall, Math. Ann. Vol. 74 (1913), Vol. 75 (1914), Comptes 
Rendics (1914}), Amer. Trans. Jan. (1914), Vol. 15; C. Jordan, 
Gours d' Analyse, 2nd ed., Vol. 2, p. 252 ; and B. H. Camp, Bull, 
of the Amer. Math. Soc. Vol. 18 (1912), p. 236. The con- 
vergence of series of Legendre polynomials has been discussed 
very thoroughly by G. Darboux, Liouvilles Journal (2), t. 19 
(1874), p. 1 ; (3), t. 4, p. 393 ; O. Blumenthal, Dissertation, 
Gottingen (1898); E. W. Hobson, Proc. London Math. Soc. (2), 
Vol. 7 (1909); L. F6jer, Math. Ann. Bd. 67, p. 76; D. Jackson, 
Amer. Trans. Vol. 13 (1912). 

§ 17. The scattering of plane homogeneous electro- 
magnetic waves by a spherical obstacle. 

The effect of small particles in scattering incident radiation 
has been discussed very thoroughly by Lord Rayleigh * who has 
used it as the basis of a mathematical theory of the blue colour 
of the sky. The action of a single spherical particle is of 
fundamental importance and so the electromagnetic theory 
of the scattering of light by a dielectric sphere has been worked 
out by Lord Rayleigh f. Prof. Love;]: and other writers. This 
theory can also be developed so as to cover the mathematical 
theory of the rainbow. 

The more general theory of the scattering of incident radia- 
tion by a spherical obstacle § with arbitrary optical properties || 
admits of some v^ry interesting applications in the study of the 
colours exhibited by metal glasses, metallic films and colloidal 
solutions or suspensions of metals. The electromagnetic theory 
of these colours has been developed by J. Maxwell Garnett IT, 
G. Mie*^, R. Gansff and Happelff, who have considered 

* Phil. Mag. Vol. 41 (1871), pp. 107, 274, 447; Vol. 12 (1881), p. 81; 
Collected Papers, Vol. 1, pp. 87, 104, 518. 

t Phil. Mag. Vol. 44 (1897), pp. 28—52; Collected Papers, Vol. 4, p. 321; 
Proc. Roy. Soc. Vol. 84 (1910), p. 25 ; Vol. 90 (1914), p. 219. 

X Proc. London Math. Soc. Vol. 30 (1899), p. 308. 

§ The work of Stokes, Camb. Trans. Vol. 9 (1849), p. 1, with later appli- 
cations, Collected Papers, Vol. 4, and of L. Lorenz, Wied. Ann. Vol. 2 (1880), 
p. 70, opened up the subject. 

II The case of small conductivity was discussed by G. W. Walker, Quart. 
Joutn. Vol. 30 (1899), p. 204; Vol. 31 (1900), p. 36. 

^ Phil. Tram. A, Vol. 203 (1904), p. 385; Vol. 205 (1905), p. 237. 

** Ann. d. Phys. Vol. 25 (1908), p. 377. 

tt Ihid. Vol. 29 (1909), p. 280 ; Vol. 37 (19121. o. 881. 
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the cases of spheres and ellipsoids endowed with the optical 
constants €, /u., a. 

The particular case of a perfectly conducting sphere was 
worked out by J. J. Thomson* * * § and has been discussed in 
greater detail by J. W. Nicholson "t*. 

The problem is also of importance in connection with the 
theory of comets* tails which has been developed by Euler, 
FitzGerald J and Arrhenius §. The pressure of light on a 
perfectly conducting spherical obstacle has accordingly been 
calculated by K. Schwarzschildll and J. W. Nicholson H. The 
more general case of a sphere with the optical constants e, fiy a 
has been treated very fully by P. Debye** 

Let us assume that the electric and magnetic forces E\ H' 
at any point of space are the real parts of vectors E, H of the 
form Ae^\ where A is a complex quantity independent of t. 
We then write 

H ± ivE = where (86) 

flO) 

and we find that the differential equations satisfied by il/ in a 
medium whose optical constants are e, /i, a can be written in 
the form 

* • fin 0 W’’ ^ ' a? <’■*•>] 

* - ffL 'iw - 1 ^ H ' 

, , av(o , , €aco^ — iixaw 

where k = - — or = . 

c 

* Recent Researche^^ p. 437. 

t Proc. London Math. Soc. (2), Vol. 9 (1910), p. 67 ; Vol. 11 (1912), p. 277. 

X Scientific Writings, pp. 108, 531. 

§ Phys. Zeitschr. Vol. 2 (1901), pp. 81 — 97 ; Das Werden der Welten, Leipzig 
(1907), p. 85. 

II Sitzungsber. d. Kgl. Bayer. Akad. d. Wiss. Vol. 31 (1901), p. 293. 

If British Association Reports (1910), p. 544; Monthly Notices of the Royal 
Astronomical Society, Vol. 70, p. 544. See also J. Proudman, Ibid, Vol. 73 (1913), 
p. 535. 

** Ann. d. Physik, Vol. 30 (1909), p. 57. 
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These equations may be satisfied by putting 
-Mr = ^ (rfl) + A^rXi 

^ r drdO r sin 0 d<f> 

1 3nrn) ^A; 3(rn) 

^ rsin^ 9r3</> 9^ 

where the function ft = i7 ± iF is a solution of (52). 

The electric and magnetic forces may be derived at once 
from these expressions by equating the ambiguous and un- 
ambiguous parts as explained in § 2. 

We shall now assume that the incident wave of plane 
homogeneous monochromatic polarised light is represented by 

(ifr, M^) = ei^rcoBe±ii^ ^ (39^^ 

the electric vector being parallel to the axis of y. 

The corresponding function fto is obtained by solving the 
equation 

02 

(rfto) + A;Vfto = sin 0 . 

We easily find that 

rfto = [cosec 6 . -f {9, <f>) H- /a (6, (f>) 

, (90). 

Choosing the unknown functions so that rfto is finite for 0 = 0 
and ^ = TT we obtain finally 

r 0 0 

rfto = 2 cosec 9 , ^Her _ ^-ikr 

JL M 

(91). 

We may now assume an expansion for rfto of the form 
— 2 an'^n{kr)P^{G0^9). 

n-X 

To determine the coefficients an we multiply by sin 0 and 
differentiate both sides of the equation. Then since 

~ [sin 9 . (cos 0)] = - 71 (n 4- 1) sin 9 . Pn (cos 9) (92), 
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the coeflScients may be determined at once with the aid of 
Lord Rayleigh s expansion* 

ikr sin 6 . = 2 (2n 4- 1) y/rn (kr) Pn (cos 0) sin 0 

w=0 

(93). 

We thus obtain 

rOo - 2 Pn (cos 0) (94). 

n=l 7? (n 4- 1) 


Now let Hi, fla be the functions from which the electric 
and magnetic forces in the scattered light and transmitted 
light may be derived respectively. The appropriate forms are 
given by equations of the type 


rUi== 2 fn (kr) Pn^ (cos 0) cos <t> 

n=l 

r Fi = 2 JSnfn (f(^r) Pn^ (cos 0) sin ^ 
»=1 



r(7a= 2 C nyjrn (hr) Pn^ (cos 0) cos <}> 


r Fa = ^ Dn'^n (hr) Pn^ (cos 0) sin <l> 

n = l 

where h, tj are the values of k, v respectively within the sphere. 

It is easy to deduce from (88) that the tangential components 
of the electric and magnetic forces are continuous in crossing 
the sphere r = a, if when r = a 

h(J, + V,) = hV, 

* Theory of Sound, Vol. 2, p. 272. The expansion was also obtained 
independently by Heine, Kugelfunktimen (1878), Bd. 1, p. 82. 
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These conditions give 

- An^n (ka) - - On-^n (ha) = ^ (ktt) 

1 4- 1 

(ha) - (Aa) = ^ (*«) 

kB„ r„ (*a) - Ai)„ (ha) = J (*a) 


...( 98 ). 


Solving these we get 

^ — {n+i 1^-2 ^ n --.2L ^n+i ^ 

ji(rH- !)(?„’ “ n(?i + l)G„’ 

O - ^'M+i ^-2 2ri, + 1 En jy _Jl .-n+i + 1 J'n 

where 

Hn = (A;a) -v/r/ (Aa) - (A;a) y^n (ho)' 

En = {ha) yjrn (ha) - (ka) yjrn (ha) 

En = U (ha) yjrn (ka) - fn (A:a) -i/r/ (i;a) [ , . .(99). 

On ~ (ha) ^n (ko) v^n (ka) yjrji (ha) 

Ln ~ v'^n (ha) ^n (ka) ^^n (ka) (ha) 

We can proven that our series all converge absolutely and 
uniformly at about the same rate as a power series of the form 


" 1 . 3 . ..( 271 + 1 )’ 
kct^ 

where x is either kr or — , The proof depends on the fact 
that when n is large we have approximately 


(x) 


r^+l 


i?s ... (271 + 1 ) ’ w - 1 3 + 1 ) ^ 


?n (^) 1 . 3 . . . (271 1) « ^ 

U/ 


^n (x) = — 1 .3 (2n — 1 ) nix‘"^~^ e~^^. 
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§ 18. Free damped vibrations for the space outside the 
sphere. 

It should be noticed that some of the terms of our series 
become infinite when either Gn^O or Ln = 0: fortunately, 
however, the roots of these equations turn out to be complex 
and so when k is real no values of k need be excluded from the 
discussion. The damped vibrations determined by the roots of 
the equations Ln=0 may be distinguished as the electric 
vibrations, those determined by equations of type (zn = 0 as the 
magnetic vibrations. Some of the roots of the equations have 
been calculated for the case of a totally reflecting sphere by 
Sir J. J. Thomson*, who finds that the roots are all complex. 
The vibrations for the space inside a totally reflecting sphere 
have been discussed by Prof. J. W. Nicholson f, those for the 
space between two concentric spheres by Sir J. J.^ThomsonJ, 
Sir Joseph Larmor§, Prof. H. M. Macdonald || and A. LampalT. 

P. Debye, who has calculated some of the roots for a case of 
a dielectric sphere, finds that the roots are complex and of two 
types. When the index of refraction is large, the imaginary 
part of a root of the first type varies very little with the index 
of refraction N and approaches a limit different from zero when 
N <X) . If on the other hand p is a root of the second type, 
Np tends to a finite real limit, viz. a root of (Ap) = 0, when 
A--> 00 and so the imaginary part of a root of this type must 
be very small when A is large. 

The vibrations belonging to the space outside a sphere 
must be in all cases damped on account of the loss of energy by 
radiation; when the refractive indices of the outside medium 
and sphere are very nearly equal, they are clearly very strongly 
damped ; thus it is only when the refractive index is large that 
some of them are durable. It is doubtful whether a substance 
exists which has a large refractive index and does not absorb 
light to a marked extent. 

* Proc. London Math. Soc. Ser. 1, Vol. 16 (1884), p. 197 ; Recent Researches, 
p. 361. 

t Phil. Mag. 1906, p. 703. 

X ^Recent Researches, p. 373. 

§ Proc. London Math. Soc, Ser. 1, Vol. 26 (1894), p. 119. 

il Electric Waves, Chapters 6-7. IT Wien. Ber. 112 (1903), p. 37. 

B. 4 
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It will be seen later that the characteristic vibrations play 
an important part in determining the size of the sphere on 
which the pressure of a given type of incident radiation has 
a maximum value. 

Prof. Love* has used the solutions corresponding to the 
characteristic vibrations to discuss the mode of decay of an 
arbitrary initial disturbance. He makes use, in fact, of the 
functions where k is one of the roots of one of the 

equations (?n = 0, = 0 ; only the sphere is treated as a perfect 

conductor. 

This method can easily be extended so as to provide us 
with a method of discussing the problem of the scattering of an 
arbitrary primary disturbance by a spherical obstacle. In this 
method we assume that the total disturbance outside the sphere 
can be represented by 

r = 2 2 2^n.m,p^n {K'd (cOS 0) COS m {(f> - (f>o) ' 

ftaeO «f|s=0 jp 

rF = 2 2 2 ^|r„ (kp'r) P„’”(cos 0) sin m {<j) - 

(®). 

where the kp8 are roots of one of the equations of the type 
Tfyjrn (ha) (ka) - (ka) (ha) = 0 ) 

(ha) yjrn (ka) ~ rjylrn (ka) yjrn (ha) = 0) “ 

The coefficients must then be chosen so that this total disturb- 
ance has the same character as the primary disturbance at its 
singularities outside the sphere and at an infinite distance, 
taking into account of course the presence of diverging waves 
from the spherical obstacle the effect of which is, however, 
negligible at infinity. 

The field inside the sphere is represented by equations 
similar to (o)), one with hp written in place of kp. The boundary 
conditions are satisfied in virtue of (100). 

If a series of type (co) should fail to represent the disturbance 
outside the sphere, it may be necessary to add terms cor- 
responding to the free characteristic vibrations. These are of 
the form (co) with the function fn written in place of yjrn and 
numbers kn determined by equations of type Ln = 0, Gn == 0. 

* Proe, London Math, Soc, Ser. 2, Vol. 2 (1904), p. 88. 
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§ 19. The case of a very small obstacle. 

When a, the radius of the sphere, is very small compared 
with the wave-length \ of the incident radiation, we may treat 
ka and ha as small quantities. We may then obtain some idea 
of the relative magnitudes of the different coefficients in our 
series by using the expansions (60) and (61). 

It is easy to see that the values of Am Bn decrease very 
rapidly in absolute magnitude as n increases. The disturbance 
radiated from the sphere can consequently be represented 
approximately by superposing a small number of partial waves, 
the effect of the others being negligible. 

Remembering that vh = 7)ky we find that when H is finite 


Hn 

En 


(P - A-) 

(2n + 3).P.3^..(2^T)2 ’ 
(n 4- 1) k^h^^~^7) {k-^ - li^) 

I2.32...(2n4-1)2 

2n4-l 



and it is easy to see that all our series converge. 

It appears from these expressions that the w-th magnetic 
wave, i.e. the disturbance due to the nth term in the expansion 
for [7i, is of the same order of magnitude as the (n -f l)th 
electric wave, i.e. the disturbance due to the (n -f l)th term in 
the expansion of Fj. This is in sharp contrast with the result 
obtained by Sir J. J. Thomson for the case of the totally 
reflecting sphere wherein the nth electric wave and the nth 
magnetic wave are of the same order of magnitude. 

The first electric wave is clearly of chief importance and 
Mie has proposed to call this Rayleigh's radiation. We easily 
find that 


' ~ 2k^ + h^ “ ’ 2 . 3" 3k^ + 2/t» “ 


ik 




.( 102 ). 


The following diagrams, which are taken from Mies paper, 
indicate the character of the electric lines of force for the first 
four partial vibrations of each type. For the magnetic waves 


d— 9 
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= 0 and so the electric lines of force are spherical curves. 
In the case of the electric waves the lines of force lie on 
certain cones and the diagrams represent the intersections of 
a sphere with these cones, the vertices of the cones being at 
the centre of the sphere^. 


§ 20. Polarisation of the scattered light. 

Let us now look for cases when the light scattered by the 
sphere is linearly polarised. It is easy to see that and 


both vanish when = 0 and = 0. 

ou do 


These conditions are 


both satisfied by </> = 0, i.e. when the observer looks in a direction 
at right angles to the electric vibration in the incident wave 
(Fig. 3). 

It appears from the figure that the component of the. electric 
vector of the scattered light, which / 

is at right angles to the direction / 

in which the observer is looking, / L \ 

is parallel to the electric vector / / \ 

in the incident wave. In a similar / / \ 

way it is found that and Mq \ 

both vanish when <^ = ± ^ , i.e. / ' — 

when the observer looks in a \j^ 

direction at right angles to the H 
magnetic vibration in the incident 

wave. The magnetic vibrations in the incident and scattered 
waves are now found to be parallel. 

The experiments of Steubingf with different kinds of col- 
loidal gold solutions have shown that when the solution is 
illuminated with polarised light and viewed in the manner 
described, there is always a small quantity of unpolarised light 
sent out from the particles, but the greater portion of the 
scattered light is polarised in the way the theory requires. 
The slight disagreement between the theory and observations 
is attributed to the fact that the metallic particles are probably 


* author is indebted to the publisher of the Annalmder P/iysifc, Herr Johann 

Ambrosius Barth, for permission to reproduce the figures on pp. 62, 69 and 64, 
t Digsertation, Greifswald (1908) ; Ann, d, Phys, Vol. 26 (1908), p. 329. 
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not all sph^tes*, that some may have developed into crystals 
perhaps of octahedral form. The mathematical theory of the 
seattering of waves has not yet been fully developed. The 
problem is, however, One of great importance in meteorological 
optics. The case of a regular distribution of atoms or molecules 
has recently been brought into prominence f by experimental 
work on the scattering of Rontgen rays by a crystal J. Approxi- 
mate mathematical theories have been given by several writers §. 


§ 21. Intensity of the scattered light. 

When Rayleigh’s radiation alone is considered, we have 


H.-O, 

pr — 9 {rV) 

^ riun^ d<f> ’ 

ifid(rV) 

r d0 ’ 

where 


* vr drdd 

E 1 d’(rV ) 

* w sin 0 drd<f> 


.(103), 


sin ^ sin ^ . .(104). 


At a great distance from the origin the radial electric force 
is of order ^ while the transverse electric and magnetic forces 

are of order ‘ Hence the intensity of the scattered light 

diminishes ultimately according to the inverse square law when 
points on the same radius are considered. It also varies as the 
square of the volume of the particle. 


* This remark is made by both Maxwell Garnett and Mie. 
t It had previously been considered Lord Rayleigh, “ On the influence 6f 
obstacles arranged in rectangular order on the properties of a medium,” PhiL 
Mag. (5), Vol. 34 (1892), p. 481; Scimtific Papen, Vol. 3, p. 19; and by 
T. H. Havelock, Proc. Roy. Soc. A, Vol. 77 (1906), p. 170. 

t Lane, Friedrich, und Knipping, Sitzungsher. der Kbnigl. Bayerinchen 
Akad. d, Witi, June 1912. 

§ See, for instance, W. L. Bragg, Proc. Camh. Phil. Soc. Vol. 17 (1913), 
p. 43 ; Proc. Roy. Soc. A, Vol. 88 (1913), p. 428 ; M. Laue, MUnchener Bet 
(1912), p. 863; Ann. d. Phys. Bd. 41 (1913), p. 989, Bd. 42 (1913), p. 397; 
P. P. Bwald, Phyt. Zeittchr. (1913), p. 466 ; L. S. Omstein, Amtterdam Proc, 
(1913) ; M. Bom u. T. v. Earr^n, Phyi. ZtiUchr. (1912), p. 297. 
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An approximate formula for the intensity is 

I == l(^a^ (cos^ 6 sin® 4- cos® <^) . . .(105), 

the intensity of the incident radiation being 

It a = 0 for the medium outside the sphere, the quantity 
k is inversely proportional to the wave-length X of the incident 
radiation. Hence when h is large compared with Jk, we have 
Lord Rayleigh s result that the intensity of the scattered light 
varies inversely as the fourth power of the wave-length. The 
short waves are on this account scattered far more profusely 
than the long ones and so we have an explanation of the blue 
colour of the sky. 

The above formula for the intensity of Rayleigh’s radiation 
indicates that there is no light of this type in a direction for 

TT 

which ^ , i.e. when the observer is looking in, a direction 

parallel to the direction of the electric vibration ii/the incident 
wave. To obtain an expression for the intensity of the light 
sent out in this direction we must take into account the second 
electric wave and the first magnetic wave. Referring back to 
the expressions for B 2 and Ai, we find that when a is neglected 
both inside and outside the sphere, the intensity of the scattered 
radiation varies inversely as the eighth power of the wave- 
length*. This corresponds to Tyndall’s “residual blue” which 
is purer than the blue seen under other concfitions. 

§ 22. The absorption of light by a spherical obstacle. 

The total energy absorbed from the incident radiation by 
a single particle consists of two parts; first of all the energy 
scattered and secondly the energy which flows into the particle 
and is transformed into heat or chemical energy. Both these 
quantities may be calculated by the following method which is 
a simplification of the one given originally by Mie. 

The flow of energy across unit area of a very large sphere 
concentric with the spherical particle takes place at a rate 
measured by the radial component of Poynting’s vector, i.e. 

* Lord Rayleigh, PhiL Mag. Vol. 12 (1881), p. 81. 
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Now if Ey E and Hy H are conjugate complex quantities, 
this component is equal to 

\ [{E^e^^ + EqB-^^) 

- + E^e-^^) -f 

Integrating with regard to t so as to obtain the mean value 
of this quantity over a period, we obtain an expression which 
may be written in the formf 

S = ^ - MeM^ - 

olv 

where My M* are the values of M corresponding to the signs 
+, — respectively in (86) and My M* are derived from them by 
changing the sign of i. 

The function fl for the outer space is the sum of CIq and Hi, 
hence we may write 

rfl =-• S Pn^ (cos 0) [une^^ Vne^ ’h tVne"^] \ 


rll = 2 Pn^ (cos 0) -i- 4* Wne^] 

»=1 I 

00 

rn*= 2 Pn^(cOS0)[Une~^ -i- Vne~^ Wne^*^] 
n=l 


...(106), 


where 


rfl* = 2 Pn^ (cos 6) [Une^ + Vne~^ + Wne'^ 

n=l 

1 271 + 1 . ,, , 1/. . 


Un = ri'(^r+T) Vn = ^{An-\- Bn) fn (kv), 

Wn = 2 {An Bfi) ^n {kv ). 

When the expression S is integrated over the spherical 
surface, the result can be expressed in the form 

-^0 4 " -^1 ^Qly 

where Iq depends only on the incident radiation, Ii depends 
only on the scattered radiation and represents the amount of 
energy absorbed by scattering, 7oi depends on both types of 
radiation and represents the total absorption of energy from 


t We assume now that <r=0 outside the sphere. 
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the incident wave. The sum of the three terms with its sign 
changed represents the amount of energy that flows into the 
sphere and is truly absorbed by it. 

Performing the integration with regard to <f> and collecting 
the different terms together, we find that the surface integral 
can be written in the form 






X ^ dw^ dw,n 


\\dr ^ dr j \dr dr 
+ 


dr dr 


kTc (Un 4- Vn) (lim 4- Vm) j • 

This expression can be simplified with the aid of the relations 






Pm^4- sin^ 6 


dPr^ dP^y^\ dd 


), 


dd dd J sin d 

2?i4- 1 


= 0 mi=^n 


711 = n 


...(108). 


The firat of these is evident since the integrand is the 
differential coefficient of a function which vanishes at both 
limits. To prove the second we make use of the formulae 

sin d “ Pn‘ (cos ^) = — cos ^ . Pn' (cos d) 


— ?i (m 4- 1) sin d . Pn (cos d) 


P„1 (cos ^) = - Pn (cos d) 
+ n (n + 1 ) cos ^ 


d 


+ m (ni + 1) cos 0 . PmPn 


= -^*^[cos0.P„>P,„^] 


...( 109 ). 
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The integral is thus equivalent to 

n (ti 4- 1) m (m 4- 1) f Pn (cos 6) P^ (cos 6) sin 6d6, 

Jo 

and so has the value we have assigned to it*. 

Our surface integral now reduces to 

Since cr = 0 outside the sphere, we have ]c = k. Also when 
r is very large we may use the approximations 

2^ + 1 


\kr — n + 1 . - v 

( 2) n(n + l) 


Uft cos 


*>•-« + 


i (An - Bn) (- i)»+> 
We thus find that /o = 0 and 


L..(110). 


/. = n (n + 1) [i** (An + Bn) - (- i)» (^n + iin)] • • .(HI), 

T —1!^ 5 n^jn + 1 )° ^ ^ Z j. W fi ^ /'119\ 

V ,„i 2n + l 

When 1 / = 1 these expressions are just half those given by 
Mie on p. 436 of his memoir. The reason why they must be 
doubled is that our expressions for the incident light give an 
intensity J. 

In the case of a solution containing N spherical particles 
per unit volume the total absorption coefficient is thus 

A=^ t n{n + l)\i^{An + Bn)-{-i)’'(An + Bn)]...{m). 


IV 


»=1 


It is interesting to study the variation of this quantity with 
the wave-length for particles of different sizes. Mie has drawn 
* Cf. L. Lorenz, Omvte* scimtifi^uety p. 626. 
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the absorption curves for particles of gold varying in size from 
20 /A/i to 180 fifjL, 



uao ^00 650 600 C60 


Fig. 4. 

It will be seen that for ruby red gold solutions containing 
very fine particles there is an absorption maximum in the 
green corresponding to a wave-length of about 525 /i/i. 

Mie has also drawn curves showing the pure absorption in 
colloidal gold solutions. 

The colours of silver particles in colloidal silver solutions 
have been discussed with the aid of the mathematical theory 
by E. Muller ♦. The particles of a silver solution show beautiful 
colour phenomena, all colours of the spectrum from the extreme 


Ann. d. Phyi. Bd. 36 (1911). 
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blue to the extreme red being present*. For other applications 
of the theory we must refer to the papers of Lord Rayleigh, 
Garnett and Mie and to Prof. Wood’s Physical Optics, 

The mathematical theory for a large number of particles 
has been developed further on certain simplifying assumptions 
by F. Hasenohrlf, A. SchusterJ, W. H. Jackson §, L. V. King||, 
A. Einstein IF and M. v. Smoluchowski**. It has been used 
recently by W. J. Humphreys ff in a study of the effect on 
climate of large quantities of volcanic dust in the upper 
atmosphere. 

§ 23. The pressure of radiation on a spherical obstacle. 

We shall now calculate the pressure of radiation on a spherical 
obstacle, following the work of Debye except in some of the 
details. The pressure is calculated on the assumption that the 
force exerted by an electromagnetic field on a unit charge 
moving with velocity v has components of typelJ 

+ (114), 

c c 

and that the equations of the field can be derived by a process 
of averaging from the electron equations 

^ ^ I fdE \ _ IdH] 

( 116 ). 

divi?==^, div.ff=0 J 

Now if 

= h - Ey^ - + h 

Xy=E^Ey+ H^Hy= F* 

X, = E^E, + H^H, = Z^ [■•••( 116 ), 

S^^c(^EyH,-E\Hy) I 

* A coloured reproduction of an ultramicroscopic picture of a silver solution 
is given by H. Siedentopf, Ber. d. Deutsch, Phys. Ges. (1910), p. 6. 

t Wien, Berichte (1902). J Astrophysical Journal, Vol. 21 (1908). 

§ Bull, of the Amer, Math. Soc. Vol. 16 (1910), p. 473. 

11 Phil. Trans. A, Vol. 212 (1913), p. 375. 

IF Ann. d. Phys. Bd. 38 (1910), p. 1275. 

** Boltzmann Festschrift (1904), p. 626; Ann. d. Phys. Bd. 25 (1908), p. 205 ; 
Phil. Mag. Vol. 23 (1912). 

ft Bull, of the Mount Weather Observatory (1913) ; Journal of the Franklin 
Institute, Aug. (1913). 

it We assume now that <r=0, e=/i=l for the external medium, so that *^ = 1. 
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.(117). 


etc., we have identically on account of (115) 

„ ^dx^dXy dX, idS^ 

^ * dx ^ ~dy dz c^'dt *** 

Hence by Green’s theorem, if (I, m, n) are the direction 
cosines of the outward drawn normal to a surface a enclosing 
all the charges in the field, 

jj(lXx-hmXy’h nXg)d(r =Jjj dy dz. ..(IIS). 

The a;-component of the force exerted by the electromagnetic 
field on the obst^xcle is thus the same as if there were tractions 
X, F, Z at each point of the surface and a volume force 

other components. 

Now when we integrate with regard to t so as to obtain the 

dt 


mean value of the pressure over a period, the term 


contributes nothing on account of the periodicity of S. The 
pressure may consequently be calculated from Maxwell’s tractions 
X, F, Z. In the present case the mean value of the pressure 
in the direction of the axis of z is derived from the quantity* 
(IZa: + mZy + nZz) = — ^ cos 0 [Eq^ -f E^'^ + Hq '^ 4- 

- e;^ - h;^] - sin 0 [e;e,' + h;h^'i 

The surface integral of the mean value of this quantity over 
a period is accordingly 


("n r2iT _ _ _ 

— I ^Q>o^0 {^E^EQ’\-E^E^'\rH^IlQ-\-H^II^ — EyEf.’—HyHr) 

^JqJo _ __ _ 

4* sin 0 (EyijQ 4“ EyE^ 4- HyHQ 4" HyH^^ sin 0 . d0 d<f}. 

When the incident radiation only is taken into account, this 
integral becomes simply 

^ fir f2ir 

— - I I cos ^ sin 0 d0d<f> 

^ J 0 J 0 

and vanishes completely. Again, it is easy to see that when r 
is large the components Ey, Hy for the scattered field are small 

^ We assume now that the surface <r is a sphere whose centre is at the origin. 
Our expression is easily obtained by writing down the expressions for j&y, Ey\ E, 
in terms of Ey\ Ee\ E^\ 
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compared with the transverse components, a product such as 
HtHq being of order may be neglected when r is large, and 
so we have only to consider the integral 

^ rrr /•2ir _ — 

— r I I d6dj>y 

0 J 0 

where E = Eo-\- Eiy The terms which depend 

only on the incident field may, moreover, be disregarded. 

To evaluate we need the values of the integrals 

in- 

By using the relations (109) we may transform the first of 
these into 

(n + 1) m (m -f 1) f Pn (cos 6) Pm (cos 0) sin 6 cos 0 d0 
Jo 

~ iodd ^ ^ 

Now 

(2n + 1) cos 6 . P„ (cos = (w + 1) P„+, (cos 0) + nPn-i (cos 0) 

(119). 

and 


(2n + 1) cos 0 . P„> (cos 6) - nP\+t (cos 0) + (rj + 1) P>„_, (cos 0) 

( 120 ); 

hence when the second integral is integrated by parts we 
obtain two integrals of the type (7G) and so we have finally 


^(ii-l)’n(n + l)’ 
(2« - 1) (2n + 1) 
^ (n + 1) (n + 2^ 
(2fl+l)(2a + 3) 


m gfe n ± 1 V 
m = n — 1 

+ 1 , 


( 121 ). 


When the second integral I 2 is integrated by parts 




}PmJ sin 0d0 = 0 


^ n(n+l) 

2n+l 


m^n 




it becomes 

...( 122 ). 
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Writing P, in the form 

P, = - ^ /o"C 


+ sin 0 cos dd6d<f> 

and making use of equations (121), (122) and (106), we obtain 


p n^{n + l){n + 2f [ di;„+A /dw„ dv„\ 

^ “ »=i (2« + 1) (2« + 3) L\ dr dr 1^) 

dxVfi^i d^Wf^ 1 ? o / \ ✓ - — V 1 

+ ~d^^^ (“n+i+ «'«+i)(Mn+Wn) + A^n+i W„ 

/ dUn / ^„+i d>ln+i\ dWn dw„+i 

\ dr dr j V dr dr ) ^ dr dr 


+ k^(Un+ v„) (m„+j + V„+,) + k^WnWn+j 


. “ n(n+ 1) r. x/dM„ , dr„\ 

-SThT L<“" ^ '■> (t + s) 


_ ,/dUn , dv„\ , _ dWn d-Wnl ^ 


We now use the asymptotic expressions for Un, Vn, Wn when 
r is large and omit the terms that depend only on the w’s. 

We also write 




2n+l 




and obtain after some simplifications 
P2=— 2 (2m + l)[an + + 

^ n=l 


^•n-l 


2n + 1 
71 (77 4- i) 




OO r 1 _ 

00 ^ ^ 2) _ 

+ -n-l^' 2 -A-A [«««»+! + an«n+i 
n=l n -f i 


"1" ^n^n+1 "h ^n^n+i] 

(123). 


This expression turns out to be negative, i.e. the pressure 
acts in the direction in which the incident light is moving. If 
the constants in the incident light are chosen so that its 
intensity is unity, the above expression must be doubled. The 
numerical value of the pressure has been calculated from the 
abo^ formula by Debye in a number of cases. When the 
radius of the spherical obstacle is small compared with the 
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wave-length of the incident light, the functions yltn (Jca), 5'n (^^)> 
Sn(^Gt) occurring in the expressions for a„, can be 
expanded in ascending powers of a. This method, however, 
fails when ka approaches unity and numerical values of the 
functions must be used. 

In the case of a totally reflecting sphere 


a = 


(ka) 
U(lca)^ 
and Debye finds that if 


= ak — 


/3 = 


27ra 


^n(ka ) 


.(124), 


where X is the wave-length of the incident light, L denotes the 
light-pressure and W = is the energy of the incident train 
of waves per unit length of a cylinder circumscribing the 
spherical obstacle and having its axis parallel to the direction 
of motion of the waves ; then 

^ = ( 125 ). 


The first term of the series was given by Schwarzschild. 
The convergence of the series is slow as p approaches unity and 
several terms of the series (123) must be taken into account. 
By using numerical values of the functions ^n(p), S’n(p) and 
their derivatives for ?i = l, 2... 5, Debye has succeeded in 

drawing a curve ^or 
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It will be seen that the pressure has a maximum value for 
a certain value of p, approximately equal to 1. When p is 

large the ratio approaches asymptotically the value 1. 

Debye compares the light pressure so obtained with the 
gravitational attraction for a spherical particle of specific gravity 
s under the influence of the sun's radiation. He finds that if 0 
is the gravitational attraction, 

L ^ 4800 L 
G Xsp * W 

and to get a numerical estimate he takes \ = 600 /t/t, s = l. 

It appears that the ratio vanishes both for small and large 
values of p : it has a maximum value of about 20 for p = 1. 

In the case of a dielectric sphere with refractive index n, 
the expansion corresponding to (125) is 


L 1 Y ^ r p" n« - 29n* -f -f 1^20 1 

+ ^ 1_^ 15' (n* + 2)(2n» + 3) ‘"J 

(126), 

and is suitable for calculations only when r?p is small. 

Debye has drawn curves for ^ in the cases n = oo , n = 2, 

n = 1'5 and n = 1*33. When ?i = 2 the curve appears to have 
three maxima and two minima between the values p = 1 and 
p = 3. • 

The greatest value of is now about 2'6 ; the following 

table indicates when the light pressure exceeds the gravitational 
attraction, the numbers po and pi give the extreme values of p 
belonging to the range in which this is the case. 


p^ n* — 29n* -f 34/i^ -f 1^20 
15* (n2 4*2)(2n»-f 3) 


The maximum light pressure is just balanced by the 
gravitStional action when n is about equal to 1*33, the value 
for water : for smaller values of n gravitation prevails. 
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In the case of an absorbing spherical particle, the equation 
which takes the place of (126) is 


w~ 


(k + 2)» + ■ 


.(127). 


When a is small the light pressure and gravitational action 
are both of order a® and their ratio tends to a finite limit, hence 
for certain types of absorbing material there is no lower limit 
in the size of a particle below which gravitation exceeds the 

light pressure. Debye has drawn a curve for ^ for the case of 

a gold particle and finds that there is a maximum value for 
p == 1*5 nearly. 

The existence of a maximum value for ^ in the cases that 

W 

have been discussed appears to be due to the fact that the 
value of p for which the maximum occurs is very nearly equal 
to the real part of the complex value of p corresponding to one 
of the free damped vibrations*. The first electric vibration 
seems to be of chief importance in determining the position of 
the maximum. 

The determination of the limiting value of ^ for very small 

wave-lengths, i.e. for large values of p, is a matter of some 
difficulty, it depends on some expressions giving the behaviour 
of the Bessel functions for large values of n and p. These have 
been found by J. W. Nicholson f and P. Debye 
If p is real and n -f ^ < p, we have when p^ oo 

e-'M ' 

bn \P) — 


fn(p) = 


(sin T„)i 

("/•-I) 


cos 


(sin T„)' 


4 


.(128), 


* Cf. Debye, loc. cit.y and the similar remarks for the case of optical 
resonance by F. Pockels, Physik. Zeitschr. Bd. 5 (1904), p. 152. 

t British Association Reports ^ Dublin, 1908, p. 595 ; Phil. Mag. Vol , ' H (1906), 
p. 195 ; Vol. 14 (1907), p. 697 ; Vol. 16 (1908), p. 271 ; Vol. 18 (1909), p. 6. 
Math. Ann. Bd. 67 (1909), p. 535. 
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is an angle lying between 0 and — for which 


cos To = 


71 -f ^ 


and /o = sin To — Tq cos Tq. 

When n-h i > p and p -► oo , we have 

-^p/o \ 




(i sin To)^ 


fn(p) = 


Jpfo 


(i sin To)^ f 
where Tq is now the root of the equation 


.(129), 


cos To = 


n-hj 

P 


whose imaginary part has a negative sign. 

When n and z arc very nearly equal the values' of 5'n (p) and 
\frn(p) can be made to depend on Airy's integral and are much 
more complicated ; for these we must refer the reader to the 
original memoirs. 


24. Other problems which may be treated with the 
aid of polar coordinates. 

The diffraction of electric waves travelling round the earth 
is a problem of some importance which has heen discussed by 
H. M. Macdonald*, Lord Rayleigh f, H. Poincar6|, J. W. 
Nicholson§ and other writers. 

The calculations are very long and depend on the use of the 
formulae to which we have just referred. Rybcynski|| has 
recently treated the problem by a method due to March and 
has taken into account the finite conductivity of the earth. As 
we have already mentioned this was done by Zenneck and^ 
Sommerfeld for the case in which the earth s surface is treated 


* Froc. Roy. Soc. Yol. 71 (1903), p. 251 ; Vol. 72 (1904), p. 59 ; Vol. 90 (1914), 
p. 50; Fhil. Trans. A, Vol. 210 (1909), p. 113. 
t Proc. Roy. Soc. Vol. 72 (1904), p. 40. 

t Rend. Palermo (1910) ; Proc. Roy. Soc. Vol. 72 (1904), p. 42. 

§ PhU. Mag. Vol. 19 (1910), pp. 276, 435, 516, 757; Vol. 20 (1911), p. 157; 
Vol. 21 (1911), pp. 62, 281; Jahrb. d. draht. Telegr. Bd. 4 (1910), p. 20. 

11 Ann. d. Phys. Bd. 41 (1913). 
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as a plane. The results of Nicholson and Poincar^ indicate 
that diffraction round a perfectly conducting surface is not 
sufficient to explain the apparent bending of the electric waves 
round the earth's surface. A generally accepted opinion is that 
the ionisation of the atmosphere by the sun's rays is a very 
important factor in producing the observed effects*. 

The diffraction of a solitary wave or pulse by a spherical 
obstacle might be discussed with advantage. The evaluation 
of certain definite integrals involving Bessel functions, however, 
presents some formidable difficulties which probably account 
for the fact that the problem does not appear to have been 
solved. 

The scattering of electric waves by a perfectly conducting 
conical obstacle has been treated very briefly by H. S. Carslawf . 


EXAMPLES. 

1. Prove that 

1 + 2 Pn (cos 6 ) Pn (cos a) = log ^SCC* ~ COSGC* 0, 0<a<(9<7r. 

(C. Neumann.) 

2. Prove that 

COB ^ sin a) == ^ I 2 (2w-f 1) t" Pn (cos B) (cos a). 

(E. W. Hobson.) 

* Cf. the discussion at the British Association meeting, Dundee (1912), and 
an article by W. H. Eccles in the Year Book of Wireless Telegraphy (1913). 
Some quantitative experiments on long distance telegraphy have been made 
recently by L. W. Austin, who obtains an empirical relation between the 
magnitude of the current received and the distance between the two stations, 
Bulletin of the Bureau of Standards^ Vol. 7 (1911). 

I* Phil Mag, Vol. 20 (1910), p. 690. 
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§ 25. The wave-equation in Cylindrical Coordinates. 

If we put a? = p cos 0, y = p sin 0, the wave-equation becomes 


0^16 \du 1 d^u 0^w 1 dH 
dp^ ^ p dp ^ d(f)^ dz^ dt‘ 


,(130). 


Two particular solutions of this equation are^ suggested at 
once by the general solution of § 5 : they are* 



r2v 

u= F 

Jo 

z -h ip cos a, 

< — - sin a 

c J 

doL .. 

....(131), 

and 

II 

s 

-f ip cos a. 

4 P • ' 

^ - Sin a 

c 

doL .. 

....(132), 


respectively. The first of these represents a wave-function 
which is symmetrical round the axis of z and which reduces to 
27 rF(z, t) when p = 0. It gives us at once the formulae 


r"P„(yiA) = 2^j^ (^ + ^>cosa)'*(^a = •••(133), 

("-s ■■■(**«• 

where n is an integerf, and many other interesting formiltae 
may be written down by simply choosing different wave- 
functions that are symmetrical round the axis. 


* The first of these is an obvious generalisation of a formula given by 
D. Ed#ardes, Educational Times^ Oct. (1904). 

t The formula is also true under certain limitations when n is not an integer. 
See Hobson, Phil, Tram. A, Vol. 187 (1896). 
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For instance, 

_2r^ 

ttJq z-^-ipco&d r + z 

1 log (z 4- ip cos 0) - log 2^0 _ J- 1 ^ 4- ^ 0 -f.fi 

TrJo - -S'o + ip cos 0 ““fi v-VZq- R 

( JJ2 = p2 4- - ^To)®, 2r > 2 : 0 ). . .(136), 

_L P*" sin A; 4- ip cos g) ^<jfcp8ina ^ si nArr 
27rjo ^4-^pcosa r 


There is another general formula for a wave-function sym- 
metrical about an axis, viz. 

= J fi — - cosh a, z — p sinh a^da . . .(138), 

where the function F is of such a nature that 

1 , . , dF 

- sinh a ^ 4- cosh a 
c ot dz 


vanishes when a = 4 00 . As a particular instance of this we 
have the function* 

fl = ^ J” F ^ cosh a) da (139), 

which may be regarded as the cylindrical wave-potential for a 
line source of strength F (t) along the axis of z. 

A peculiarity! two-dimensional propagation of waves 

is the existence of ,a “ tail to the disturbance when F (t) is zero 

for ^ < 0 and ^ > r, for if ^ > r 4 - - we have 

c 

n= J F^t — ^ cosh CL^doi (140), 

where p cosh 5 = c (i — r). It is clear that this expression for 
fl does not generally vanish. The wave-function (139) is thus 
essentially different from Euler’s wave-potential for a point 
source, viz. 



* Cf. H. Lamb, Hydrodynamics ^ pp. 281, 600 ; V. Volterra, Acta Math, 
t. 18 ; Levi-Civit^, Nuovo Cimento (4), t. 6 (1897). The formula is a particular 
case of a more general one given by Dr Hobson in 1891. 

t This was discussed by 0. Heaviside, Phil. Mag, (5), t. 26 (1888) ; Elec- 
trical Papers^ Vol. 2. See also Lamb, Hydrodynamics, p. 282. 
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V 

for in this case fl is zero for ^ > r -f - , provided the source is 
only active when 0 < ^ < t. 

V. Volterra* has obtained a number of elementary wave- 
functions of the form 

M = j ^’(s) (141). 

He finds that F must satisfy the differential equation 
s=(l-s“)^+s(2w-«“)^ + »(w- 1)^ = 0 ...(142). 

If we try to solve (130) by means of a function of the form 
u = Wp^ cos m (<^ — (^o), 

where W is independent of we find that W must satisfy the 
equation 

0^Tf,2m + 10ir d^W Id^W ^ 

w * r ^ T.J - ? w - " <“*>■ 

Solutions of this which are independent of z may be derived 
from the following formulaef, in which m > — 

IT = J f(t — ~cosaj sin^"^ a , da (144), 

W = I F cosh sinh-"* r) .dt) (145). 


There are, of course, certain limitatiohs concerning the 
behaviour of F (t) at infinity. The first formula enables us to 
determine the value of W when its value is known for p = 0. 

§ 26. Elementary solutions of Au + k'^u = 0. 

The differential equation 


dht 19^, 

0p2 ^ p dp 


1 dhi dht ,, 


possesses elementary solutions of the types 

Jm {p e cos m {<f> - (f>,) 


.(146) 

(147), 


Km ip Vv- ^*0e*‘** cos m {^-<po) (148), 


* icta Math. Vol. 18. 

t E. W. Hobson, Proc. London Math. Soc. Vol. 22 (1891), p. 431. The 
first solution is due to Poisson. 
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where X, h and m are real or complex arbitrary constants. The 
first solution may of course be generalised into 

M= r J^{p^¥n^)e*''^f{h)hdh (149), 

J a 


and a similar remark applies to the second. If we wish to 
express a given wave-function in the first of these forms, the 
following inversion formula due to Hankel is particularly 
useful*. 


If 

then 


Jo 

f{t)= [ Jra (^0 ^ {^) 


(150). 


Let us use this formula to express in the form (149) 

when m = 0. Since the representation should be valid for 2 ? = 0, 
we find on putting A?* 4- = X^, that 

jk(i rao 

— =1 *To (Xp)y'(A) XdX 
p Jo 

k>0. 


Hence we obtain Sommerfeld’s formula | 





«/o(Xp) 


XdX 


(151), 


the upper or lower sign being taken according as ^ > 0. A 
more complete proof is obtained by applying Hankels inversion 
formula to the equation 


/ Jo (Xp) 
J 0 


gifcr 


pdp 

r ” Vx* -~k^ 


- x^ 


I 

[...(152), 


which is established in Prof. Lamb’s paper. 


* See Gray and Mathews, Treatise on Bessel Emictions (1895), p. 80 ; 

N. Nielsen, Handhmh der Theorie der Cylinderfunktionen (1904), p. 366. 

t Ann. d. Physik, Bd. 28 (1909), p. 683. Some analogous formulp.e are 
given by H. Lamb, Proc, London Math. Soc. Ser. 2, Vol. 7 (1909), p. 140 ; 

O. Heaviside, Electrical Papers, Vol. 2, p. 478 ; N. Sonin, Math. Ann. Bd. 16. 
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A few more formulae will now be written down to illustrate 
the method of generalisation by integration with regard to 
a variable parameter 

gikR 

Ko (p V V - ic^) cos \ (z- b) = ij cos X (a - 6) da 

(X^>k\ R^ = p^+(z-ay), 

If m is zero or a positive integer and n is a positive integer, 

,-i;i ■ r (» - m + 1 j /, < V> >•■‘1^. »>«. 

i - 2f (d'm + 1)/, ‘ 

(Hobson.) 

§ 27. The propagation of electric waves on a semi- 
infinite solid bounded by a plane surface*. 

In this problem the surface of the earth is regarded as an 
infinite plane and the waves are supposed to be generated by an 
antenna, of which one portion is vertical and the other horizontal. 

Let us assume that the electric and magnetic forces are the 
real parts of vectors of the form if respectively, 

then if M = H ivEy where ^ we may satisfy Max- 

~ fJLQ)'' 

well’s equations by putting 

ilf = rot n ± ^ grad div H ±icll (lo3), 

where AH + i^H = 0, -h ificoo- (154). 

To imitate the action of the antenna we shall place two 
vibrating doublets at a point at distance a from the plane. If 
one of these vibrates vertically and the other horizontally, we 
may put for the primary radiation no = (Paj, 0, P^), where 

pikR pikR 

P* = -5^, R^ = p> + (z-af...Xv^/, 

and the axis of z is vertical. We write B with a negative sign 
to signify that the horizontal branch of the antenna is drawn 
in the negative a?-direction. 

♦ A. Sommerfeld, Ann, d, Phy$, Bd. 28 (1909), p. 665 ; H. v. Hoerschlemann, 
Jahrb, d, draht, Teleg, Bd. 5 (1912), pp. 14, 188 ; Dmertationy Munich (1911). 
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A convenient expression for IIo is obtained by using Sommer- 
feld's equation 

0<z<a (156), 

where 1=^ — k\ Appropriate functions IIi, Ha for the re- 

flected and transmitted disturbances are obtained by putting 

111= S cosw<^ I Fn{\)dX 

n=0 Jo 

^<a...(157), 

00 r* 

112= S C08 n<f> Jn(Xp) On(X) d\ 

n«0 Jo 

where m — and h is the value of k in the second 

medium. The functions (\), On (X) are vectors with com- 
ponents [/n(X), 0, irn(X)l [ffn(XX 0, Xn(X)] respectively. 

We can satisfy the condition that the tangential components 
of the electric and magnetic forces should be continuous at the 
surface of sepaT:ation of the two media by putting for ^ = 0 

ITo + 111 = Ilg 

rot (Ho -f- Hi) = rot Da 

idiv(n,+n,)=ldivn, 

where p> has been taken to be unity for both media. 

Substituting the integral expressions for Ho, IIi, ITain these 
equations and equg^Jing to zero the coefficients of functions 
of type Jn{Xp) in the resulting integral equation, we obtain 
the system of equations 


/„ e-^ - j e-'“ = go 


A 

+ -j e~^ = Xo e"™*. 


■ Ifo e~^ — Be~^ — mgo 

- Ifn = mgn e“™“, 

1 

0 

1 

— ^ Y e~”*® 

-^aXoe , 


(:^e lO'-lyfr.e ^ 
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where the last equation has been simplified with the aid of the 
relations /n = 5 ^n = 0 (n>0) which are a consequence of the 
previous equations. 

Solving these equations we eventually find that if 

n, = (Q*, 0, = 0, 

B,— 


Qz^A 


Q—l (z+a) 


hH — k^mXdX 
hH -f k^m I 


j^Jo(Xp) 

+ 2B c, 




hH + k^m 


f J,(\p)e^ 

J 0 

+ 2B cos <i> [ Jj (Xp) 
Jo 


(h^-k^)X^dX 
(I + m) (hH 4- k^rn) * 


The directed effect '' depends on the presence of the terms 
involving cos <f) in the expressions for and Rz. Now when 
<7 = 00 for the second medium, /t = oo , and these terms vanish 
altogether ; hence the possibility of directing the energy of the 
radiation sent out from the bent antenna is due to the im- 
perfect conductivity of the earth. Von iHLoerschlemann has 
given a numerical discussion of the above formulae but the 
investigation is too long to be inserted here. 


§ 28. Propagation of electromagnetic waves along a 
straight wire of circular cross-section*. 

Let us consider the symmetrical case when the electric 
force at any point is in a plane through the axis of. 

* H. Hertz, Electric JFaves ; J. J. Thomson, Froc. London Math. Soc. 
Vol. 17 (1886), p. 310; Recent Researches^ § 259; A. Sommerfeld, Ann. d, 
Physikf Bd. 67 (1899), p. 233; Gray and Mathews, Bessel Functions^ Ch. 13; 
M. Abraham, Encykl. d. Math. Wiss. Band V. 2, Heft 3 (1910), p. 526; 
J. Lp^or, Proc. of the 5th hit. Congress of Mathematicians, Vol. 1 (1912), p. 206. 
The problem considered in this last paper is chiefly that of alternating currents, 
viz. the forced alternations of flow produced by a uniform periodic electric force. 
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and the magnetic force is in circles at right angles to this 
plane. The field equations are then of the types 


p dp^''“*'’ c dt 
fjL dH^ _ dEp ^ dEz 
c dt dz dp 

These may be satisfied by putting 
pdp' 


dz 


(159). 


E. 


(p^?) 

, E, 

V dp J 

> p 

€ a^n 

_o-an 

c dtdp 

c dp 


dpdz I 


.(160), 


where IT satisfies the equation 


€/i3^n orfldYl 


an 


^ 1^/ 

dz^ pdp\ dp 


= An 


.(161). 


Putting n = u, we find that Aw -f k^u = 0, where 




€fjL(o^ ~ ifiaxr 


We now assume that for points outside the wire* 
u = A Kq (p 

and that for points inside the wire 

li = Jq {ip Vv — h^)j 

where h is the valup of k inside the wire. These assumptions 
are made for the purpose of determining the periods and rates 
of decay of electric waves that can travel along the wire and 
maintain their own field. The first solution is chosen so as 
to make the flow of energy negligible at an infinite distance 
from the wire so that the system is self-contained; and to 
ensure this it is necessary to suppose that the real part of 
is positive. The second solution is chosen so as to 
make the electric and magnetic forces finite on the axis of 
the wire (/> = 0). 

Let p = a be the equation of the surface of the wire, then 


* We follow here the work of Sommerfeld in which it is supposed that there 
are no conductors outside the wire. Sir J. J. Thomson allows for the p^''sence 
of external conductors by supposing the dielectric surrounding the wire to be 
bounded by a cylindrical conductor having the same axis as the wire. 
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the continuity of the tangential components of the electric and 
magnetic forces requires that when p — a 


pKo{p^X^-k^) pJo(ip^X^-h^) I 


( 162 ), 


Av ^ Ko{ps/\^- — Bvi~ J^{ip^\^ — h^) 

op op 

V and Vi being the values of the quantity <T-\-im€ at points 

outside and inside the wire. The elimination of A and B gives 

rise to a transcendental equation for the determination of X. 

The total current flowing along the wire is 


J = 27r ( a-Ezpdp 

J 0 


= - 27 rae-^^^pB ~ Jq {i p V V - A*) for p = a. 

On the other hand the electric force at the surface 
of the wire is 


hz — ^ 

pdp 


pBe"*^ _ (^ip V A* - h*) 


for p = a 


= - J5 (X* - h^) 0^^ Jo {ia VX^ - 
Hence, if we put 

e,=j\R+'^L-. 

where R and L are the resistance and self-induction of the wire, 
we get 


iJ- 


icD j _ X^ — A* 


~ J o { x) " 
_xJ„'{x) 


x=TaVA*-/i“ 


The roots of the transcendental equation have been discussed by 
Sommerfeld who used a method of successive approximations. 
The values of X are complex as the waves which travel along 
the wire are damped owing to the imperfect conductance of the 
wire. It appears that when the disturbance does not penetrate 
far into the wire, the damping is small and so the v^^ity of 
propagation is very nearly equal to the velocity of light. 
When, however, the field does soak into the wire to some 
extent, the damping is of course considerable and so the wave 
travels with a velocity a little less than that of light. In the 
first i5ase the real part of is large, in the second case it 

is small. 
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In Lecher’s arrangement* there are two conjugate parallel 
wires between which the waves travel, consequently the field is 
not symmetrical round the axis of one of the wires. This case 
has been discussed by G. Mief with the aid of bi-polar 
coordinates, and the lines of force have been studied by W. B. 
Morton f. The latter also considers the case of n parallel wires 
passing through the corners of a regular polygon §. 

The mathematical analysis for the case of a curved or twisted 
wire has not yet been fully developed. The important case of 
a spiral wire has, however, been discussed by H. C. Pock- 
lington|| and J. W. Nicholson If. The latter gives numerous 
references to the literature of the subject. D. Hondros** has 
recently discussed the propagation of some types of unsym- 
metrical waves along a single wire. The electromagnetic 
theory of an electric cable has been given by Sir Joseph 
Thomson ff and F. Harms The latter considers the case 
when the outer conductor of a cable is replaced by air. 

§ 29. Other problems which may be treated with 
cylindrical coordinates. 

The diffraction of electromagnetic waves by a cylindrical 
obstacle has been discussed by Lord Rayleigh §§, W. Seitz ||||, 
W. IgnatowskylFII, P. Debye*f , C. Schaefer*^ , and J. W. Nichol- 
son*§. Schaefer has made an extensive study of the case of a 

I 

* Ann. d. Phys. Bd. 41 (1890), p. 850. See also D. Mazzotto, II Nuovo 
Cimento (4), t. 6 (1897), p. 172. 

t Ann, d. Phys. Bd. 2 (1900), p. 201. The case in which the capacity 
is small is discussed at length by J. W. Nicholson, Phil. Mag. Feb. — Sept. (1909). 
The current is supposed to flow along one wire and return along the other. 
t Phil. Mag. Vol. 50 (1900), p. 605 ; Vol. 4 (1902), p. 302. 

§ Phil. Mag. Vol. 1 (1901), p. 563. 

II Proc, Camb. Phil. Soc. Vol. 9 (1897), p. 324. 

H Phil. Mag. Vol. 19 (1910), p. 77. 

d. Phys. Bd. 30 (1909), p. 905; Dissertation, Munich (1909). 
tt Proc. Roy. Soc. Vol. 46 (1889), p. 1 ; Recent Researches, p. 262. 

XX Ann. d. Phys. Bd. 23 (1907), p. 44. 

§§ Phil. Mag. Vol. 12 (1881), p. 81. 
nil Ann. d. Phys. Bd. 16 (1905), p. 746. 

HIT Ann. d. Phys. Bd. 18 (1906), p. 495. 

*t Phys. Zeitschr. (1908), p. 775. 

*x Ann. d. Phys. Bd. 31 (1910), p. 462. 

*§ Proc. London Math. Soc. Ser. 2, Vol. 11, p. 104. 
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dielectric cylinder and his results have been tested experimentally 
by Grossmann. This is one of the few cases in which the 
influence of the material properties of the obstacle has been 
taken into account in the mathematical treatment of a diffrac- 
tion problem; the necessity of doing this has been clearly 
indicated by some of the experimental results. 

P. Debye has studied the diffraction problem with reference 
to the theory of the rainbow. This had been done for the case 
of the sphere by L. Lorenz long ago*. 

Electrical vibrations in regions bounded by cylinders of 
various shapes have been studied by Sir J. J. Thomsonf, Lord 
RayleighJ, Sir J. Larmor§, R. H. Weber ||, A. KalahneU, and 
J. W. Nicholson**. The latter has also calculated the pressure 
exerted by a train of plane electromagnetic waves on a perfectly 
conducting cylinder ff. 


EXAMPLES. 

1. An infinitely long metal cylinder of specific' conductivity o* 
and permeability ft, bounded by the surface r=/2, is surrounded by a 
dielectric of specific inductive capacity f. A train of waves, in which the 
electric force is perpendicular to the cylinder and to magnetic force and 
if undisturbed would be represented by the real part of 
is passing in the dielectric. Prove that the magnetic force inside 
the cylinder and the part Hi of the magnetic force outside representing 
the scattered wave, are given by the real parts of 

= 2 {kr) cos 2 {hr) cos 

m m 


where -47r/xorfa), and 


4Trcr 


KJm {hli) - 

j„: m, 


* Oeuvres scientijiques ^ pp. 405 — 502. See also Gans and Happel, loc. cit. 
(p. 44). The most recent paper on the rainbow is by W. Mobius, Ann. d. Phys. 
Bd. 33 (1910). . 

t Recent Researches, p. 344. 

X Phil. Mag. (5), Vol. 43, p. 125. 

§ Proc. London Math. Soc. (1), Vol. 26, p. 119. 

li Hahilitationsschrift, Heidelberg (1902) ; Ann. d. Phys. (4), Bd. 8 (1902), 

p. 721. 

1[ Ann. d. Phys. (4), Bd. 19 (1906), pp. 80, 879 ; Bd. 18 (1905), p. 92. 

** ^hil. Mag. Aug. (1905), May (1906). 

ft Proc. London Math. Soc, Ser. 2, Vol. 11, p. 104. 
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for positive integral values of w. The constants here have reference to the 
electromagnetic system of units. 

(Cambr. Math. Tripos, Part II, 1905.) 

2. Plane electromagnetic waves represented by 

fall upon the perfectly conducting cylinder p=a. Prove that in the 
scattered field 

(P. Debye.) 

3. The wave-potential for a circular ring of point sources is given by 

0 =J" M\p)J(,{\a)\dK, 2 > 0 . 

(A. G. Webster.) 

4. The wave-function 

Q=e^^‘jy-^Jo {p •J¥+^) cos a\ dX 

is zero at points on the plane 2=0 which lie inside the circle p^<a^; for 
2=0, p^>a^ its value is (p^— a*)”"^ cos (Sonin.) 

5. If /?^=p2 + ( 0 -f asinh w)2, ^=6^*, the integral 

j jQ(\p)\d'\ 

represents the function outside a paraboloid of revolution 

whose focus is at the singularity 2 = —a sinh 24, p = 0, and which passes 
through the circle p = a, 2 = 0. It is zero inside the paraboloid. 


6. The circuital relations in cylindrical coordinates are 


dHz 


p 0£’p 

c dt 


p dEz ( tj \ 


p^_ dHp_ dH. 
c d( dz ^ 9p ’ 


Op ■■ ■ ■ 0<p oz 

and similar equations in which E is replaced \rj —H and H by E. 

7. If Jr= or cos o)^ — y sin o)^, F=4r sin ca^-hy cos o)^, Z—z—vt^ where v 
and ft) are constants, the function Q^F{Xy F, ^) is a wave-function if F 
satisfies the partial differential equation 

A /'i - aA "Vf +’^-(x + r^-l\ 

V c2 S-P"*" V ) SF^"^^ c* ^^dXhY'^c^ \r ^dv) 


/_ d^E . 
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Obtain particular solutions of the form 

Q ^ X(Z-H) +<»(♦- „«) 

f V^\ 

where p^ = \‘^ I ^ ) + ~~^' and \ A^B are arbitrary constants. 

8. An oscillatory current is induced on a circular wire of radius a 
excited by a uniform electric force acting on its surface from the 

surrounding medium. Obtain expressions for the inductance and resistance 
of the wire per unit length when the wire is regarded as straight and no 
disturbing conductor is near. 



CHAPTER V 


THE PROBLEM OF DIFFRACTION 


§ 30. Multiform solutions of the wave- equation*. 

The wave-functions required to solve many of the boundary 
problems of Mathematical Physics are not single- valued func- 
tions of X, y, t in an ordinary space. We may, however, 
regard them as single-valued functions in a Riemann’s space. 
This is a simple generalisation of the Riemann’s surface of the 
theory of functions of a complex variable f ; every plane 
section of the Riemann’s space is in fact a Riemann’s surface. 
Instead of branch lines and branch points we have branch 
membranes and branch curves. Thus in the physical problem 
of the diffraction of light through a circular hole in a screen, 
the boundary of the shadow of the screen is the branch mem- 
brane and the edge of the hole the branch curve. 

We shall commence by finding a multiform solution of 
the equation 


d^u d^u 
da^ dy^ 


4 * k^u = 0 


(103). 


The fundamental solution u = (x cos a + y sin a) _ ^ikp cos - a) 

of period 27r and can be expanded in the form 

= jQ{kp) 4- 2 S i^Jn{kp) cos 7i(<^ - a). ..(164). 

1 


* This theory is due to A. Sommerfeld, Math, Ann. Bd. 45 (1894), Bd. 47 
(1896); Zeitschr. fiir Math. u. Phys. Bd. 46 (1901); Proc. London Math. Soc. 
(1), Vu7. 28 (1897), p. 417. It has been developed by H. S. Carslaw, Proc. 
London Math. Soc. (1), Vol. 30, p. 121 : (2), Vol. 8, p. 365 ; Phil. Mag. Vol. 5 
(1903), p. 374: Vol. 20 (1910), p. 690; Fourier's Series and Integrals^ Ch. 18; 
W. Voigt, G'dtt. Nachr. (1899) ; E. W. Hobson, Camh. Phil, Trans. Vol. 18 
(1900), p. 277. Different methods have been used by H. M. Macdonald, Electric 
Waves^ Appendix D; K. Schwarzschild, Math. Ann. Bd. 55 (1902), p. 177; 
Proc. London Math. Soc. Vol. 26 (1895), p. 156; C. W. Oseen, Arkiv far mat. 
Bd. 1 (1904), Bd. 2 (1905). 

t See Harkness and Morley’s Theory of Functions (1893), Ch. 6* 
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A solution of period 2m7r may evidently be constructed by 
writing down a series of the form 

«> ^ 

2 CbnJn i^p) COS ~ (<f> — Ct), 

^ = m ^ 

where the an s are suitable constants. The solution that seems 
the most natural extension of (164) is 

00 !L 

Pm (p, <f>, </>o) = *4 (^p) + 2 2 1”* (kp) COS ^(<f> — .(165). 

To sum this series when m-2, we transform the terms for 
which n is odd by means of the equation 


rt 


^ n (fcp) = 
2 


(kp)' 


‘r(i>r(!y 


f,- 

ii + ly ^ 


Qikp cos a 


2 ; 


sin" OL . d(x. 


fv -> C 08 * tt + tkp cos a 


sin a . da 


Summing the two series separately we find that 

F, (P, 4>, <#.,) = +/(</) - <^„) +/(<^o - 4>), 

where 

with T= cos 2 . = (- 2 O* ^ ~ • 

Now f^''^ 2 

hence we may write 

Jf’a (p, <^, ^0) = ~ 


where t = “J^kp cos ^ (<j>~~ ^o)- 

With the aid of the function we can solve some problems 
on the diffraction of plane electromagnetic waves by a semi- 
infinite plane bounded by a straight edge. Let us consider 

6—2 
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the case of a totally reflecting screen*. If the electric force 
in the incident wave is parallel to the edge of the screen, the 
electric force u — Ez for the total disturbance must vanish over 
both faces of the screen and must satisfy the differential equa- 
tion (163). These conditions are fulfilled by taking 

Ez^F2(Py(t>y<f>o)-F2(py(f>y--(f>o) ( 167 ). 


This value of Ez also satisfies the right conditions at infinity. 
To prove this we must find an asymptotic expression for F^ 
when p is large. 

Now when r > 0, we have the asymptotic expansion f 


j e dX ^ ' 


2iT“ (2tT7 


while when t< 0, f has a similar asymptotic expansion with the 
. -00 

sign changed. This means that when cos ^ (<^ — </>o) > e > 0 we 
have 

F, (p, 4>, <#>.) ~ cos ^ ^ _ 

V Zirkp cos — <]f)o) 


while when cos ^ (<^ — </>o) < ^ < 0 there is a similar asymptotic 
expansion in which the first term is missing. It thus appears 
that the electric force in the geometric shadow vanishes at 
infinity to the order 

If the magnetic force in the incident wave is parallel to the 
edge of the screen, the magnetic force u — Hz must satisfy the 

• • • 

differential equation (163) and be such that ^—0 over both 

(j(p 

faces of the screen. The conditions are fulfilled by putting 

Hz — F.2{p, -f F^ip, (f>, - (f>o) (168). 

Prof. H. M. Macdonald has shown that the solution of 
a problem concerning a perfectly absorbing body can be made 
to depend on the solution of two allied problems “ A per- 


* The incident waves are supposed to come in a direction for which 
<l>=:w-\-<pQ. An approximate solution of this problem was given by H. Poincar6, 
Acta Mathematical Bd. 16, p. 297 ; Bd. 20, p. 313. 
t Bromwich’s Infinite Series, p. 328. 

X Phil Trans. A, Vol. 212 (1912), p. 337 ; Proc. London Math. Soc. (2), 
Vol. 12 (1913). 
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fectly absorbing body may be regarded as a body which is 
incapable of supporting either electric or magnetic force ; 
hence if C is the electric current distribution on the surface 
of the body when it is supposed to be perfectly conducting, 
and 0' is the magnetic current distribution on the surface of 
the body when it is supposed to be incapable of supporting 
magnetic force, the superposition of these two distributions 
gives the electric and magnetic current distributions on the 
surface of the body when it is perfectly absorbing and the 
amplitude of the incident waves is doubled.” 

Now if we suppose our screen to be incapable of supporting 
magnetic force, the boundary condition is that the tangential 
component of the magnetic force should vanish. When the 


electric force is parallel to the axis of z 


^ d(l> 


the screen. Hence 


must vanish over 


E, = <f>o) + F,(p, (/), - (^o) (169). 

The solution for a totally absorbing screen is thus simply * 

E,^F,(p,cf>,<f>o) (170). 

Similarly, it can be shown that when the magnetic force 
is parallel to the axis of Zy the solution for a perfectly 
absorbing screen is 

H,^F,(p><f>><f>o) (171). 


Prof. Lambt has discussed the case of perpendicular incidence 
with the aid of the parabolic substitution 

^ = p^cosl<f>, v = P^ 
x = y-2^r], p = + 

The curves ^ = const., r}~c nst. are confocal parabolas, = 0 
is the screen. A solution of Maxwells equations is obtained 
by writing 



dEx _ du 
dt ^ dy' 


where 


* W ^oigt, G'ott. Nachr. 
screen. 

t Proc. London Math, Soc. 


Hy = 0y 

dEy du dE;; _ ^ 
di ^ dx* dt 
d“U ^ dhi _ 1 
dx^ dy^ c^dt^ 


.^^73), 


(174). 

(1899), p. 1, discusses the case ol an absorbing 


(2), Vol. 4 (1907), p. 190 ; Vol. 8 (1910), p. 422. 
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Starting with Poisson’s wave-function* x = 4 ^ p)> 

let us put 

du 

0*=^- 

Transforming to the coordinates t;, we obtain 

= 

Solving this partial differential equation by Lagrange’s method 
and adjusting the complementary function so that the boundary 
024 

condition = 0 for ?? = 0 is satisfied, we obtain 
dr) 

r^+n 

W = Jo f(ct + y-^^)d^+j^ f(ct-y-^^)d^ 

+iiF(ct+y)+^F(ct - 2 /). ..( 176 ), 
where Fh an arbitrary function. If the boundary condition is 
u = 0 for 7}=^0 the sign of the second term must be changed. 
It is easy to verify that each of the integrals represents a wave- 
function. 


Let us now put /(x) = F' {x — v^) dv, 

and make the substitution f == o- cos a, v — or sin a, then after 
a little reduction we find that*)* 


u • 


1 

-F{ct + y) + sec*^ a] da 

rr J 0 

n 

1 f 2 

d" " j ^ ~ (p ““ y) a] da x<0, y^O 


TT 

1 f 2 

u = F {ct-\-y) + F {ct —y) ^ + 2/ ~ (p + y) sec*-^ a] da 

ttJ 0 

TT 

1 /*¥ 

— — j Flct—y— (p—y) sec^ a]da 

TT 

1 r2 

u 1 F[ct-\-y — {p -{-y) sec* a] da 
TTJo 

TT 

1 

-h - Flct — y — {p — y) sec* a] d 
TT jo 


x>() 

y>0 


X >0 

y>0 


.( 177 ). 

* Journal de VFxole Polytechnique, Cfth. 19, t. 12 (1823). Seo also V. 
Volterra, Acta Math. t. 18 ; Hantzschel, Reduction der Potentialgleichung ^ Ch. 1. 
t These formulae are not given in Prof. Lamb’s paper. 
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Each of the integrals represents a wave-function provided 
F{± <X))z=:0 and 

TT 

tan oiF' [ct ±y'-(p ± y) sec^aj — » 0 as ^ "2 • 


In these circumstances we have the solution of the diffraction 
problem for the case when the initial disturbance is represented 
by u=^F {ct y\ the magnetic force being parallel to the 
axis of z. By suitably choosing F we can deal with the case 
of a solitary wave. 

A new method of solving the problem of diffraction by 
a straight edge has been given recently by Oseen*. 

Problems connected with a wedge have been treated success- 
fully by Sommerfeld and other writers by using a certain type 
of contour integral. The fundamental solution of (163) is now 

w = f ^ (178), 

IttiJ da V ^ ^ 


ia l<f>o 

where v = e ” — e ” and the path of integration is %, simple con 
tour which starts from 00 1 -h 7 and goes to 00 i -f y witho’ 
crossing the real axis. The quantities 7 , y are subject to ^ 
inequalities 

27r > 7 > TT, 0 > 7 ' > — TT. 

This function is multiform and of period 2 n 7 r, but 
7i-sheeted Riemann’s surface with the origin as branr 
and the line </> = - (tt — <^o) Rs branch -secti^pn, it is 
With the aid of this function a number of 
problems may be solved. 

Thus in the case of a perfectly conducting pr 
27r — a if the electric force in the incident waves 
the edge of the screen and is represented by 
of the expression 

[ct4-p cos (<f> — 

it can be shown by an extension of the methr 


for the total disturbance the electric force 
the expression t 




ikct 


I 


Q%kp cos ^ 1 q 

'^<6 dc 


* Arkiv for matematik, r 

t Macdonald, Electric Wi ^ 
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where 

2 /; = cos ^ — cos — <^o)» = cos ^ — cos ~ <^o)> 


and the path of integration is the same as before. 

In the associated problem when the magnetic force in 
the incident wave is parallel to the axis of z, the magnetic 
force for the total disturbance is the real part of the expression 

H, = J e’’*" f ~ log (ww,) d^. (180). 

These solutions and the solutions of analogous problems have 
been discussed by W. H. Jackson*, H. M. Macdonald f, F. 
Reiche J, A. Wiegrefe §, and other writers. 


§31. Elliptic coordinates 1 1. 

If we put 

(T = cosh ft) cos y = sinh eo sin ^ (181), 

the differential equation (163) becomes 

^ (cosh^ - cos^ y) M = 0 (182). 

dft)* 

elementary solutions are now of the form 

u = E{(o)F{x)> 

E and F satisfy the equations of the elliptic cylinder IF 


~ + ()t^cosh“®+p)£: = o| 
^-(*=cos>x+i>)i^=0j 


(183). 


solutions of these differential equations have been 
tly by Prof. Whittaker**. 


Math. Soc. Ser. 2, Vol. 1 (1904), p. 393. 

913), p. 430. 

37 (1912), p. 131, 

^2), p. 449. 

inn. Vol. 1 (1869) ; Mathieu, Louville's Journal, Ser. 2, 
ein, Hoppe's Archiv (2), t. 14, p. 170 ; R. C. Maclaurin, 
17 (1898), p. 41. 

Heine, Handbuch der Kugelfunktionen ; Lindemann, 
^"chel, Zeitschr, f. Math. u. Phys. Vol. 31, p. 25 
urnal, Ser. 2, Vol. 13 (1868). 
e (1912). 
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Elliptic coordinates are appropriate for the solution of 
problems connected with the scattering of electromagnetic 
waves by an elliptic cylinder*. 

W. Wien has suggested f that the problem of the diffraction 
of light through a straight slit in a screen J may be treated 
with the aid of elliptic coordinates by regarding the screen as a 
limiting case of a hyperbolic cylinder. 

H. Weber § has shown that when the elliptic and 

parabolic substitutions are the only transformations which lead 
to elementary solutions of the equation (163). For further 
properties of this differential equation we may refer to Pockels, 
Tiber die partielle Differentialgleichung Au + k“U = 0, Teubner, 
Leipzig (1891), and to Lord Rayleighs Theorif of Sound. 

§ 32. Other diffraction problems. 

The diffraction of light and electric waves by a grating of 
wires is a problem of importance, but the mathematical treat- 
ment is very difficult and the theories that have been given so 
far are of an approximate character. Sir J. J. Thomson|| has 
discussed the theory of Hertz's grating H which consists of 
a number of parallel equidistant metal wires. When electric 
waves whose wave-length is large compared with the distanc 
between the wires fall normally on the grating, they p? 
through if the electric force is at right angles to the wires 
are reflected if the electric force is parallel to the wires. ^ 
Lamb ** has considered the case of a grating which consi 
parallel strips of metal ; his theory has been supported I 

* See, for instance, K. Aichi, Proc. Tokyo Math. Phys. Soc. (2), 

(1908) ; B. Sieger, Ann. d. Physik (4), Bd. 27 (1908), p. 626. 

t Jahreshericht d. deutseh. Math. Perein, Bd. 15 (1906), p. 42. 

For this problem see K. Schwarzschild, loc. cit. ; Lord Ra 
Mag. Vol. 43 (1897), p. 259 ; Scientific Papers, Vol. 4, p. 283 ; P 
A, Vol. 89 (1913), p. 194. An interesting experimental result ho 
recently by P. Zeeman, Amsterdam Proc., Nov. 28 (1912), p. 59 

§ Math. Ann. Bd. 1. See also Hantzschel, Reduction der ‘ 
p. 137. 

II Recent Researches (1893), p. 425. 

IF Collected Works, Vol. 2, p. 190. For some recent r 
H. du Bois and H. Rubens, Ann, d, Phys. Bd. 35 (1911 
Phys. Review, Jan. (1913). 

** Proc. London Math. Soc. Ser. 1, Vol. 29 (1898) 
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experimental work of Cl. Schaefer*, J. Langwitz*, and G. H. 
Thomson f. 

Lord Rayleigh I has given an approximate electromagnetic 
theory of the action of a grating on waves of light and this 
theory has been extended by W. Voigt § so as to take into 
account the properties of the material of which the grating 
is made. Voigt's theory has been tested experimentally by 
B. Poganyll who gives an account of previous experimental 
work on the subject. 

The diffraction of light through a circular hole in a screen 
is a problem of interest to mathematicians which has yet to be 
solved IT. A promising method of attack is to regard the screen 
as the limiting case of a hyperboloid of revolution of one sheet. 


Examples. 1. 
ikct 


Prove that 


^ 2 (pj </>o) = j A 


ikcTj 

e dr 


) pcos((/)-<^o) 


2. If + ty = a cosh (co )» ^ r ) = \kae ^ » 

fi — — 

2 2 

'here the summation extends over all even integral values oinMv is even 
i over all odd integral values of 7i if v is odd. 

(J, H. Hartenstein, GrunerVs Archiv (2), t. 14, p. 170.) 


32a. The introduction and elimination of discontinuities. 

/ave-functions with singular lines or with singularities 
Ung along straight lines with the velocity of light may 
nes be employed with advantage in the solution of 
on problems. To illustrate the method to be adopted 
consider the diffraction of waves of sound by an 

Phys. Bd. 21 (1906), p. 587. The theory is developed from a new 
'Y C. Schaefer and F. Reiche, Ibid. Bd. 32 (1910), p. 677 ; Bd. 35 

' (1907), p. 365. 
c. A, Vol. 79 (1907), p. 532. 

911). II Ann. d. Phyt. Bd. 37 (1912), p. 257. 

lutions have been given by G. G. Stokes, Gainb. Phil. 
'•enz, Vidensh. Selsk. Skr.^ Copenhagen (18901 * H. A. 
Vol. 6 ; A. Grimpen, Dus. Kiel (1890) ; A. E. H. Love, 
(1901). 
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infinitely thin semi-infinite plane bounded by a straight edge*. 
We shall suppose that the waves are sent out from a stationary 
source and that the screen acts as a perfect reflectorf. 

Let the axis of z be taken along the edge of the screen and 
the axis of x in the plane of the screen at right angles to the 
edge. Let P be the source of sound and Q an arbitrary point 
on the edge of the screen. If the disturbance issuing from P 
were reflected according to the laws of geometrical optics, the 
total disturbance would be discontinuous in crossing two semi- 
infinite planes, each of which is bounded by the edge of the 
screen. The first of these planes is a boundary of the geometrical 
shadow, when continued across the edge of the screen it passes 
through P. The second plane is the boundary of the geometrical 
shadow for the optical image of P, viz. Pj. 

To obtain the correct solution of the diffraction problem we 
must add to the disturbance just described a second one having 
discontinuities which will annul the above-mentioned dis- 
continuities, the new disturbance must also be chosen so that 
the boundary condition is satisfied at the two faces of the 
screen. We shall now show that the required disturbance can 
be built up by superposition from elementary disturbances 
with singularities along lines such as PQ and P^Q produced. 

Let R be the distance of an arbitrary point {x, y, Zy t) 
from Q, then if (0, 0, are the coordinates of Q and c is the 
velocity of sound, we know that a function Sf type 

' 1 ^ / R X -h i]l \ 

w r-f+p/ 

satisfies the wave-equation. Let us choose the arbitrary function 
f in such a way that the expression becomes infinite along the 
line PQ produced and returns to its initial value when the 
point Xy y, is rotated twice round the edge of the screen. 
This last condition is added so as to enable us to sa^sfy the 
boundary condition. 

* This problem has been solved by H. S. Carslaw, Proc. London Math. Soc. 
Vol. 30 (1898), p. 121. A transformation of his solution suggested the method 
describea^ere. 

t This assumption is usually justifiable. In Prof. A. G. Webster’s experi- 
ments on the reflection of sound from the ground {Phys. Review y Vol. 28 (1909) » 
p. 65) it was found that the reflection is more than 90 7o- 
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Now if we write 

where u and (f> are real quantities, and use Eq, Uq, <^o to denote 
the values of E, u, <f> for the point P, it is easy to see that the 
function 

1 / E \ 

T^Fit I sec 4 r<<> - <i)o + i + t/o)l 


^ P “ j sec ^ - (/)o + i (ti 4- Ho)] 

satisfies the requirements, for it is periodic in (p with period 
47r and is infinite along PQ produced, where 

u=—Uo, ^ — (f>o + 7r. 

We now imagine sources corresponding to wave- functions of 
this type to be associated with each element df of the edge, 
and suppose the strength and phase of the source at Q to 
depend on its position relative to P in such a way that 




P + Pc 


where f(t) is the strength of the source at P at time t 
In this way we obtain an integral 


P 4" Pc 


% 

j sec ^ [(/> - </>o + i (w + '^o)], 


which will be shown to be discontinuous in the way required 
as the point (x, y, z) crosses the boundary of the shadow for P. 

In a similar way we can construct an integral 

- & /". m4‘- 5 »+♦.-«“ + “.)! 

which can be shown to be discontinuous in the way required as 
the point {Xy y, z) crosses the boundary of the geometrical 
shadow for Pi. 

Now let r, be the distances of the point a?, y, z respectively 
from P and Pj, then the velocity potential V of the total 
disturbsCiice is given by the following expressions in the different 
regions of space 

F=-F,-F, in-Sf,. 


ff.1- 
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The space Si is bounded by the screen and the limiting 
plane of the geometrical shadow for Pj, S 2 is bounded by the 
limiting planes of the shadows for P and Pj, S^ is bounded 
by the screen and the limiting plane of the shadow for P. 


dV 

The boundary condition is that ^ should be zero over the 

two faces of the screen and it is easy to verify that this 
condition is satisfied. To show that V is continuous for the 
whole of the space outside the screen and vanishes at infinity 
when the function f is finite, we shall transform the integrals 
Fi, V 2 to the forms given by Prof. Carslaw. To do this we put 
u + Uq=: b, then if p is the distance of a point from the axis 
of z, we have 


pe^=:z- ?+ P, poe^^ == ^0 - f + Po, 

R^^p^^{z^ tf, Po* = p,^ + (z, - f)*, 

= P - ^ 4- Poe~^* = Po- ^0+ 

PPq cosh (li + - 0 {zo - 0 + PPo, 

+ Po^ - ^oY + 2ppo cosh (w 4- = (P + Po)^ 

db = dti 4- dt/o = - df 4- . 

Hence it follows that 


- s/'. STJe. - i - ♦. + »). 


db P4"Po' 

.pTPo*^^ T 
On substituting the expression for P 4«Po in terms of b we 
obtain an integral which is equivalent to the one given by 
Prof. Carslaw. To see that it is discontinuous* we use and F_ 
to denote the values of the integral for <^ = tt 4* </)o + € and 
= TT 4- — € respectively, where e is a small quantity. The 

difference between these quantities may be regarded as a 
contour integral and can be evaluated by Cauchy’s theorem 
We may write 




I 


d( 


4f7ri J P 4* Pq 




R + R, 


cosec 


Taking the residue for | = 0, i.e. <^ = tt + 6 = 0, we get 

1 R + R,\ 1 r^ 






A more careful proof is given in Prof. Carslaw’s paper. 



94 


THE PROBLEM OF DIFFRACTION 


[CH. V 


for the conditions <^ = 7r4-^o» b = 0 imply that the radii 
R and R^ are in one straight line and so give r when added 
together. 

It is now clear that the integral Fj possesses the right type 
of discontinuity and a similar remark holds for the integral V 2 . 
The method can no doubt be modified so as to give solutions of 
other types of diffraction problems, the chief difficulty arises 
in the choice of a function which will satisfy the boundary 
conditions. At any rate the method suggests an interesting 
type of boundary problem in which the desired wave-functions 
have specified discontinuities instead of being continuous 
everywhere. This type of problem ought to be studied more 
completely. 

In the general problem of the diffraction round a moving 
object of the waves issuing from a moving source, the wave- 
functions that are derived by the methods of geometrical optics 
have discontinuities at a certain boundary which is the locus of 
points travelling along straight lines with the velocity of light. 

The points in question start from certain points of the 
moving object and move along tangents to the surface of the 
object, their paths being in fact continuations of the paths of 
particles that may be considered to have been emitted from the 
source. Indeed, if we imagine the source to emit particles in 
all directions as it moves about, the particles which just graze 
the moving object will, when they continue their rectilinear 
motion with the velocity of light, form the boundary at which 
the discontinuities arise. 

In Chapter viii we shall obtain a class of wave-functions 
with singularities moving along straight lines with the velocity 
of light. These functions seem to be just the ones that are 
required for the building up of wave-functions with dis- 
continuit/ies of the type just described. The problem of forming 
in this way the functions which will enable us to complete the 
solution of the diffraction problem is one which awaits solution. 



CHAPTER VI 


TRANSFORMATIONS OF COORDINATES APPROPRIATE FOR 
THE TREATMENT OF PROBLEMS CONNECTED WITH A 
SURFACE OF REVOLUTION 


§ 33. Spheroidal coordinates. 

Problems in which there is symmetry round the axis of z 
can often be treated with the aid of a substitution of the form* 


p+iz= /(a + f/9) (184). 

Taking a, /3, 4> as orthogonal coordinates, we have 

+ dy^ + dz^ = (da’‘ + d/3») + p’d(f>^- ■ -(185), 


and equations (18) of § 8 become 



= ± kM, 


= + AjiT/fl ' 

~ j ^ 


= ± hM^ 

55 ■•'Sis (5*') ■'■ 5 * 

II 


(186), 


where 


7-2 = ^ («lA) 
d (/), z) ■ 


* This substitution has been used in other branches of mathematical 
physics by C. Neumann, Theorie der Elektricitdts- iind Wdrme-V ertheiUing in 
einem Binge (1864); E. Mathieu, Cours de physique mathimatique (1873); 
A. Wangerin, Berliner Monatsherichte (1878) ; Hantzschel, Beduction der 
PotmtialrJuichung ; Michell, Mess, of Math. (1890) ; Basset, Hydrodynamics ^ 
VoL ,2, p. 8 ; F. H. Safiford, Amer, Joum. Vol. 21; Archiv der Math, Bd. 13 
(1908), p. The important developments on which the following analysis is 
founded are contained in papers to which we shall refer presently. 
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These equations may be satisfied by putting 

-asT), 

where H + i F is a solution of the partial differential 
equation 


1 /( 

The problem of finding the periods of free electrical oscilla- 
tions on a conducting spheroid is of considerable interest 
because a straight rod of circular cross-section can be regarded 
as approximately equivalent to a prolate spheroid whose major 
axis is relatively much longer than the minor axis. This 
problem has been treated very fully by M. Abraham* * * § , 
R. C. Maclaurinf, M. BrillouinJ, F. Ehrenhaft§ and J. W. 
Nicholson||. The effect of a spheroidal obstacle on a train of 
waves has been studied by K. F. HerzfeldlT. 

For prolate spheroids the appropriate substitution is** 

z-\’ip = a cosh (a -f tyS) (189), 

or p = a sinh a sin y9, z = a cosh a cos 

d ~ (cosh^ a — cos^ /S) (190). 

The partial differential equation is now 

4- — coth ® ^ ^ (cosh^ a — cos* /S) ft 

(191), 


1 /dp 9ft dp 9ft 

( 4* (-\“A --V A 


+ J,ft = 0...(188). 


and there are elementary solutions of the form il = A(a) B (/3) 
where A and B satisfy the differential equations 


d^A ,, dA 


(A;* a* cosh* 2 -\-X) A 


d^B dB 

— cot ^ ^+(\ + l(?a^ cos’ I3)B = 0 


> ...(192). 


* Disiertation, Berlin (1897) ; Ann. d. Phtjs. Bd. 66 (1898), p. 435 ; Math. 
Ann. Bd. 52 (1899), p. 81. 

t Cambr. Phil. Tram. Vol. 17 (1898-9), pp. 41—108. 

X Propagation de Vilectriciti (1904), Ch. vi. 

§ Wiener Berichte (1904), p. 273. 

II Phil. Mag. (1906). If Wiener Berichte (1911), p. 15^7. 

** Cf. Heine, Crelle, Bd. 26 (1843), p. 185 ; Kugelfnnktionen, Bd. 2, § 38 ; 
Lamb’s Hydrodynamics, p. 132. 
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These equations are discussed in some detail in the papers 
to which we have just referred, they may be reduced to a special 
form of an equation obtained by Prof. C. Niven* in a study of 
the conduction of heat in ellipsoids of revolution. 

For oblate spheroids the appropriate substitution is 

p ^iz= a cosh (a -f ^/3) (193), 

giving p = tt cosh a cos z=^ a sinh a sin yS, 

f - ; (cosh^ a — cos^ /9). 
d (a, /3) 

The surfaces /3 = const, are now hyperboloids of one sheet, 
the surface yS = 0 can be regarded as the surface of a screen 
which is pierced by a circular hole of radius a. 

The partial differential equation for is now 


-h ^ - tanh a ^ -f tan yS ^ (sinh^ a -f sin^ y5) fl = 0, 

cot (7p“ OOL op 

and there are elementary solutions of the form fl = (a) 5 (yS) 

where 


d?A 

dc? 

d^B 


d A 

— tanh CL . ~T- 4- (X + dk^ sinh^ a) = 0 
da ^ 

+ tan/3. ^ +(a'*A;*sin";8-\) = 0 j 


...(194). 


When n is independent of t, k = 0 and the elementary 
solutions are of the form 


n^fP„{^)d^iPn{v)dv, ^ = cosha, 7; = cos ^ ...(195) 
for prolate spheroids, and of the form 

n=fPn{^)d^JBn{v)dv> ? = isinha, 7? = sin ^...(196) 

for oblate spheroids. In either of these solutions a function Pn 
can be replaced by Qn- The corresponding solutions of Ifaplace s 
equation are of the type 

F= [APn (^) + BQn m [GPn iv) + PQn (v)]- *(197), 

where A^B, (7, D are arbitrary constants. 


* Phil. Tram. (1880), p, 138. 


B. 


7 
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Exanvples. 1. Prove that 

[(cosh a cos jS — cos y)* + sinh* a sin* ^ 

= J 2 (2w+l) Qn (cosh a) (cos P„ (cos y). 

0 

(C. Neumann.) 

2. Prove that a function B O) which satisfies the differential equation 
(192) and is ZiOro for /3=0, /S= 7 r is a solution of the homogeneous integral 
equation 

B(6)=fiaiB^d( R(ik)^0, 

J 0 

where n is determined by the condition that the integral equation should 
possess a continuous solution which is not identically zero. 

(M. Abraham.) 

3. If A (a) be defined by the equation 

A (a) = sinh* o J e ' ®®®'' “ ®“* ^ £ (/3) sin ^ rfjS, 

it satisfies the differential equation (192). A second solution of this 
equation is given by 

I(a)=sinh*a A(^)8inhfrf^ Riik)>0, 
and is suitable for the representation of divergent waves. 

(M. Abraham.) 


§ 34. Paraboloidal coordinates. 

If we write 


z + ip ,== (Oo + i^ofy «o^ = - a, = /3 
so that the transformation is 

z = -a-^, p = 2V^a;tf (198), 

the differential equation (143) becomes* 


9a> 0a ^0/3® ^"^'^^^0^8 0t® “ 


and is satisfied byf 


W = A (a) B (^) 


* Cf. H. J. Sharpe, Quarterly Journal^ Vol. 15 (1878) ; Proc. Camb. Phil, 
Soc, Vol. 10 (1899), p. 101; Vol. 13 (1905), p. 133; Vol. 15 (1909), p. 190; 
H. Lamb, Proc. London Math, Soc, Ser. 2, Vol. 4 (1907), p. 190. 

t The existence of elementary solutions for the paraboloid and certain 
other surfaces is established in BOcher’s Die Peiheimitwickelungen d^, Potential- 
theorie. See the table on pp. 266-7. 
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if 


a^+{m + \)^-{h-¥a)A=0] 


doL^ 




'“1 


.(199), 


where h is arbitrary. 

Putting 2ikn = ik {vi + 1) ^ h we find that the differential 
equations are satisfied by putting 

A = (2ika), B = (2ikff), 

where Fyn^{s) satisfies the differential equation* 

s'^^j,+(m+l-s)^ + nF=0 ( 200 ). 

When 71 is a positive integer, one solution of this equation 
is furnished by Sonin’s polynomial f Tm^ (s), which may be 
defined with the aid of the expansion 

J^_ 00 

(1 + = - r (?« + »i + 1) (s). . .(201). 

n=0 

A few properties of this function are given here for the sake 
of reference. 


(^) — tT 


r (m + 71 -h 1) \7i r (m 4- /i) I n — 1 1 1 

5^-2 


r (m -i- — 1) I — 2 1 2 


I e-*s”‘ r,„'‘ (s) r,/ (s) rfs = 0 


, v=7=n 


r (?i -f 1) r (m -f + 1 ) 

^,r„«(s)=2’,::;:(«) 


1 / = r? 




^[6-2 ’™«(s)] = s”'-"2’:-.(s) 


...( 202 ), 

...(203), 

...(204), 

...(205), 


* This is a slight modification of Weiler’s canonical form for an iquation of 
Laplace’s type, Crelle^s Journal^ Bd. 51 (1856), p. 105. The equation is discussed 
for real values of m and ii by 0. ^ch\ 6 mi\ch,H()here 7 i Analysis ^ Bd. 2 (1874), p. 517. 

t Math. Ann. Bd. 16. Further properties of the function are given by 
L. Gegenbauer, Wien. Ber. (1887), p. 274, who proves that the roots of the equation 
( 3 ) =:«C, •considered as an equation of the »th degree in 3 , are all real, 
positive and unequal. This is a generalisation of the result obtained by 
Laguerre i(5^ the case m = 0. A geometrical proof has been given by Bocher, 
Proc. of the Amer. Acad, of Arts a 7 id Sciences^ Vol. 40 (1904). 


7—2 



100 


SONIN’S POLYNOMIALS 


[CH. 


w- r t r(. 4. 


Equations (201), (202) and (203) were given by Abel* * * § and 
Murphy f for the case m = 0: the polynomial is then equivalent 
to the polynomial of TchebycheffJ; and Laguerre§ which occurs 
in the theory of interpolation and also in the theory of continued 
fractions. When m = ±^, the polynomial can be expressed in 
terms of the polynomial Un{oc) discussed by Tchebycheff|| and 
HermitelT, or in terms of the function of the parabolic cylinder, 
discussed by Weber**, Whittakerff and others;^;]:. 

The above analysis indicates the existence of a wave-function 
of the form 

n = Tm^ {2ika) (2ik/3) . . .(208). 


This function can be expressed as an integral of the form 
used in § 5, we have in fact the equation 

(kp)^ Tm^ {2ika) TrrT (2ik^) 

(~1)" 


27rr {m 4-ri 4- 1) 


l^2n 

)Jo ^ 


gkpe^Y rpn jr_ ^k (p COS 7 + izy\ e~^^ydy 

(209), 


from which the required representation can be immediately 
derived. In this formula m is either zero or a positive integer. 
The convergence of a series of terms of type (208) in which 


* Mimoires de mathimatiqxie par N. H. Abel, Paris (1826) ; Oeuvres^ Sylow 
and Lie, t. 2. 

t Camhr. Phil. Trans. (1833). 

t M4m. de VAcad. de St P^tersbourg (1860). 

§ Bull, de la Soc. math, de Fran^e^ t. 7 (1879) ; Oeuvres de Laguerre, t. 1, 
p. 428. 

II Loc.^cit. See also Sturm, Liouville^s Journal, Vol. 1. 

IT Comptes Bendus, t. 68 (1864), p. 93. The Hermite functions have been 
generalised by Curzon, Proc. Loudon Math. Soc. Vol. 13 (1914), p. 417. The 
generalised functions are intimately connected with the functions considered here. 

** Math. Ann. Bd. 1 (1869), p. 1. 

ft Proc. London Math. Soc. Ser. 1, Vol. 36 (1903), p. 417. 

Baer, Dies. Cilstrin (1883) ; Hantzschel, Zeitschr. filr Math. Bd 33 (1888) ; 
Adamoff, Annales de St Pitersbourg, t. 5 (1906) ; G. N. Watson, Proc. London 
Math. Soc. Ser. 2, Vol. 8 (1910), p. 393. 
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n takes different integral values can be partially discussed with 
the aid of the equation 

(ix) Tm” (- ix) 

1 fi . «(«t + M+l) 

r(jn-l)r(m+l)J L ^ l(OT+l) 2 (m+ 2 ) 

n (w — 1) (m + 71 + 1) (m -f + 2) 

1 . 2 (m + 1)^ (m 4- 2)^ (m -h 3) (m + 4) 

which shows that the modulus of Tm^(ix) increases with x. 
Hence if a series of terms of type (208) converges absolutely 
for any given value of a, it converges absolutely for all smaller 
values of a. 

For a fuller discussion of the convergence it would be useful 
to have an asymptotic expression for Tm^{s) when n is large. 
Suitable asymptotic expressions have already been found for 
the case m = ± ^ by Adarnoff and Watson. 

The differential equation (200) has been studied for general 
values of m and ti by Pochhammer*, Jacobstahlf, Whittaker t 
and Barnes§. It usually possesses two distinct solutions which 
can be expanded in power series converging for all finite values 
of s. If, however, m and n are positive integers, there is only 
one solution which can be represented by a convergent power 
series in s, the other may be defined by the equation 

Urn’' («) = Y (n+ 1 ) ^ . .( 211 ): 

it contains a logarithmic term. For negative integral values of 
n we may adopt the definition 

=1^5 (212). 

It should be noticed that when |a| is large, Um^{2ika) has 
an asymptotic expansion of which the first term is 

The solutions of type 

Qik iz+ct) ±m JJji {2ika) Trrt^ (2ik^) 

• « 

* Math. Ann. Vol. 36, p. 84 ; Vol. 46, p. 584. f Ihid. Vol. 56, p. 129. 

X Buu'^Amer. Math, Soc, (1904). 

§ Cambr. Phil. Tram, Vol. 20 (1906), p. 253. 




( 210 ). 
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are consequently suitable when a is large for the representation 
of waves diverging to infinity in the positive direction of the 
axis of z. 

It may be worth while to mention here that the functions 
(s) both satisfy Gegenbauer s difference equations 

( s ) = ( m + «) f:-^\ (s) + f:-;i (s) j 

n (to + n) FZ {s) - {s - (to + 2m - 1)} F^^ (s) + F"~^ (s) = 0 

(m - 1) F^ (s) = (6- - (TO + m - 1)1 F:;: (s) - Kt («) 

(m - 1) {s) = {s - (to + 1)1 F:i\ {s) - s (.9) 

(213). 

The function {s) also satisfies an equation analogous to 
(204). 

§ 35. Relations between different solutions. 

Many useful formulae may be obtained by expanding known 
wave-functions in series of elementary wave-functions of type 
(208) and by identifying our elementary wave-functions with 
certain definite integrals which are known to represent wave- 
functions. For instance, we have the expansion 1 tan ^ o) 1 < 1 

^iicz cos sin 

= tan ^ j sec* — S {— ly n I (m n) I 
X tan*^ ^ Trn^ {2ika) {2ik^) (214) 


which enables us to represent a plane wave with the aid of 
a double series of solutions of the form (208). 

Further identities may be obtained by deriving wave- 
functions from Cunningham’s solutions’*^ of the equations 


du _ d^u 
dr dx^ 0?/* 


du 

3 t dx^ ’ 

The first equation possesses the polynomial solutions 

^:h) 


^Tl 


H)- 


XT 


9t-l 


( 215 ). 


( 216 ) 


* Proc. Roy, Soc. Ser. A, Vol. 81 (1908), p. 310. See also Wera Lebedeff, 
DtM. Gottingen (1906) ; Math. Ann. (1907). The first result is given 8y Appell, 
LiouvilWi Journal^ Ser. 4, t. 8 fl892), p. 187. 



RELATIONS BETWEEN DIFFERENT SOLUTIONS 


and also the solutions 
_ ^ 


n+2 
XT 2 


...(217). 


" -UW’ 

The second equation possesses the polynomial solutions of 
typo 

T” p™ 2 to” sin TO (^ - ^o) (218) 


and also the solutions of type 


e'*" Tto” sin m(^ - .,,(219). 

Wave-functions may be derived from these solutions by the 
method of § 13. 

We add here a few relations which are obtained by ex- 
pressing the solutions thus formed in terms of old solutions. 

^ _ 2 cos w] 


00 Vt 

= S (- 1)” r (m + n + 1) p Tto" (?) Tto" {y)^. . .(220), 

jrTO'‘(s) (221), 

• m-—n 

= Wil 

2»1'+1X2»4»«^^ (Xp) 

= r (n + 1 ) r (»! + m + 1 ) p’« Tto" (£) dr 

(223). 

The proofs of these are left to the reader. 

Prof. G. D. Birkhoff has remarked to me that the differential 
equation (200) can be regarded as a limiting case of the hyper- 
geometric equation when two of the singularities coincide at 
infinity*, consequently many properties of the solutions can 
be derived from known properties of hypergeometric functionsf. 
It shquld be noticed that when W is independent of t there 

* Cf. B6cher, Die JReihenentivickelungen der Poteritialtheorie, p. 137. 
t This^ethod was used in a particular case by Kummer, Crelle^s Journal^ 
Bd. 15 (1836), p. 138. 
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are elementary solutions of equation (143) of the form W = AB 
where 

’ (2iv'A/3)™ 

We thus obtain elementary solutions of Laplace’s equation 
of the form 

Jm (X Va) /„ ( XV/S) cos m ((^ - <^o) (224), 

where m, </>o are arbitrary parameters. 


§ 36. Toroidal coordinates. 

If we put 

a) — pcos(j), 2 / = psin<^, , 2 r= f cosh f sinh ci)...(225), 

. . .a* — iylr 

p + ^5'= a coth — 2 “ > 

^ sinh <7 a sin yjr /226) 

^ cosh cr — cosi/r ^ ^ COsh o- — cos 

the wave-equation becomes 

d ( sinh o- sin yfr du) d ( sinh cr sin yfr du ) 
do- |(cosh cr — cos yp-y do] dyp ((cosh o — cos ypy dyp) 

sin i/r d^u sinh o d^u _ ^ 

sinh cr (cosh o — cos ypy d<f>^ sin yp (cosh o — cos ypy dco* 

(227). 

This is satisfied by 

u = F(o) 0 (yp) (cosh cr — cos yp) cos m{<^ — (f>o ). . .(228) 


if 


c<«ch,r (sinh -{«<»+ 1) + ^-^} F. 0 

coseo t ^ (sin + ||) + {„ („+ 1) - F-o' 

(229). 


Hence we obtain wave-functions of the form 

(cosh o — cos yp) (cosh < 7 ) (cos yp) cos m{<p’- (/>«) 

(230). 

Other solutions of the wave-equation may be obtained by 
replacing the functions P^"", P, * by Q/. 

Many useful formulae may be obtained by e'gpanding 
particular wave-functions in series of wave-functions of type 
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(230). The expansion of unity, for instance, gives rise to 
Neumann’s expansion 

co.h.Lo..>. -.!/"'•+') «" f) ...(asi). 

It should be noticed that when we make the substitution 
(225) the wave-equation becomes 


d^y 1 du 
dp^ p dp 


^ ^ 0...(232) • 


it thus possesses elementary solutions of the forms 

Kp (Xf ) (kp) cos m (<^ - (^o) (233), 

Jp (^X^) J,n (kp) cos m ((^ - ^o) (234). 

The expression of solutions of type (230) in terms of the 
solutions just found leads to some interesting identities. Thus 
we have the equation 

Kp (\0 {\p) J„ {\a) XP+»‘-"+‘ d\ 

I r + m + 1) r (p + 1) r (m + 1) 

F(n:+iT'* 


_ ^n—p—7n—2 , 


X (cosh <r - cos -f (cos f) (cosh <r) 

p> — l, m> — ly 1 (235). 

Many important formulae connected with Bessel functions 
are simply particular cases of this one*. It should be remem- 
bered that 

^4 (^) = \/£ 4 

The corresponding integral in which Kp (Xf) is replaced by 
J-p (X^) can also be evaluated in terms of Legendre functions, 
but the formulae are more complicated. The case p = — m is 
discussed by Macdonald i*. 

It should be noticed that if we write 


cosh (a — ft)) = i cot yfr, 
sinh (a — ft)) = i cosec yf/, 


cos (i3 — (f>) = coth cr 
sin (j3 — (f)) = ±1 cosech <7 



three relations of 


type 

0(07, y) 


■ ^ 8 («, /Q) 

c d(s, t) 


,(238) 


* See, for instance, the formulae given by H. M. Macdonald, Proc, London 
Math. Soc. li^r. 2, Vol. 7, p. 147, and by the author, ibid. Vol. 12, Abitracts. 
t Loc. cit. p. 142. 
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are satisfied and so the functions a, /9 can be used to obtain 
an electromagnetic field by the method of § 5. It is easy 
to verify that the function 

u — (cosh a — cos /8) (239) 

satisfies the wave-equation, / being an arbitrary function. 

We add here a few formulae for (cosh <t), Qn^ (cosh a) ; 
these and other formulae will be found in the memoirs of 
Dr Hobson and Dr Barnes to which we have already referred. 

Pn^ (cosh a) = -ptt- \ coth^ ^ 

1(1 — 7 ) 1 ) 2 


xPj— n, n+l; 1 — m; — sinh^ - 

027n 

^ n — p-2<r\-m ^-(n+Do- 

X P - m, 1 4- n — m ; 1 — 2m ; 1 - 

cr > 0, 


QrT (cosh 0 -) = (- 1)™ V-n- ^ (1 ~ 

V 1 (?i 4- 

X P {m 4- n 4- m 4- 1 ; ^ 4- f , o’ > 0. 

Various asymptotic expansions for these functions are given 
by the authors jusfc named and by Dr Nicholson*. 

It should be mentioned that the solutions of the wave- 
equation that have just been obtained are not directly useful 
for the treatment of the boundary problems of mathematical 
physics. They may, however, be used to construct useful 
solutions of the equation Au 4* k^u = 0 by means of various 
artifices. If, for instance, we multiply one of our wave-func- 
tions by and integrate with regard to t between ^ and x) , 
the resulting function will often be a solution of Av 4- k^v 0. 
This may be illustrated by taking the wave-function 
u = Jo (ixl;) Jrn (xp) cos m (<^ - (f>o) 
and using the formula 

Too 

Jo (X Vc ‘^^2 _ ^‘ 2 ) QXkct ^ Jc 2 ^ 

z CVA:- — X^ 




^ > 0 
X^>¥ 

..r.(240). 


British Association Reports, Winnipeg (1909), p. 391. 



VI] 


TOROIDAL COORDINATES 


107 


The integration can be taken between other limits in 
certain cases; for instance, the integral 


z‘^+ ip+a)'^ 



J Kp + - (p - 


represents the solution of Av -\-k^v = 0 corresponding to a 
circular ring of sources. In this case our wave- function u 
is a constant multiple of cosh a — cos yjr. 

The theory of electrical oscillations on a conducting anchor 
ring has been treated by H. C. Pocklington*, W. McF. Orrf 
and Lord Rayleigh without the use of toroidal coordinates ; 
the results are of course only approximate. 

§ 37. Solutions of Laplace's equation. 

If we put 

_ a sinh a _ a sin i/r / 9 A 9 \ 

^ cosher— cos cosh cr — cos * 

the angle -i/r may be interpreted as the angle which two fixed 
points A, B whose coordinates are z = p = ± a, subtend at 

PA 

a point P (p,z); the quantity o* may be interpreted as log . 

The surfoces yfr = const, are spheres having a real circle (p = a, 
-3^ = 0) in common, the surfaces a = const, are anchor rings. 

If we use the toroidal coordinates a, yfr, Laplace’s equation 
becomes§ 


da- (cosh cr — cos ^jr do-} d-^ (cosh a — yfr j 

1 dhi „ 

^ ^ — Q 

sinh a (cosh cr — cos 'v/r) d4>^ 


^ Proc. Comb. Phil. Soc. Vol. 9 (1897), p. 324. 

+ Phil. Mag. Vol. 6 (1903), p. 667. 

X Proc. Roy. Soc. Ser. A, Vol. 87 (1912), p. 93. See also C. W. Oseen, 
Phys. Zeitschr., Dec. 1st (1913) ; Arhiv for Mat. A,st. och Fysik, Bd. 9 (1913). 

§ B. Riemann, PartielleDiferentialgleichimgen, Hattendorf’s edition (1861) ; 
C. TSeummini^ Theorie der Elektricitdts- tind Wdrme-Vertheilung in einem Ringe, 
Halle (1864); W. M. Hicks, Phil. Trans. (1881), p. 609; A. B. Basset, Amer. 
Joum. Vol. 1?^, Hydrodynamics^ Vol. 2. For an alternative method see F. H. 
Safford, Annals of Mathematics^ Vol. 12 (1898), p. 27. 
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and possesses solutions of the form 

P^_j(cosh <r) 

u = (cosh a — cos yjr)^ cos n (^|r — yfto) COS m (<^ — (f>o) 

Q^_j(cosh cr) 
(243), 

which are suitable for the treatment of problems connected 
with the anchor ring, circular disc and spherical bowl 

For problems connected with two spheres bipolar coordinates 
may be used ; the appropriate substitution isf 
asinyjr ^ asinhcr 

^ cosh a — cos yfr * cosh a — cos yfr ’ ' 

The surfaces a ~ const, are now coaxal spheres with imaginary 
common circle. The radius of the sphere cr = o-q is a | cosech ctq | 
and the distance of its centre from the origin is a | coth ctq | . 

The ratio of the distances of a point from the limiting points 
of the system of coaxal spheres is and the angle between the 
radii from these points is 

The appropriate solutions of Laplace's equation are now of 
the typej 

u = (cosh cr — cos [A cosh (yi -f cr + P sinh (n + cr] 

X cos <jf>o) [fPn^ (cos y[r) + (cos ^/r)]. . .(245). 

It should be noticed that when we are using toroidal 
coordinates the function 


‘/[< 


u = (cosh <r — cos ^[r) 1 </> ± i log tanh — J cos ^ ('^ — y/ro) 

(246) 

satisfies Laplace s equation and that when we use bipolar 
coordinates the corresponding solution is 

(cosh cr — cos '^)^/ j^</) ± i log tan cosh ^ (o* — cto). . .(247). 


« 

* See for instance E. W. Hobson, Camhr. Phil. Trans, Vol. 18 (1899) ; C. 
W. Oseen, Arkiv for Matematik, Bd. 2, No. 5 ; H. C. Pocklington, Phil. Trans, 
A, Vol. 186 (1895), p. 603. 

t W. Thomson (Lord Kelvin), Liouville^s Journal (1847). 

X O. B. Jeffery, Proc. Roy. Soc. Ser. A, Vol. 87 (1912) ; G. R<«De*.n, Phys. 
Revieic (1912) ; G. Darboux, Bull, des Sciences math. t. 31 (1907), p. 17. 
Another method of dealing with problems connected with t'.f5 spheres is 
described by A. Guillet and M. Aubert, Journal de Physique^ t. 3 (1913). 
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It should be mentioned here that other simple ^solutions 
may be obtained by using the formulae 

P\ 1 (cosh (t) - - — A / — — cosh na. 

^ V TTSinho- 

EXAMPLES. 

1. If with the notation of § 36 we write 

f -f tp = a cosh (a + ip) 
the wave-equation becomes 

2 coth 2a|^-f-2 cot 2/3 (cosech^a+cose^^)^^^ 

^ 2 ^, '»«v 

-H (sech2 a - sec^/S) g-- = 0. 

Hence show that there are wave-functions of the form 
u — A(a)B O) cos m{(t)- (/iq), 

where a, ky m and 0o a-re arbitrary constants. 

2. Prove that if p+iz=f{a-\-i^) the wave-equation becomes 

0% 02^ 1 /0p ^ 1 / 1 02i^ 1 dH\ _ 

0a2 0/32 /j \0a 0a 0/3 0/3/ J^yp^dcf)'^ dt‘^ ) ’ 

and obtain elementary solutions of type A (a) B (/S) when 

gq-iprsa cosh (a-hi/3). Notice that the solutions of equation (188) are not 
wave-functions, they are analogous to the stream-line functions of hydro- 
dynamics. 
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§ 38. The method of Stieltjes*. 

Wave-functions which are homogeneous functions of x,y,z,t 
may be studied with the aid of the substitution 


x^s cos 6 cos s sin 6 cos % 

y — s cos 6 sin <^, ict = s sin 0 sin x 

The wave-equation in these coordinates has the form 



(248). 


d^u Sdu Idht 1 d^u 

ds^ s ds ^ 30- cos^ 6 d(f)^ 

1 cot 0 — tan 6 du 


= 0. 


6‘2 dO 

Putting cos 20 = fjL, we find that there are elementary solutions 
of degree 2n of the form 

(249), 

if 


d 

d/jb 




2 ( 1 +^) 


^'-1 = 0 

2(1-/^)] 
(250). 

This equation is satisfied by 

(1+ (1 - rtf H- 1+ ‘ i-") 

(251), 

with the usual notation of the hypergeometric function!. 


* Coiuptes Reiidiiis, t. 95 (1882), p. 901 ; Liouinlle'» Journal, Ser. 4, t. 5 
(1889), p. 55. See also Tisserand, TraiU de mecaniqae celeste, Paris (1889) ; 
H. Bateman, Proc. London Math. Soc. Ser. 2, Vol. 3 (1905), p. Ill 

t This gives a polynomial if either m + 1 + ^ - or - - n is zero or a 

negative integer. 
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It should be noticed that if we write 
cos 6 = sec a, sin 6 — i tan a, sin % = cosh w, cos — % sinh Wy 
fjL = sec^ a + tan^ a, 
we obtain real wave-functions of the form 

tt = ^ e'~^ @ ( M ) cos m (<^ — <^o) (252), 

the variable in the hypergeometric function is now — tan^ a. 


When m=^p the equation for @ is the equation satisfied by the 
associated Legendre functions. We thus obtain wave-functions 
of the form 

u = (cos 26) + x) (253). 

Comparing this with the elementary solution of Laplace's equa- 
tion in polar coordinates (/’, 6, </>), we see that if /(r, 6y (f)) is a 
solution of Laplace’s equation f (s^, 2^, -f %) is a wave-function. 
We may thus derive wave-functions from harmonic functions; 

1* 

in particular, the fundamental harmonic function — gives us the 
fundamental wave-function ~ = r:; • We have 

5 - + 2/2 -f 2:2 - ^2 ^2 

already remarked in § 13 that Lord Kelvin’s method of inver- 
sion may be extended to wave-functions, it is easy to see that 
the result is an immediate consequence of the fact that the 
differential equation (250) is unaltered A^heii — (?i + l) is 
written in place of 

It is easy to see that there are (n + 1)^ linearly independent 
polynomial solutions of degree n, for a general polynomial of 

degree n contains — (/i + 1) (fi -f- 2) (/i -f 3) coefficients and when 

this is operated on with Q the vanishing of the resulting 

polynomial of degree n — 2 gives ^(n— l)n(n-f-l) cdhditions. 

The difference between these two numbers is (n -4- 1)^. 

A polynomial solution of degree n is given by the integral 

^ ~ a + y sin a -f iz)^ {x sin a — y cos a — da, 

A set of {n + 1)^ linearly independent polynomials is obtained 
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by allowing m and p to take the values (0, 1 ... n). A better 
set of solutions is obtained by using the integrals of type 

f [(a? — iy) i (-3^ — ct)]^ 

Jo 

X [{x -1- iy) -f i 4- doi. 

Polynomial solutions may also be obtained by differentiating 

the fundamental wave-function ^ and using generalised inver- 

s 

sion*. The polynomial solutions were first discussed by 
Cayleyf. WaelschJ has recently studied them from a new 
point of view. 

Example. Prove that when is a positive integer 
1 3 2% — 1 1 r”" 

k k ~ ^ J 0 0 ^ ^ ^ cos d(f) d^. 


§ 39. The method of Green§. 

Homogeneous solutions may also be investigated with the 
aid of Green’s substitution 


a; = 8 sin a sin /3 cos <^, 
^ = 5 sin a cos yS, 


y = s sin a sin ^ sin <f>] ,^54) 

ict = s cos a 1 ’ 


The wave-equation now becomes 


d^u S du 2 du 1 d^u 
ds^ s 9^ 9a 9a‘^ 


1 I 1 -PO 

sin^ a sin 0 9y8 \ 9/9/ sin* a sin* ^ 90* ' * ' 


and possesses elementary solutions of the form 
u = s^A (a) B (/3) cos m (0 — 0o), 


* Cf. F. Didon, Annales de V:^cole Normale (1), t. 5 (1868), p. 229; t. 6 
(1869), p. 7; t. 7 (1870), pp. 89, 247 ; P. Appell, Rend. Palermo, t. 36 (1913) ; 
K. de F4riet, Comptes Rendus, Nov. 17th (1913). 

t Lio-i&ilU^s Journal, t. 13 (1848) ; Phil. Trans. Vol. (165) ii. (1875), p. 675. 
See also Hermite, Oeuvres, t. 2. 

Deutsche Math. Vercin, Bd. 19, p. 90. 

§ Carribr. Phil. Trans. Vol. 5 (1835), p. 395 ; Collected Papers, p. 187 ; 
Cayley, loc. cit. See also Heine, Handbuch der Kugelfunktionen, Bd. 1, p. 449 ; 
Crelle, Bd. 60, 61, 62 (1862 — 1863) ; E. W. Hobson, Proc. Londoh Math. Soc. 
Ser. 1, Vol. 24, p. 67, Vol. 25 ; F. G. Mehler, Progr. Danzig (1864) ; Crelle, 
Bd. 66 (1866), p. 161; C. Neumann, Zeitschr. Math. Phys. jod. 12 (1867), 
p. 116 ; V. Giulotto, Gior. d. Mat. 39 (1901), p. 162. 
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1 ^ 
sin/9 rf/S 
d^A 


■ 


doC^ 

We may thus take 

5 = Cl (cos - C2 (cos 
A = Vcosec a |^6i + ^ (cos a) +62 + 

where Ci, Cj, 61, 62 are arbitrary constants. 


....(256), 


§ 40. Wave-functions of degree zero. 

If fl is a wave-function of degree — 2, the formulae 




an an 


a.c 


dy * 

an 


.-r an 

Aa; — — - 

c a^ 


■ ct 


jr dil 


„ en en 


dx * 


Ey=y'^-2-+ct 

y c dt 

Ei = -^ + ct - 
c dt oz 


m 

dx 

djl 

dil 


-...(257) 


give a solution of Maxwell’s equations. 

A homogeneous wave-function of degree — 2 can, of course, 
be derived from a homogeneous wave-function of degree zero. 
If F is an arbitrary function of two variables subject to suitable 
restrictions, the integral ^ 

r.r sin a — y cos a — ( 


fin 

n= F 


cos a — ct 

7 -^>« 

nn a-\-iz 


doL (258) 


! X cos a y sin t 

represents a wave-function of degree zero, and when this is 

multiplied hy \ sl wave-function of degree — 2 is obtained. We 
s 

add here a few particular wave-functions of degree zero : 


X i iy 


). 


tan 


-1^ 


tanh--, log^^— ^ r 259 ). 


\z ±ct! ' X 

Electromagnetic fields which are derived from this type of wave- 
function of degree - 2 may be generalised by writing x-^ (t), 
y—riirf, ^(t), t-T instead of x, y, z, t respectively and 
integrating jround a closed contour in the complex r plane. 
The integrals thus obtained can generally be evaluated by 

B. ® 
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means of Cauchy’s theorem. Many of the results given in the 
next chapter are suggested at once by this method and may 
be thoroughly established by a method of direct verification. 

It is worthy of note that if we write r in place of ct in a 
wave-function of degree zero the resulting function is a solution 
of Laplace’s equation Ai6 = 0. A general solution of Laplace’s 
equation of degree zero can be derived at once in this way from 
the first of the solutions (259). We thus obtain Donkin’s 
formula* 



A similar result is that if 


n {x, y, z, Wy t) 

is a homogeneous function of degree ^ satisfying the equation 

dx^ dy^ dz^ ’ 

and 8 be written in place of Wy the resulting function is a 
wave-function. Now if f(XyyyZ) is a solution of Laplace’s 
equation, the function 

'^w — cz \w — ct' w — ct' w — ctj 


satisfies the requirements, consequently we may conclude that 
the function 




— h f ■ 

V 5 - CC \s-ct' 


» y ^ \ 

s — ct' s — ct) 


(261) 


is a wave-functionf. Other wave-functions may be derived 
from this by generalised inversion or by interchanging the 
variables Xy y, Zy ict 


* Phi^. Trans. (1857). This solution may be obtained at once from Jacobi’s 
theorem that if p, q, r are three functions of u which satisfy the equation 
p® + -f r2 = 0 and u is defined by the equation au = x'p{u) + yq (u) + zr (w), 
then an arbitrary function of w is a solution of Laplace’s equation, Werkcy 
Bd. 2, p. 208. See also Forsyth, Mess, of Math. (1898). 

t This result is obtained in another way by Pockels, Vhev dij partielle 
Differentialgleichung Au + k'^u = 0. Teubner, Leipzig (1891). 
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ELECTROMAGNETIC FIELDS WITH MOVING 
SINGULARITIES 


§ 41. An electromagnetic field with a simple singularity 
or electron, first model of a corpuscle*. 

We shall now derive a family of wave-functions from the 
fundamental wave-function 1/5^, where 

52 = [a; - ^ (t)] 2 -h [ 2 / - 7 ; (t)P + (t)P - c"" - Tp. . .(262) 

and T is a variable parameter, which is at first independent of 
Xj y, z, t Using a method invented by Prof. A. W. Conway f 
we consider the integral 



taken round a closed contour in the plane of the complex vari- 
able T. If this contour contains only one root r of the equation 
.9^ = 0, the value of the integral is ^ 




where 


V = (T)(a; - + T)' (r) (y - tj) + (r) (z-^)- c® (t - t) 

(263) 


and T is the root in question. 


* I have ventured to use Johnstone Stoney’s term “ electron ” to denote the 
simple point singularity and Sir Joseph Thomson’s term “ corpuscle ” to denote 
the elementary charged particle which has been discovered by experimental 
work. * • 

t Proc. London Math, Soc. Ser. 2, Vol. 1 (1903). Integrals over complex 
paths had beei-*employed previously in electromagnetic theory by Sommerfeld 
a.nd other writers. 


8--2 
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This function v vanishes when = y = < = t and 

so the wave-function 11 has a singularity which moves along the 
curve r represented by 

y==v{T), ^ = (204). 

If, moreover, the velocity of this singularity E is always less 
than the velocity of light, it is easy to see that v does not vanish 
for any real values of {x^ y, 5 , t) other than those just mentioned. 

When the velocity of the singularity E is always less than c, 
there is only one value of t less than t for which the equation 
52 = 0 is satisfied : x, y, z, t being supposed to be given. 

To prove this we surround each point E on F by a sphere of 
radius c (i — t) having E as centre ; then it is clear that each 
sphere lies entirely within the neighbouring one corresponding 
to a smaller value of t, provided T<t and + drf + d^‘^<c^dT^. 

This shows that one and only one of these spheres passes 
through a given point of space and so there is only one value of 
T< t for whieh the equation 

[« - f (t)J* + [y - ■»? (t)P + [z-^ = c^(t- t)®. . .(265) 

is satisfied*. 

Now let a point Q (Xy y, Zy t) move with a velocity less than c 
along a curve G and let us consider the variation of r with L 
As t increases from t tot + dt the radius of the sphere associated 
with each point E will increase by cdt and since Q moves a 
distance less than cdt in the interval dt, its new position will 
lie within the new sphere associated with the time t. Conse- 
quently the new position of Q lies on a sphere associated with a 
greater time r. 

Hence if Q moves in any manner with a velocity less than the 
velocity of lights r increases with t. 

Things are quite different when the velocity of E is greater 
than c. The spheres then have a real envelope and there may 
be more than one sphere through a given point in space, also r 
may sometimes decrease when t increases. 

In this case, however, v vanishes for values of Xy y, Zy t other 
than x — ^y y = rjf 2 — ^ = t, and so H has 00 ^ singular lines 

* Cf. A. W. Conway, loc. dt, ; H. Bateman, Manchester Memoirs (1910) ; 
G. A. Schott, Electromagnetic Radiation (1912). The theorem is due to Li^nard. 
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through each point E. These singular lines form a right 
circular cone whose axis is along direction of motion : each 
singular line is described by a singular point that travels with 
the velocity of light. 

When the velocity of the singularity is less than c we can 
obtain a solution of Maxwells equations having the moving 
singularity by using the potentials* 




^r(T) 

47ri/ 


_ eV(T) 
47ri/ ’ 


4>Trv ' 


4> = 


ec 


It is easy to verify that they satisfy the relation 

div ^ ^ = 0 

c dt 


Asirv 

.(26G). 

.(267). 


When the electric and magnetic forces are calculated from 
these potentials with the aid of the formulae 


H = rot A, 
it is found thatf 


E=^ 




H.. 


e 9 (t, a) 
4f7r 9 (y, z) ' 




e 9 (t, cr) 
ilTC d (x, t) 


.(268) 

(269), 


where 


<rv = r' i<»-^) + v"(y-v)+ K" (2 - ?) - (r^ + v'^ + n + C* 

' (270). 

It is clear from these equations that the magnetic force is 


* A. Li(5nard, Uedairage ^lectrique, Vol. 16 (1898), pp. 5, 53, 106. See also 
E. Wiechert, Arch. nierlamlaUes (2), Vol. 5 (1900), p. 54; K. Schwarz schild, 
Q'dtt. Nachr. (1903). The potentials are usually written in the form 
e , e[v] 


<!> = . 


47r 




, where the square bracket indicates 


r '• 

that the quantity enclosed is to be calculated at time T = t--. Cf. H. A. 

Lorentz, The Theortj of Electrons, p. 50. To obtain a model of a corpuscle we 
must write de instead of e and integrate over a small region. 

t These expressions for the components of E and H were communicated to 
me by 'Mr R. Hargreaves in 1909 ; they should be of some historical interest 
in connection with the general theory of § 5. This was, however, the outcome 
of some inaJl^endent work. Cf. Proc. London Math. Soc. (2), Vol. 10 (1911), 
p. 96. 
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perpendicular to the electric force and also to the radius 
from the effective position of E\ for we have the relations 


t 

dr 

v^=y-v 

1 

II 



( 271 ). 


It also follows from these relations that r satisfies the char- 
acteristic equation 



This is to be expected because, as Jacobi has remarked* for the 
case of Laplace’s equation, the argument t of an arbitrary func- 
tion occurring in the solution of a partial differential equation 
must satisfy the partial differential equation of the character- 
istics f. 

To prove that there is a constant charge e associated with 
the singularity of our electromagnetic field we shall calculate 
the integral of the radial component of E over a sphere having 
the singularity as centre. We have to evaluate the integral 

ce f f /da- dr da- dr da- dr 1 da- St\ , ^ 

4}7r JJ V9a? do! dy dy ^ dz dz dt dt) ’ 


which is easily transformed into 


47r JJ 


dS, 


Transforming the axes so that the axis of z is in the direction 
of motion of the singularity, we may put 

pz=r{v cos 0 — c), dS = r* sin 6 d0d(f)y f -h 4- f 
and our integral becomes 

ce /*"■ (c* — v^) sin 6 d0d<f> _ 

^ttJoJo (vcoQ0 — cy ~ 


* WerkCy Bd. 2, p. 208 ; Crelle^s Journaly Bd, 36 (1848). 
t For the general theory of characteristics see Hadamard, Propagation des 
Ondes (1903), Chapters vn. and vni. ; J. Coulon, Comptes Rendusy t. 128 
(1899), p. 1386 ; A. V. Backlund, Math, Ann, Bd. 13 (1878), p. 41^; J. Beudon, 
Comptes RenduSy t. 124 (1897), p. 124. 
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Prof. E. T. Whittaker* has calculated potentials F, 11 from 
which A and can be derived by using (7). He finds that 

n-(0,0,Sf), r=(0,0,i\r), 

where 

4f7rS = e 8inh~^ , 477^ = — e tan”^ - — ^ 

47rK = — 6 log V(a^ ~ ^)- + (j/ — 

(273). 

It may be verified without difficulty that the functions are 
wave-functions. This result is a particular case of the 
following general theorem. 

Iff{x, y, z) is a homogeneous function of degree zero satisfy- 
ing Laplaces equation Att = 0, the function 

n Hr), y-^(T),^~r(T)] (274) 

is a wave- function. 


* 

§ 42. The electromagnetic field due to a moving doublet. 

Let us now derive an electromagnetic field by superposing 
two electromagnetic fields of the type just described wherein 
the singularities move along the two neighbouring curves 

y = v(r), •s = ?(t), 

^=?(n)+ e«(Ti), y = '»7(T,) + e/3(Ti), ^ = ?(t,) + 67 (t,), 

e being a quantity whose square may be neglected. 

If Ti is defined in terms of x, y, ,3^, t by the equation 

[x-^ ( ti ) - ea + [y-v ( ti ) - 6/3 ( t ,)]“ 

+ \z-K (Ti) - 67 = C= (< - Ti)* 

and Tj = T + €0, we easily find that 

j/(9 + a (a; - I) + ^ (y - 7?) + 7 (^ - ?) = 0. 

Also if Vx is the quantity corresponding to v, we hav% 

= 1/ + 6 [Ov<T + a' (a? — ^) -f- 0' {y —v)4-y (z — ^ 

— — — 7 ?'] = 1 / + € [0v(T 4- p], say. 


No\/ if 


A '=- 


e ire'll) + €a'(Ti)} 


47n/i ’ 


* Proc» London Math, Soc, Ser. 2, Vol. 1. 
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we find that 




■ [vcl' ( t ) + pd^"(T)-p0<r^\T)l 


ec 


47ri^^ 


[p + vQ(i\ 


POL 


But 

p^' + P0^'' - P0a^' = n' ( 2 / - - 7 n' - O ~ 

+ c^OL — r)'n- <T {a{t — t) — n{y — 7)) + m {z — f)}, 

where ^ = / 85 '' — 7 ^', — cl ^\ n = a?;' — /3f'. 

Hence we may write 


<#>=- 


47r 

ec 

47r 


L9y 

d 

dx 


© 

(“)■ 


dz V 


0 y 


0 /a\ 
i^y 0 ^ \z^/ 

©' 


...(275). 


The electromagnetic field derived from these potentials is due 
to a moving electric doublet. It should be noticed that we 
have the relations 


.(276). 


Za-f my8 + n7 = 0] 
nirj' n^' 0\ 

We can write down by analogy the potentials for an electro- 
magnetic field due to a moving magnetic doublet. They are 


a^ = 


<^ = - 


47r 

ac 

_u - 

47r 


i/ 

Idy \ 

\d_ 

dx 


:© 

©■ 


where 


d 

dz 

dy 


L 

z \ v , 

U y 


©■ 


dt 
9^ V I' 


....(277), 


So + ^070 = 0 ] 

,v'4-^or=o] 


.(278), 


Zoao+mo/3o 

-f mo^ 

and Ofo; 7 o> ^o> ^ 0 , are functions of r. When a, / 8 , 7 , Z, m, 
n are functions of r which are not connected by the relations 
(276) the potentials (275) can be used to construct an electro- 
magnetic field which must be regarded as that due to an 
electric doublet and magnetic doublet which move together. 


§ 43. Electromagnetic fields in which singularities are 
projected from a moving point or curve and travel with the 
velocity of light. • * 

We shall now develop some mathematical analy^s of con- 
siderable interest whose physical significance is not yet fully 
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understood. At first sight it seems appropriate for a discussion 
of an emission theory of light in which waves in the aether are 
either produced or guided by small particles which move in 
straight lines with the velocity c. After further study I have 
thought it may be useful for a discussion of the question Has 
the aether a structure ? 

This question has already been raised by Sir Joseph 
Larmor* and Sir Joseph Thomsonf. The latter has, indeed, 
developed a theory in which the aether has a kind of atomic 
structure of which the elements are Faraday tubes};. In the 
most recent form of the theory it is assumed that the electric 
and magnetic forces are zero outside the tubes and that a 
certain amount of work is performed when one corpuscle 
crosses a tube of force attached to another. In an application 
of the present analysis to Sir Joseph Thomson’s theory the 
aim would be to build up his discontinuous electromagnetic 
fields from electromagnetic fields with certain types of sin- 
gularities, making use of discontinuous definite integrals. To 
illustrate the possibility of doing this it will be sufficient to 
mention the definite integral 

jj. r da dal 47r ^ 

Jo [.-S' + liTcosa -h zy sm a rj r 

= 0 ^ < 0 . 

The electrostatic field derived from the function V is zero on 
one side of the plane ^ = 0 and has the character of the field 
due to a point charge on the other side. It should be noticed 
that the integrand is a potential function which becomes infinite 
along the line z — 0, oj cos « + 3/ sin a = 0, and as a varies this 
line sweeps out the plane of discontinuity of our electrostatic 
field. 

To generalise this result we must endeavour to •solve the 

* Aether and Matter (1900), p. 188. The question as to whether the aether 
is continuous or discontinuous is discussed by H. Witte, Ann, d, Phys. (4), 
Bd. 26 (1908). 

t Presidential Address, British Association Bcports, Winnipeg (1909). 

X Becent Researches on Electricity and Magnetism', Electricity and Matter', 
Phil. Mag, V%. 19 (1910), p. 301, Oct. — Dec. (1912). See alsoN. R. Campbell, 
The New Quarterly (1909); Phil. Mag, Vol. 19 (1910), p. 181. 
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general problem of finding electromagnetic fields whose 
singularities lie on moving curves*. 

A partial solution of this problem may be obtained by con- 
sidering first of all the field represented by equations (10) 
of § 5. We may obtain a suitable pair of functions a, ^ by 
solving the equations 


'0a /Bay ^ 

,dx] \0;y j ^ \dz) & V9^/ 



0a 0/3 0a 0^ 0a 0/3 _ 1 da d0 
dx dx ^ dy dy dz dz dt di 


for clearly 


'0a 0y3 da d0\^ 

^y dz dz dy) 

K^y) \dz} J LV0y/ \dz) J \0t/ dy dz dz ) 

\0i/ \dx) J ic^\ dt ) \dx) J \c^ dt dt dx dx) 

— — ^ /0a 0y8 ^ 0a 0^Y 
\0iZ7 dt dt dx) ‘ 


Two other equations can be obtained in a similar way and so it 
follows that if we make a suitable choice of an ambiguous sign 
which is involved V in the definitions of the functions a and 
the equations (10) will be a consequence of equations (279). 

A more general electromagnetic field is obtained by 
multiplying the components of M in equations (10) by an 
arbitrary function /(a, /3). Since the components of M are 
necessarily solutions of the wave-equation it follows that, if 
g = l/f, an expression of the type 

0(a,^) 

9 (y. ' 

is a solution of the wave-equation (8). 




The aim must be not only to obtain a complete generalisation o^ Green’s 
equivalent layer which will be applicable to the case of a moving surface, but to 
obtain, if possible, an analysis of the electric and magnetic currer# sheets which 
are required. Of. p. 29. 
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The electromagnetic field which has just been obtained 
generally has singularities at space-time points for which 
/(a, /3) is zero. Let us write f=:u{x, y, z, i) -f iv (x, y, z, t) where 
u and V are real when x, y, z, t are real ; then the points for 
which /= 0 lie on the moving curve, defined by the equations 
u =0, ?; = 0. Now it follows from (279) that / is a solution of 
the equation 



and consequently 



du dv dii dv dudv ^ 1 du dv 
dx dx dy dy^ dz dz & di * 


Now let F{u,v) = 0 be the equation of a moving surface which 
always contains the moving curve, then if 


dFdv 

dFdu 

dFdv 

dF du 

du dx 

dv dx * 

^ du dy 

dv dy ' 

dF dv 

dFdu 

__ dF dv 

dFdu 

du dz 

dv dz ' 

^ du dt 

■ dv U ’ 


we have 



consequently 



This means that the component velocity of the surface .F = 0 in 
the direction of the normal at {x, y, z, t) is less than the velocity 
of light, it is equal to the velocity of light only in an excep- 
tional cas:^ Since this is true for any surface that always 
contains the moving curve it follows that the curve can be 
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regarded as moving with a velocity less than that of light. It 
should be understood that the curve generally changes in shape 
as it moves but this is not necessarily the case. 

It often happens that the moving curve / = 0 can be 
regarded as made up of the instantaneous positions of a series 
of points which move in straight lines with the velocity of 
light. To see this let us suppose that a function <f) (a, 0) can 
be found such that is a real function of Xy y, z, t It is evident 
from (279) that 

dx ~dx ^ cy dy^ dz dz & dt dt ’ 
and this means that if a point starts from {Xy y, Zy t) and moves 
with the velocity of light along a straight line whose direction 

cosines Z, m, n are proportional ^ ^ function 

/ will remain constant aloug its path, and consequently if the 
point once lies on the moving curve /= 0 it will always lie on 
this curve. It should be noticed that the function </> and its 
first derivatives with regard to Xy y, Zy t all remain constant along 
the path of the moving point. 

The case in which no such function (j> exists may be of 
importance in future developments of the subject ; this case 
has not yet been discussed in any detail. 

Two methods of solving equations (279) are known, but 
they are not really distinct. In the first method the functions 
a, yS are defined by equations 

(a, +[y-r, (a, /3)P +[z-i; («, [< - t («, 

I (a, /3) [« - f (a, ^)] + m (a, 0)[y-v (a, /3)] 

+ n (a, /3) [^ - f (a, jS)] = C‘‘p (a, yS) [( - t (a, /3)], 

(280), 

where y, f, r, Z, m, n, p are arbitrary functions satisfying the 
relation 

Z^ 4- 4- jol 

The functions a, may evidently be replaced by two other 
functions a\ /8' defined by equations such as 
a'^F{cCy/3)y /3'=:(?(«,^); 

consequently we may without loss of generality introduce a 
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further relation connecting the functions Z, m, Uyp. The 

relation we shall choose is 

7 


, 0^ 0i; 0f o 9t 



When OL and jS are defined in this way we can obtain the simple 
specifications of two types of electromagnetic fields ; one type 
is obtained by writing 




.(2b2), 


_ 

and two similar equations. A second type of electromagnetic 

field may be derived from the potentials* 

. _lSg(a,^) _mSg(a,$)^ 

* FS-QR’ ^ F8-QR 


where 


r)Sg(i,R) 

' FS-QR’ 


- ^ c pSg (a, 0) 
FS-QR 


.(283), 




^ dl3 3/3 3/3 3/8 


„ , 3f 3v 3? 

R = l ^ + m -- + n~ ■ 
da da da 








dm , 9 dll . 

5g (»-.))- 3^ 


It is easy to see that these potentials satisfy the relation (267). 
To prove that they are wave-functions we remark that 

P^R ^ - ’'>! 

pm 0 (of, /8) 07^ 0 (a, 0p 0 (of, ^)1 

Now it has already been proved that an arbitrary function 
of a •and /9 multiplied by one of the Jacobians represents a 

* We c£^ also call these potentials L^;, L,, A and derive an electro- 

magnetic field from them by the method of § 4. 
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wave-function, hence A^; is a wave-function. In a similar way 
it can be shown that the other potentials are wave-functions. 

It should be noticed that the electromagnetic field specified 
by the potentials (283) is conjugate to the field given by (282). 
To prove this we observe that 



and there are similar expressions for A jj, A <!>. Hence for the 
electromagnetic field specified by the potentials (283), we have 

M ' _ ^ (^’ ”) . ^ jQ) . i 9 (u > «) , i 9 (i>, jg) 

* d {y, z)~ cd{x,t) ~ c d (x, t) ' 

The relation (3) is now seen to be satisfied in virtue of two 
equations of type 

0 (a, ^) 0(u,a) 0_(a,jS) d(a^0) d(u, a) 

0 ( 3 /, z) 0 (O), t) 0 {z, x) 0 (y, t)'^d (x, y) 0 {z, i) 

4 . ^ ^ «) . S (a, yQ) 0 (^, a) ^ ^ 

0 (a;, 1 ^) 0 (y, ^) 0 (y, t) 0 (^, ^) 0 (z, t) 0 (^, y) 

(284). 

It should be remarked that the vectors J?, JI in both fields 
generally become infinite when FS — QF = 0. This equation 
is certainly satisfied when (x, y, z, t) lies on the moving curve 
defined by the equations 

^ = y = 77(a, /3), ^=5'(a,/3), f = T(a, ^)...(285). 

In some cases this curve may reduce to a moving point, 
as for instance when t;, r are independent of ^9. 

It is evident that the quantity PS — QR is usually zero for 
space-time points which do not lie on the moving curve (285). 
If, however, we regard Z, m, 71, p as complex functions of the 
type <j> (a, /C) -f (a, yS), the equation PS = QR will generally 
give rise to two distinct equations connecting x, y, z, t, when 
we equate real and imaginary terms on both sides. Hence 
all the real singularities that are defined by PS = QR will 
generally lie on one or a number of moving curves.' 

So far we have said nothing about the choice of ^ suitable 
pair of roots of equations (280). In general we cannot expect 
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them both to be real and it is difficult to lay down rules which 
will enable us to pick out just one a and just one /3 in all cases. 
To proceed further we must consider some particular examples ; 
before we do this, however, it will be worth while to point out 
that if we assign given complex values to a and /8 the equations 
(280) will generally determine two real points x, y, z, t, but 
in special cases they may give oo * space- time points which can 
be regarded as the consecutive positions of a moving point. 
Thus we have an interesting specification of the real points 
in space by means of two complex quantities. If we assign 
a complex value to the corresponding space-time points 

2/> ^ generally lie on a moving curve which travels with 

a velocity not greater than that of light. Hence in the para- 
metric representation of Xy y, Zy t in terms of the complex 
quantities a, y8 the loci a = const., /3 = const, are generally 
moving curves which may alter in shape as they move but 
never travel with a velocity greater than that of light. 

Matters are somewhat different if a or sortie function of 
a and /3 is always real when Xy y, Zy t are real. This case will 
now be illustrated by a particular example. 

Examples. 1. Prove that the ratios of the Jacobians 

3 (g, fi) d (g, /jj) 

3(y, 0)’ d(z,x)' 

are functions of a and /3. 

2. Prove that the ratio Q/S depends only on c^and /3. 

,3. Obtain the general solution of equations (279) by taking .r, y, g, 
as new independent variables. 

§ 44. Projection of singularities from a moving point, 
second model of a corpuscle. 

Let us now suppose that f, 97, r are independent of /3 and 
that T == a. We may then define a uniquely by restricting 
it to be real and introducing the inequalities 



To obtain a single value of y8 we may assume m, n, p 
to be' linear functions of /3. Consequently we may put 
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where U, li, etc. are functions of a which satisfy the equations 

+ V = c’po’ 1 

+ + V = c’p,’‘ \ (288). 

Ui + + w«Wi = c'poPx ) 

It is easy to verify that the wave-equation is satisfied by 
a function of type 

thus we have a generalisation of the theorem of § 41. 

We shall assume that ^o> Pa are real and that some 

or all of the quantities Pi are complex. It is easy 

to see that if we assign a real value to a and a complex value 
to yS the corresponding space-time points (a?, y, .sr, t) can be 
regarded as the successive positions of a point which starts 
from the point 

a: = n«). 2 /’=v(«). ■ 2 '=?(a). < = « (289) 

and moves with the velocity of light along a straight line 
through this point. There is clearly just one line through 
this point for each complex value of fi and vice versa. If we 
consider all the points in space at a particular time t we can 
specify each point uniquely by a real parameter a and a real 
or complex parameter yS. 

Let us now consider the electromagnetic field which is 
specified by the potentials 

= = .I> = ^ (290), 

where I — — Iq, m = ySmj — mo, n = ySni p — ySpi —^o and 

f is an arbitrary function of a and y8. These potentials are 
derived from (283) by putting Q = 0. 

After a long calculation we find that the component of the 
electric fofce along the radius from (f, a) to (^, p, z, t) is 


-^L 


To obtain an electromagnetic field in which there* is a 
constant electric charge associated with the singularity 
(?) ?) «)) we assume that po = =/= 1 and that 
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^ ^ + Wl 


dri 

'da 


??_ 

da~' 


7 


- (^iw M 

\da) vaoy \da) 




...( 291 ). 


The expre.ssion for the radial electric force then becomes 




P2 

Comparing this with the expression for the radial electric 
force in the case of an electromagnetic field with a simple 
singularity (^, rj, r), we see that there is a constant electric 
charge 47r/c associated with the singularity (f, r), a). 

It should be mentioned that the second of equations (291) 
is a consequence of the other equations satisfied by Iq, nio, no. 
To prove this we take the axis of in a direction parallel to 
the velocity of the singularity (^, r), f) at time a. We then 

have for this instant ~ ^ - 0. If, moreover, we choose the 

fia da • 

axis of y in such a way that tii = 0, we may satisfy the first 
of equations (291) by writing 

~ = c cos ffy li- c sec dy nil == io tan U = c cos 

da 

nio = 0, no == ± c sin 

and then it is clear that the second of equations (291) is 
satisfied. » 

Let us now write U — c cos 0, mo = 0, fio = + c sin 
= c sec 6y mi = — ic tan 6, iii = 0, x — ^ = X, y •^7j=::Y, 
z ~ f = Z, t — T=T; then it is easy to see that if 

Y-n,Z, 8 = T -l,X --m, Y-n,Zy 

U = T — lo X — Y - no Z, U = T—loX ■^moY’-rioZ, 

we have cos^ 6 88 = UU, U U ^ — 2P. 

% 

Now since ^8=- U it follows that the potentials (290) become 
infinite when IT = 0, i.e. when 

l/Q n,o 1 

When a is given these equations are satisfied by a point which 
starts at (f, r^y a) and moves with the velocity of light along 

B. 9 
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a straight line whose direction cosines are proportional to 
Wo, This line makes an angle 6 with the direction of 
motion of the point (^, r?, a). 

The electric and magnetic forces in the electromagnetic 
field derived from the potentials (290) are generally infinite 
ioT U —0 and so our field possesses a number of singular points 
which are projected from the moving point t;, a) and travel 
along straight lines with the velocity of light. It should be 
remarked, however, that if we retain only the real parts of the 
potentials (290), the projected singularities disappear as soon 
as the singularity Tjy C moves in a straight line with uniform 
velocity and 1, m, n are independent of a. The field then 
becomes identical with that derived from Lienard’s potentials. 
To prove this we shall show that, on the above assumptions, 
the field derived from the potentials 




P ^ 




m — r)' 


A,o = R , <E)» = R (292), 


is everywhere null, R being used to denote the real part of a 
quantity following it. In the first place we remark that 
we now have 


R = Zf ' -f mrj' + Tif' — c^p = ~ 


Hence 


2A 


.3a 

% 




_s d {U\^Sd /U\ „ 3 ,, p, 

P3a; (-Sf/‘^P3a; (<Sf) 

U d . S8 ^ d ,, p, 

CT 3 , U - 3 ,, p. 


Now since U U 2P = 0, it follows that U/P is a function of 

U/U and so J./, — - 4>® are the derivativ/*3 of a single 

c 
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function, consequently the electromagnetic field derived from 
these potentials is everywhere null. 

Summing up our results we can say that when the 
conditions (291) are satisfied, the electromagnetic field derived 
from the potentials (290) contains a point charge which moves 
with a velocity less than that of light ; attached to this point 
charge there is a certain curve which becomes the locus of 
a series of moving point singularities whenever its form differs 
in any portion from a straight line or its direction changes. 
The form of the curve at any instant is subject to the condition 
that the points of the curve can be regarded as having been 
projected from the moving charge at different instants, the 
direction of projection being partially determined by the law 
cos 6 zzzvjc where v is the velocity of the point charge, and 6 is 
the angle between the direction of projection and the direc- 
tion of motion of the point charge. 

We may now obtain a new model of a corpuscle by con- 
sidering an aggregate of elementary fields of the type just 
described, the point charges and exceptional curves being 
nearly coincident. If we write de for the charge associated 
with one of the elementary fields we may obtain a field in which 
the electric and magnetic forces are finite by a suitable process 
of integration. According to this idea a corpuscle has a kind 
of tube or thread attached to it. When the motion of the 
corpuscle changes a wave or kink runs alorlg the thread ; the 
energy radiated from the corpuscle spreads out in all directions 
but is concentrated round the thread so that the thread acts as 
a guiding wire. This theory of radiation is in some respects 
similar to that given by Sir Joseph Thomson in his theory 
of the Rontgen rays*. It is in accordance with his idea that 
the energy may be concentrated round certain points of the 
wave-front. ^ 

The following figure indicates roughly the changes in 
the form of a tube which always lies in one plane and is 
attached to a corpuscle performing a simple harmonic motion ; 
it is seto that a type of progressive wave travels along the tube. 

* Electriciiy and Matter ^ London (1904) ; Phil. Mag. Vol. 19 (1910). 

9—2 
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In Sir Joseph Thomsons theory of the Rontgen rays the 
kink in the tube of force becomes longer and longer as it 
recedes from the charge. A similar remark applies to the 



oscillations of the thread attached to our point charge. This 
phenomenon may be due entirely to the fact that the tube of 
force and thread extend to infinity. If we suppose that the 
tube or thread does not extend to infinity but ends at some 
other point charge, the circumstances of the motion will be 
different. If in this case we treat the thread as a singular 
line of an electromagnetic field and suppose that it is given by 
an equation of the form 

/(«,y3) = o 

where a and /9 are functions which satisfy (279), we must 
conclude that there is no function of type F (a, /S) which is 
a real function of x, for if this were the case the moving 

thread would be the locus of points travelling in straight lines 
with the velocity of light and would consequently extend to 
infinity. 

Electromagnetic fields with moving point charges joined by 
singular curves which do not extend to infinity have not yet 
been obtained, but I think there is some hope of deriving them 
by the general methods of § 43, when the quantities a, ^ are 
both complex. 

Examples. 1. Discuss the properties of the electromagnetic fields 
that can be derived from the potentials 
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respectively and determine the lines of electric and magnetic force. 
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2. Prove that if a and ^ are defined as in § 44, the electromagnetic 
field specified by the equations (282) is conjugate to the field specified by 
the potentials (290), or to the field specified by Lidnard’s potentials. 

§ 45. Electromagnetic fields with singularities moving 
with velocities greater than that of light. 

Some of the preceding analysis holds and provides us with 
solutions of Maxwell’s equations when the velocity of the 
primary singularity (^, r), r) is greater than that of light, but 
in the case of a field specified by potentials of type (266) 
a transition from a velocity less than that of light to a velocity 
greater than that of light does not seem to be physically pos- 
sible on account of the occurrence of infinite values of the 
electric and magnetic forces in the critical case. Moreover, it 
is difficult in the general case to give a rule which will enable 
us to pick out just one root of the equation (265). An interest - 
ing type of field may, however, be obtained by a process of 
summation over some of the roots of the equation*. 

The case of infinite velocity is of some interest, for then 
we obtain electromagnetic fields with singularities along a 
fixed curve at a given instant of time. The following example 
indicates that the case in which the primary singularities are 
imaginary may be associated with another case in which they 
are real. 

Consider the two equations ^ 

(x—a cos clY + (y — a sin aY + 
x^-\-y^ + {z — ai cosh ^Y — sinh )8)S 

and write in analogy with (263) 

1 / == a cos a (y — a sin a)— a sin 0 L{x — a cos «)= a (y cos a — a? sin a), 

Vi = ia sinh -- ia cosh /3) — iacosh /8 {ct — ia sinh 

= ia (z sinh ^ — ct cosh fi). 

We evidently have * 

4i/2 = 4a^ (iP* + y^) — (x^ -h z^ — a*)^ 

(2^ - cH^) 4- (^ + y' + 3^' - cH^ - 
hence* • i/* 4- = 0. 

* A more“^omplete discussion is given in G. A. Schott’s Electromagnetic 
Radiation. 
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It is easy to verify that the function 

is a wave-function and that the functions a, ^ are suitable for 
obtaining electromagnetic fields by the method of § 5 ; they are, 
in fact, the functions considered in § 36. 


§ 46. Second solution of the fundamental equations. 

The fundamental equations (279) are also satisfied when 
the functions (a, /3) are defined by the relations 

Z (a, m (a, y -f n(a, /9) t + 0] 

X (a, )8) ic 4- (a, y +N{(i, ^)z - c^P (a, yS) ^ + (?(a,^) = O) 

(293) 

where Z, X, etc. are arbitoxy functions satisfying the relations 
P 4* 4 71^ = &p\ X" 4 -MH = c^P\ 

IL -\-mM ■\-nN = c^pP (294). 

If now 


dm 

da 

dn 

cH^/ + 
0a 

0a’ 

dm 

dn 


^<7 

■0yS’ 

dM 

dN 

».3P 

0G 

da 

da 

oa 

+ 0a’ 

dM 

djsr 

o.SP 

dG 

3/8 



'^d/3’ 


we can derive an electromagnetic field from the potentials 

* = 

pv AW — pv Aw — fJLV 


where f is an arbitrary function of a and /3, and 


,SL , dM dN „ dP 


XHT — flV 

...(295) 


'...(296). 
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To verify that these potentials satisfy the wave-equation and 
the relation (267), it is sufficient to remark that 

/--ij _ _ 9^ 9 ^ dN d {a, 0) dP d («, \ 

^ ^ 0/3 3 (a;, y) 0/3 d (x, z) d/S d (x, t) 

^idMd(a,0) idNd(ct,0) . 3P3(a,/8)[‘**^ '''^- 
” c 9/8 3(^,0 “ c 0/3 0 (y, t) ~ 9/8 0 (y, z)j 

When similar expressions are obtained for Ay, Ag, 4> it is clear 
that the relation (267) is satisfied. Other types of electro- 
magnetic fields are obtained by writing a instead of /8 in 
equation (296) or by writing I, m, n, p in place of Z, M, N, P 
in (295) and (296). If A^^, Ay, Ag, <I>' are the potentials 
obtained in either of these ways we have clearly 

A^ A^ + -f .4, - dxD' = 0 (298). 

I have noticed that this relation is often satisfied by the 
potentials of two conjugate fields. 

Another type of electromagnetic field may be obtained 
by the method of § 5. 

Some particular cases of the preceding theorems may be 
deduced by contour integration. To illustrate the method let 
us suppose that I, m, ... are functions of a parameter a 

which satisfy the relations (294). If we regard these quantities 
as independent of x, y, z, t, the contour integral 

1 [ F{Lx^' My Nz oi) ^ 

”” ^iri J lx -h my -^nz — c^pt g 
will represent a wave -function. Now let us suppose that the 
contour encloses only one root of the equation 

xl (a) -f- ym (a) 4- zn (a) — cHp (a)-hg(a) = 0 .. .(299), 
and that the numerator is finite and single- valued within the 
contour and on its boundary, then by Cauchy’s theorem the 
value of the integral is generally 

(300) 

where /8 = (a) + yM (a) -h zN (a) — cHP (a) -f 0 (or). . .(301) 

and 1 is*defined in the same way as before. 

We ha^e then the result that the function (300) satisfies 
the wave-equation. This is a particular case of the general 
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theorem of § 43 and is a generalisation of a theorem due to 
Forsyth*. In Forsyth’s work the functions L, M, N, P are 
assumed to be the derivatives of m, p with regard to a. 
The functions a, ^ are evidently particular cases of the func- 
tions that have already been defined and so may be used to 
construct an electromagnetic field by the method of § 5. An 
interesting electromagnetic field may also be derived from the 
potentials 



A 


_ wf (a) 



\ 


...(302); 


it is easy to verify that the relation (267) is satisfied. The 
case in which m, n, p, g are all real functions of a is un- 
interesting because then our potentials become infinite for 
00 * space-time points which lie in oo ^ planes. When, however, 
I, m, n, p, g are complex functions of the type <f> (a) + iyjr (a), 
the singularities of the electromagnetic field generally lie on 
a moving cur\ie. 

It should be remarked that when a and 13 are defined by 
the equations (299) and (301) a function of the type 

Xl=F(a, /S) 

satisfies the wave-equation, and at the same time satisfies the 
differential equation 


\dz) cAdt) ' 

This is a generalisation of a theorem due to Forsyth f and 
Jacobi J. A solution of Maxwell’s equations may be derived 
from the potentials 

Aj,® = mo(a), A/=no(a), = C2)o(a) ...(303) 


if 


jdlo . dmo . 9?io 
9a 9a 9a 


^« = o 


da 


dl 

'da 


dm 

da 


^ 0 + niQ— + ^ 0 ^ - c^po 


dn 

'^a 


= 0 . 


* Messenger of Mathematics j Vol. 27 (1898), p. 138. The tKeor6m is 
obtained in another manner by Forsyth, 
t Loc. cit, 

X Werke, Bd. 2, p. 208. 
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One way of satisfying these equations is to put 
mo = m, = Po—p> 

An interesting type of wave-function may be obtained by 
a generalisation of a method due to Schottky 
Let 


ao (cl) da, (a) = /m (a) da, Cq (a) = / 7/ (a) da, 

^0 (^) = ip («) da, eo (a) = jg (a) da, 
then ft = (a) + (a) 4- zCq (a) — c^i5do (a) + Cq (a) 

is a wave-function whose derivatives with regard to Xy j/, 2^, t 
are all functions of the single variable a. The function a is 
supposed to be defined by equation (299). 


Example. If Z=l — sa, w = s-4-a, w=i(«-a), ^ = 0 ^, where « 

is a constant, the function X vanishes when (.r, y, 2 , lies on the moving 
curve 


l -«2 

l+> ’ 


25 , , 


S (52 — 1?2) — 25 


Z~^V. 


A point for which v is constant moves in a straight line with the 
velocity of light. 


§ 47. A wave-function with a fixed curve of singu- 
larities. 

Let T be defined in terms of Xy y, by the equation 
ar = xp (t) 4- yq (t) 4- zr (r) 
where p* 4- 4- ^2 = 0 and a is a constant. 

Let a — xp (t) — yq' (r) — zr (r), 

then ft = i f(6y r) is a wave- function “t*. 

A solution of Maxwells equations may be derived from the 
potentials « 

A.^lf{e,T), Ay^lf(e,r\ A,^'-f{0,T), <I> = 0. 

* Berlin. Sitzungsherichte (1909), p. 1152. 

t This ♦esult was derived from a theorem given by Prof. Forsyth, loc. cit. 
If we put d =/3, r=o the functions a, can be used to obtain an electromagnetic 
field by the m^hod of § 6. These functions are, of course, particular cases of 
the functions defined at the beginning of § 46. 
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To obtain a real electromagnetic field we must retain only 
the real parts of these expressions. 

It can be shown by putting 

(u, v) -f iPa v\ q^Qi (Uy v) + tQa {Uy v)y 

r = Pi {Uy v) 4- ^‘P 2 {u, v)y tt = J.J {Uy v) 4 iA^ {Uy v)y 

that the electromagnetic field has generally a fixed curve of 
singularities. In the special case when 

P = g=2T, r — i(l + T% ^i = A, A^=^ky 

the fixed curve is the circle \ ^ A. 

§ 48. Cylindrical wave-fanctions with moving singu- 
larities. 

If we define t in terms of x, y by the equation 
? (t)? +{y-v (■>•)]“ = c" (« - t)», 
and define I (t), m (r), p (r), so that 

(r) -h m?/' (t) = c^py 
p 4- == &p^y 

the function v I {x - m {y — r}) — c^p {t -- t) 

1 

is such that f (r) 

is a wave-function. 

In particular, if ^ = 0, i = 1, m = i, = 0 we obtain 

Poisson’s wave-function 



Another interesting result is that if F {x, y, t) is a homo- 
geneous function of degree ^ which satisfies the wave-equation, 
the function 

• n = F(x-^,y-T],t-T) 

•also satisfies the wave-equation. 

If <r = f ' (t) (a: - ^) + v' (t) {y-rf)- c" {t - t), 

the function 


1 
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satisfies the wave-equation only when f ^ ^'2 _ ^2 When 
iz is written in place of ct these wave-functions can all be 
regarded as solutions of Laplace’s equation. 


EXAMPLES. 


1. Let 


x^xfir)- Z=zf{r)- jy{s)C(s)ds, 

y=yf{r) - JV' (») r, («) ds, i= C/(t) - J^/' (*) ; 

where r is defined in terms of a?, y, z, t by means of equation (265) and the 
inequality Prove that if an electromagnetic field H) is such 

that 


1 0r /I 0r . 0r 


.0r 






' dz 




with similar equations, where M^H±iE, then an electromagnetic field in 
the variables Jc, y, I, t can be found such that we have identically 
(2/1 z)-\-Myd (2, x) 4 - M^d (a?, y) + icM^^d (^, t) 

+ ic Myd (y, t) + icMgd {z, t) 

= Mj^d (y, z)-\-Myd (2, x) 4 - M^d (.r, y) __ __ 

+ icMj^d {x^ t) + icMyd (y, if) + icM^d {zy t), 
where d{y, z) denotes dybz-dzhy and dx, bXy etc. are two independent 
sets of increments of the variables. * 

2. Prove that if/'(r) is always positive the variable t increases with t 

Show also that if (|, 7, r) correspond to (f, 7;, f, r), the point (|, () 

moves along a curve f with a velocity less than that of light and that the 
velocities at two corresponding points of r and r are the same in magnitude 
and direction. 

3. Prove that the conditions imposed upon the electromagnetic field 
in Example 1 are all satisfied in the case of the field specified by the 
potentials (266). Hence show that the transformation transforms the 
field of an electron moving along the curve r into the field of an electron 
moving along the curve r. 

4. » Prove that the conditions of Example 1 are also satisfied for any 
electromagnetic field of type (282) where a and ^ are defined as in § 44. 
Hence show -Jhat a field of this type is transformed into another field 
of the same type associated with the curve r. 
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5. Prove that the conditions of Example 1 are equivalent to only two 
conditions which imply that the field {E, H) is conjugate to any electro- 
magnetic field of type (282). 

6. Prove that electromagnetic fields of ty^ies (302) and (303) can lie 
transformed into fields of the same general types with the aid of the 
transformation 

^ I {s) f{s) dSf 0=2 4- n {s)f{s) ds^ 

1/4- J m{s)f{s)d8, 


y 


7. If r be defined in terms of y, 0 , t by the equations 

[x - I (a, t)]2 4-Ij/ - V («) ^)f + [2 - f (a, r)P 6 (a, r)p, 


it satisfies the partial difieren>ial equation 


(^-7 


8. If 


and 


a*=pcos<^, y=psin0, 


sinh = - 


2ap ’ 


the function [x^ 4- 4- (2 ± a)^] ^/ [<#> + 

satisfies Laplace’s equation. 


9. Prove that if in*the last example we write 
+ - he''* 

and integrate with regard to u between 0 and 27r we can obtain a potential 
function which is zero outside the tube \<f> + iylr\ = k. 


10. Particles are projected in certain directions from the different 
positions of a moving electron and travel along straight lines with the 
velocity of Ifght. Prove that if the law, according to which the direction 
of projection varies with the velocity of the electron, be suitably chosen 
**the particles will at each instant form a line of electric force in the 
electromagnetic field due to the moving electron. 
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MISCELLANEOUS THEORIES 

§ 49. Kirchhoffs formula and its extensions. 

An important solution of the wave-equation is embodied in 
Kirchhoffs formula* which is usually interpreted as the mathe- 
matical expression of the principle of Huygens. This formula 
has been extended by Lovef and Macdonald J so as to give 
a representation of an electromagnetic field outside a surface 
in terms of the electric and magnetic forces tangential to 
the surface. In Macdonald s formula it is the time derivatives 
of E and H that are so expressed. Tonolo§ has given a 
formula in which E and H are expressed in terms of their 
surface values. The formulae are given in examples 3 — 5 at 
the end of this chapter. 

When the surface is a sphere Kirchhoff s formula reduces to 
the formula of Poissonll (Ex. 5) which enables us to find a wave- 
function which satisfies the conditions 

n y, z\ y> 

Poisson’s formula may be used to derive the theorem^ that the 
mean value of a wave-function u over a sphere of radius cr at 
time t is equal to the mean value of u at the centre of the 

* Berlin, Ber, (1882), p. 641; Wied. Ann. Bd. 18 (1883); Ges. Ahh. t. 2, 
'y 22. Simple proofs of the formula have been given by Beltrami, Rend. Acc, 
Lin,: Rom. (5), t. 4 (1895) ; Larmor, Proc. London Math. Soc. Ser. 2, Vol. 1, p. 1 ; 
Love, Ibid. p. 37 (1903) ; Lamb, Hydrodynamics y 2nd edition (1906), p. 477 ; 
H. A. Lorentz, The Theory of Electrons, p. 233 ; E. Laura, II Nmvo Cimento 
(1913). 

t Phil. Trans. A, Vol. 197 (1901). > 

X Electric Waves, p. 16; Proc. London Math. Soc. Ser. 2, Vol. 10 (1911), 
p. 91; Phil. Trans. A, Vol. 212 (1912), p. 295. This theorem gives an 
analytical specification of a generalised Green’s equivalent layer. See p. 29. 

§ Annali di Maternatica, Ser. 3, t. 17 (1910). 

II The details of the calculation are given by Love, loc, cit. A simple proof 
of Poisson’s formula is given by Lamb, loc. cit. p. 471. 

II Cf. Ravleigh’s Sound, appendix, and H. Bateman, Anier. Journ. (1912), 
where some OTher theorems of a similar kind are given. 
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sphere during the interval t — T to ^ + t. The function u is 
subject to the conditions in KirchhofTs theorem. 

When the function u is independent of z, Poisson’s 
formula reduces to Parseval’s formula for a cylindrical wave- 
function. Vol terra* has extended Parse val’s formula so as 
to obtain a two-dimensional analogue of Kirchhoff’s formula. 
His formula indicates that the propagation of cylindrical 
waves is essentially different in character from that of 
spherical waves. In the three-dimensional case the value of 
a wave-function u (x, y, Zy t) at a point {x^ y, z) at time t 

is completely determined by the values of u and ^ over a 

concentric sphere of radius cr at time t — t. In the two- 
dimensional case, on the other hand, the value of u (Xy t/, t) at 
a point {Xy y) at time t is not determined by the values of u and 

^ over a concentric circle at time t — r. To find u (Xy y, t) we 


must know the values of u and ^ over a series of such circles 

ot 


in which the radius ct varies from 0 to some other value ct,. 
The essential difference between the two cases may be attri- 
buted to the fact that in the three-dimensional case the wave- 


function for a source is of type while in the two- 

dimensional case it ^is of type ^ ^ wave 

does not leave the region undisturbed after it has passed, but 
has a tail or residue f. 

When w is a periodic function of ty Kirchhoff’s formula may 
be replaced by the simple formula of Helmholtz In this case 
there is an analogous formula for cylindrical wave-functions, the 
function KQ(ipk) taking the place of 


§ 50. Green's Functions. 

« The solution of a problem in which a periodic wave-function 
is to be determined from a knowledge of its behaviour at 

* Acta Math. t. 18 ; Lectures at Clark University (1912), p. 38. 
t See Lamb’s Hydrodynamics y p. 474. 

X See also J. Hadamard, Bull, de la Society math, de France y t. 28 (1900), 
p. 69 ; J. Larmor, Proc. Lojidon Math. Soc. (2), Vol. 1 (1903), p. li. 
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certain boundaries can be made to depend on that of an 
auxiliary problem, viz. the determination of the Greens 
function*. 

Let 0 2/, ; .^’i, 3/1, Zi) be a solution of Au ‘\-k^u = 0 with 

the following properties; It is to be finite and continuous, 
as also its first and second derivatives, in a region bounded 
by a surface S, except in the neighbourhood of the point 
where it is to be infinite like cos A;r/47rr, when 
r ^ 0. At the surface Sy G satisfies some boundary condition 

such as (1) w = 0 or (2) = 0. 

Adopting the notation of Plemelj f and Kneserj we shall 
denote the values of a function <\> (^, f) at the points {xy y, z)y 
(xiy 2/1, Zi) respectively by (f> (O)and (f> ( 1 ). The Green’s function 
is then denoted by the symbol (r(0, 1). The importance of the 
Green s function depends chiefly on the following theorem. 

Let </) be a solution of 

A</> + + f{xy y, Z)=^ 0 y 

which is fini ^e and continuous, together with its first and second 
derivatives, through the interior of the region and satisfies the 
same boundary condition as (r (0, 1), then 

= Jjjf (0)0 (0,1) <l^d(/ch. 

This theorem is proved by applying Green s theorem to the 
region between a small sphere 2, whose centre is at (^i>yi,2'i), 
and the surface S. For since 

JJI (<f>AG - GA<f)) dxdydz = f f ^ ^ 



we obtain the required relation by making 2-^0 and using the 
boundary conditions. 

* This function was first used by Green in the solution of a problem of 
electrostatics, Essay on the application of mathematical analysis to the theories 
of electricity and maynetismy Nottingham (1828) ; Math, Papers y p. 31. 

t Monatshefte fur Math. n. Phys. (1904) and (1907). 

it Die Integ ralgleich ungen und ihre Anwendnng in der mathematischen Physiky 
§ 31, Brunswgk (1911). 
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If g{%0) is the Greens function for the same boundary 
condition as (?(2, 0)but for k — a, we must also surround the 
point (2) by a small sphere when we apply Green's theorem 
with <f> (0) — g (2, 0), We then obtain the equation 

gi2,l) = 0 (2, 1) (2, 0) G (0, 1) dxdydz. 

This important relation indicates that ^(2,1) is the solving 
function of an integral equation of which G (0, 1) is kernel and 
vice versa. The theory of integral equations tells us that when 
(? (0, 1) is given there may be certain singular values of 
for which g(2, 1) is not finite. These are the values of for 
which the homogeneous integral equation 

<f> (1) = (a- — IJJ ^ W ^ (^> 1) dxdydz 

possesses a continuous solution <^(0) which is different from 
zero. Formula (2) indict. tes that for such values of the 
differential equation A(^ + k^^^0 possesses a solution satisfy- 
ing the boundary condition and the other conditions imposed 
on <^. The solutions of this type are of great importance in the 
theory of sound and have been discussed by many writers*. 

If we put /(O) = (o-^ — A;2)^(0, 2) and proceed as before, 
Green’s theorem gives 

g (1, 2) = (7 (2, 1) -(k^-a‘)Jllg (0, 2) G (0, 1) dxdydz. 

Putting (7 = and comparing this with the previous equation 
we get 

Hence the Green’s function is a symmetric function of the 
coordinates of the points 1, 2. When the boundary condition is 

^ = 0 this result is equivalent to Helmholtz’s theorem f. 

Since G is a real symmetric function when A; = 0 it follows 
from the general theory of integral equations that there is 


* See especially Lord Rayleigh, r^or^o/<SoMncZ, Vol. 2; Pockels,Z)ig^r»ritcZZe 
Differ entialgleichung Au-]-k’^u = 0 ; A. Sommerfeld, Encyklopddie der Math, 
Wise. Band ii. 1, Heft 4, p. 540. 

t Cf. Rayleigh’s Sounds Vol. 2, p. 131. 
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at least one real singular value of cr® ; that all the singular 
values are positive may be deduced at once from the 
equation 


<fx^ dxdydz 4*^4* dxdydz 


^ III [CtJ Ci)' + Cif] 

The Green’s function is usually obtained in practice by finding 
a suitable expansion in terms of elementary solutions of the 
equation A?/ -f k~v — 0. This method is explained in Heine’s 
Kugelfunktionen and many examples of Green’s functions are 
given for the case k = 0. The general case has been discussed 
at length by A. Sommerfeld* who also obtains a number of 
definite integrals which represent Green’s functions. These 
expressions lead to interesting generalisations of Fourier’s 
theorem. 


The problem of electrical oscillations in a cavity has been 
discussed by Weylf. With the aid of a generalisation of the 
Green’s function, viz. a Green’s tensor, he obtains a number 
of inequalities satisfied by the periods of vibration. 

The Green’s function for the equation ^u+k'^u = 0 can 
theoretically be found when the corresponding Green’s function 
for the equation Aii — 0 is known. Considerable progress has 
been made in the theory since the appearan^pe of Heine’s work 
and so a few references to recent literature will be useful J. A 


* Phys. Zeitschr. Bd. 11 (1910), p. 1087 ; Jahrcshericht tier deutsch. math. 
Vercin, Bd. 21 (1913). 

t Math. Ann. Bd. 71, p. 441 ; Crelle, Bd. 141 (1912). 

For the determination of special Green’s functions see E. W. Hobson, 
Cambr. Phil. Trans. Vol. 18 (1899), p. 277; H. M. Macdonald, Ibid. p. 292, 
Proc. London Math. Soc. Vol. 26 (1895), p. 161 ; A. G. Greenhill, Proc. Cambr. 
Phil. Soc. Vol. 3 (1880) ; J. Dougall, Proc. Edinburgh Math. Soc. (1900) ; 
H. S. Carslaw, Ibid. (1912), Froc. Loudon Math. Soc. (2), Vol. 8, p. Il65 ; C. W. 
Oseen, ArHv for mat. Bd. 2; C. Neumann, Leipziger Berichte, Bd. 58 (1906), 
Bd. 62 (1910) ; W. Burnside, Proc. London Math. Soc. Vol. 25 (1894), p. 94. 
For the general theory H. Poincar4, Ii€7ui. Palermo, t. 8 (1894), p. 57; 

Earemba, Ibid. t. 19 (1905) ; E. B. Neumann, Studien iibcr die Methode 
von C. Neumann und G. Robin zur Lbsung der beiden divertnufgahen der 
Potentialtheorie, Leipzig (1905) ; D. Hilbert, Gbtt. Nachr. 1904) ; M. Mason, 
Newhaven Math. Colloquium (1910) ; E. Picard, Ann. de V^cole Normale (1906), 
p. 509. * 
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good account of the developments up to 1900 is given in 
Sommerfeld’s article in the Erwyklopadie der Matheniatischen 
Wissenschaften. 

§ 51. The transformation of the electromagnetic 
equations. 

The transformations which can be used to transform any 
solution of the wave-equation into another solution or any 
electromagnetic field into another belong to a group which 
is characterised by a relation of the form* 

-f dy'^ -\-dz^ (dx^ -f dy^ 4- dz^ — c^dt^). 

The linear transformations belonging to this group are of great 
importance in the modern theory of relativity f ; two of the 
non-linear transformations have been mentioned in § 13. 

In addition to these transformations there are other trans- 
formations, involving arbitrary functions in their specification, 
which can be applied to certain types of wave-functions, and to 
certain type?, of electromagnetic fields. There are often two 
families of wave-functions to which a given transformation can 
be applied, when the transformation is of a suitable char acter ; 
each of these families may be defined by a linear relation which 
exists between the wave-function and its derivatives, sometimes 
between the derivatives alone. Some idea of the theory may be 
derived from the examples. It also happens that there is often 
a family of electromagnetic fields to which a given transforma- 
tion can be applied aud this family is defined by means of two 
linear relations between E and Hy which can be interpreted to 
mean that the field is conjugate to some definite electromagnetic 
field or family of electromagnetic fields determined by the 
transformation. In some cases these last fields are self-con- 
jugate and the transformation is applicable to them also. 

* H. Bateman, Proc. London Math, Soc, Ser. 2, Vol. 7 (1909), Vol. 8 (1910) ; 
E. Cunningham, Ibid, Vol. 8 (1910). 

t For this see A. Einstein, Ann, d. Phys. Bd. 17 (1905) ; Laue, Das 
Belativitdtspr inzip y Brunswick (1911) ; E. Cunningham, British Association 
Reports (1911) ; H. Minkowski, Gott. Nachr. (1908) ; E. B. Wilson and G. N. 
Lewis, Proc, Amer. Acad, of Arts and Sciences, Vol. 48 (1912), p. 38Q ; J. Ishiwara, 
“Bericht iiber die Relativitatstheorie,” Jahrbuch der Radioaktivitdt und 
ElektronikfBd, 9 (1912), pp. 560 — 648; L. Silberstein, The Theory of Relativity, 
Macmillan and Co. (1914). 
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The fact that the condition of conjugacy between two 
electromagnetic fields often implies the existence of one or 
more transformations depending on arbitrary functions, may 
be regarded as of some philosophical interest. 


MISCELLANEOUS EXAMPLES. 


1. Show that the most general periodic solution (valid for all 
space outside a given closed surface) of the wave-equation is 


1= J-J- (ta-i)} dS, 


where P and Q are arbitrary functions, r is the distance from the element 
of surface dS to the point where ^ is estimated, and yfr is the angle between 
r and the outward drawn normal. Show further that the necessary and 
sufficient condition that the value of given by the same analytical 


00 

expression, should vanish for points inside the surface, is that 1*= . 

(Cambr. Math. Tripoi, Part II, 1904.) 


2. Let H be a function which satisfies the wave-equation and is such 
that its differential coefficients of the first order are continuous functions 
of A’, y, 2 , t within a region bounded by a closed surface S. If either O or 

^ be given for points on the surface & there is only one function Q which 
on 

reduces to a given function /(a;, y, z) for t^t^. 

(A. E. H. Love, Proc. London Math, Soc, Ser. 2, Vol. 1, p. 42 ; 
J. Hadamard, Bull, de la Soc. Math, d? France^ t. 28 (1900).) 


3. If throughout a specified region of space and a specified interval of 
time u and its differential coefficients of the first order are continuous 
functions of .r, y, 2 and of if also the differential coefficients of the 
02 ^ 02 ^^ 

second order such as integrable, then a solution of 


the equation 


Qu-^<r{x, y, Zy t)—0 


which is valid for this region is given by the fonnula 




An 

1^ 

An 


//{Ms©- 

1 rdul 
r [dn] 

[jf ^ 



1 dr r^ll 

r dn Jl 


dS 


where denotes the normal to dS drawn 

into the spec^aed region and the integration is taken throughout this 

in o 
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region and over its boundary. The function o- is supposed to be finite and 
integrable and a quantity within square brackets is calculated at time 

(G. Kirchhoff.) 

4. If an electromagnetic field is such that the specified region does 
not contain any charges or convection currents and the value 

of the vector Jfat (^oj yo) ^o> ^o) is given by the formula 

cos no; — [Mn] COS ra; - cos m} ^ 
cos nx — \M^ cos rx - \_M^ cos m 


. r 1 • r ly- T ^ dS 

-f I [if,] cos — I [if yj cos nz} — . 

r 

6. Prove that in the same circumstances 

. 0 ir« 1 a rrr 3 /v\ a /3M 


(A. Tonolo.) 


where a=/x[i/,] — v [My\ etc. and (X, /i, v) are the direction cosines of the 
normal drawn into the region bounded by S, 

(H. M. Macdonald.) 

6. If w is a wave-function independent of z and periodic in t like 

«(^o.3^o)=^ / {^<>(*^) £ ^0 W} d,. 

7 . A wave-function which satisfies the conditions 

is given by the formula 

“=0^ (</)+«#, 

where J, g denote the mean values of /, g respectively over the surface of 
a sphere of radius ct having the point x^ y, z as centre. 

(S. D. Poisson.) 

0^Z4 'ey'll/ 

8. If t4*&ati8fiesg^ ^ “ 0f2 finite second derivatives within 

a suitable domain 

1 r d(T 9 /’ V 

2w J V (^2^ al “ y' 


p*=(a:-ar,)*+( 3 ^-yi)*, 


where 
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and cr denotes the area within the circle cut out on the plane t by the 
cylinder 

(Parseval and Volterra.) 

9. Prove that if a transformation of variables from {x, y, t) to 
(Xj y, z, t) is such that 

_ c^c[l 2 -s X 2 (^^2 (iy 2 ^ — c^dt^) + {Idx ■^mdt/-{-ndz — opdt) 

{lQdx+inQdy+nodz'-cpQdt)f 
where P + m^ + n^~p\ used to transform an 

electromagnetic field {E, H) into another electromagnetic field {E^ H) 
with an identical relation of the same type as that used in Ex. 1, Ch. viil, 
if the two conditions embodied in the relation 

i/x ( + My (71/0 ~ '^ 4 ) + Mg (^mo - hm) + i M^ (Z^t?o - l^p) 

+ i My ( 771^0 - i Mg {npQ - riop ) = 0 

are satisfied. 

Prove that the conditions can also be thrown into the form 

(Wo + TTlTTlo + TlTlo + VPo) M^ + 1 (tI/Jq + '•^P) My ± I (tTI^q + '^^OP) Mg 

~ Iq {IMy^ 4“ inMy 4- nM^ 4 - 1 (^o-^x "h 'm^My 4- ti^M^ 

and similar equations. 

10. In the last example if ^o=^> Po= “"Pj conditions 

are satisfied if Poynting’s vector is in the direction (^, tti, n). Show, in 
particular, that the transformation 

x±ni = F{^^^, z=/(r + c<)+^(r-c<), d.=f{r-\-ct)-g{r-ct) 

can be applied to an electromagnetic field in which Poynting^s vector is 
along the radius from the origin and that in the resulting electromagnetic 
field Poynting’s vector is parallel to the axis of z, A>pply the transformation 
to the electromagnetic field derived from the functions a, ^ given by 
equations (13), § 5. 

11. If a transformation of coordinates is such that 

dx^ 4 - dy^ \d^- c^di^ — {dx*^ 4 * 4 - dz^ - c^dt^\ 

where X is a function of x^ y, z, t', there is also a relation of type 
{x - xqY 4- (y - yo)^ + (-2 - (t - 

= XXo [(•*’ - +{y- Vo ? + (2 - «o)* - - <o)*]- 

(J. Liouville and S. Lie.) 

12. Prove that the differential equation 

is covariant for a transformation of the type considered in Ex. 11, 

I (S. Lie.) 
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13. Prove that if ri^ r are functions of t such that 

{(My — J/, + ictMx) ■“ i icrMy) dy 

+ (lyi/* - fifj, ± icrM^ dz + ic {i^Mx'\‘r\My + dt 

is an exact differential and € is a quantity whose square may be neglected, 
the value of if at the point y=y 4 'f 7 , f = ^+er may 

be calculated by assuming that the integral form 

Mx [dybz-- dz hy) My {dzhx— dx hz) + if, (dx by — dy bx) 

+ icif* {dx bt - dt bx) + icMy {dy bt - dt by) + icM^ {dz bt - dt bz) 
is an invariant for the infinitesimal transformation, it being supposed 
that the function if satisfies equations (10) of § 5. 


14. Let a transformation from the coordinates (^, y, /, t) to (A’,y, s, t) 
be such that 

d^2 4. d:y'2 - c^dtf^=^K {dx^+dy^ + dz^ - c^dt'^) 

0 ' 

+ {Idx^^-mdy + ndz — <?'pdt'f 

VA 

— ~ {I dx -^mdy ^ndz- c^pdt) {Vdx + m'dy 4 - n*dz - c^p'dt) 

V A 

0 

H — {I'dx 4- w!dy + v!dz - c^p^dty^ 

VA 

where O' — 

^-rrim! ’{■nn' -~c^pp\ A = 4>2~ee'; 

and O', I, m, w, jt?, l\ n\ p' are functions of x^ y^z^t\ then if 6 satisfies 
the equations 


^ de^ de^ de . 




it also satisfies 


0^2 02:'2 c 2 0^'2 * 


In the preceding equation we have 
dR dQ dL 
dy dz dt * 

_d^ ^ ^ 

^ dx dy dt ’ 

F=mn' — m'7if — n% R^lm' - Vm^ 

L^lp' ~ I'py if = mp' — m'jo, np' - n'p. 

Show that in certain cases a function of type is a solution of the 
wave-equation in consequence of the two equations imposed on 6. Discuss 


dz dx"^ dt ^ 

__0Z ^ 0_^ 

^ dx dy dz ’ 
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the case of the transformation of Ex. 1, p. 139: also the case of a tran 
formation which leaves the functions JT, F, Z of Ex. 7, p. 80 unaltered i 
form. 


15. Prove that if a function V satisfies the equation 

dW 

vanishes at infinity and has continuous derivatives except at points of tb 
curve O' where its normal derivative is discontinuous in such a way that 

the symbol Q being used to denote the coordinates of a point then 

ds 

where 

(Levi-Civita, Nwovo Cimento, 1897. 

16. Prove that 

1 

V(^ - ^ 4 . (y" Vy' + (7 - f )2 






- dXdfidvy 


' a “ -f -b 

and deduce that the integral j j J satisfies Poisson’* 


equation A F + (^, y, 2 ) = 0. 


(J. Weingarten. 


17. Prove that the equation = + 2 is satisfied by 

<*> = «■' {/(<-<]-'•)+/(< -<l + r)} f(r)/a(^)^ 

where ^2 ^ - r)^ - r* = (.r - + (y - yo)‘^ + (^ - ^o)*» 

and /o(^) is the Bessel’s function t/o(*^)* ^oy Voy h arbitrary 

(M. Brillouin, Comptei Rendus^ 1903. 


constants. 


18. Prove that a solution of the differential equation 




is given by 


I ” (X) C 08 (( Vx* - 1 ) + (X) ^ ’ 
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hence obtain a solution of the equation putting 

O TJ 

F= U{x^ t) e'~K If U-f (x), == g (x) when t = 0, the functions F, G may be 

determined by the equations 

/(^)= J" e'** if (X) d\ g {x)^ G (X) d\, 

with the aid of Fourier’s theorem. 

(H. Poincard, Comptes RenduSy 1893-4.) 


19. Prove that with the conditions of the last example 

• rai+tr dG 1 

where G{x,i; t^T)=Jt,>J{t-Tf-{x-^)\ 

(Laplace (1779) ; E. Picard, Bull, soc. math, t. 22 (1894), p. 2.) 


20 . 


Prove that a solution of the equation A4> = 


W 


+ 2 


04 > 

0 ? 


is given by the following extension of Kirchhoff’s formula : 


Itt F(^i, Xqj yoj 


:.)-///[« 


dV _ yd^ 

ht 3<J(=o 


dxdydz 


rer=fiL 

0n 

^ 0^ dr J 


dS 

t-0 


r 

T /0 V dr 

dV\ 

Vdr 

1 


=o+i 

r \ dt dll 

^dn)^ 

r^ dn 

J <-(i- 


,/dV 

dV dr\ 


dG 

aG'x-i 

-L L" 

\0a 

dt dn) 

+ ^dn[ 

dt 

SrVJ 

-yo)‘“ + (2 

-«o)^ 

e^^{t- 



1 d 
dTe 


and /q {&) is the Bessel’s function {id). 

The first integral extends over the volume enclosed by both the sphere 
/•=^i and the surface 2; when this sphere cuts the surface the second 
integral extends over the part of the spherical surface inside 2, the last 
two integrals extend over the part of 2 which lies inside the sphere. The 
normal n is supposed to be drawn into the region of integration. If 2 is a 
closed surface and {xq^ yo, bes outside, the region of integration is the 
space outside 2 and inside the sphere. 

(M. Brillouin, Comptes Rendus, 1903.) 


21. Let a, /3, o) be defined in terms of x^ y, z by the equations 
y2^=rC08aC0s(/3 + ^(i)), 7^2^ = cos a sin (j3 4-^0)), y^jj^sinacoso, 

where 72=1 —gin and Jb is a constant. A solution of Laplace’s 
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equation Au=0 is then given by w=yi^(a,3) provided F satisfies the 
partial differential equation 

. g. d^F sin^a+^^cos^a 027^ _ « 3 „ . _ . 

sin 203-2 + 4 ~ +2cos 2o^ 7 7^sin2a=0. 

Oa^ sin 2a 0/3^ da 4 

(U. Amaldi, Rend. Palermo (1902), p. 1.) 

22. Prove that the following transformations of coordinates lead to 
biftary potentials, i.e. to solutions of Laplace’s equation of type ??) : 

(1) y=»7, 

(2) 07=^008^', y=|sini', z^r)\ 

(3) ^rsTlcosi*, y=^sinf, z=r]--m(; 

(4) sin ^cos);, y f sin ^ sin 7, z=(cos^; 

(5) :r=|sin7e ^cosf, y=^sin7e ^sinf, z=^$cosrje^ ; 

where m is an arbitrary constant. The differential equations satisfied by 
F in cases (3) and (5) are 


d^F Id^F 1 9 ^’ 1 0 ^ 

(1+^2 cosec^ + + cosec* ^ ^ cot , ^ = 0 

respectively. In the other cases the differential equations are already 
familiar. 

(T. Levi-Civit^, Turin Memoirs^ (2) t. 49 (1900).) 
02 F 02 F 02 F 

23. If the differential equation ^ ^ satisfied by an ex- 

pression of type V^yf{6) where /is an arbitrary function, 6 must satisfy 
the differential equation 


/^Y /0^Y _ 

\dx) "w;* 


Prove that if we write /i<i^=cos a<f.r+sin acfy-c?^, 5?=^ cos a +24, 
y=^8ina+v, there is a relation between a, w, v. Discuss the cases in 
which a is a function of 6 and a constant respectively, and obtain the 
general value of y in each case. 


24. If in the last example u—f{a, 6\ v=^(a, 6) and we write 


9/ 

^ = rsina, 


-= — tcoso, 


g^=<rco8e, ^=<r8m,, 
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027 027 ^ 


when a, d, t are taken as 


[<^‘=os(a-0 ^g^+<rsin(a-*)^^ 


the new form of the equation 
independent variables is 

dv . . , a . dvi 

+ 1 [ ^ 

Prove that this equation can only possess a solution of type V=yF{3\ 

0T 

with ^arbitrary, if g^=0) in this case 6 is defined by an equation of 
type 

[^-|(<))P+[y-,W=[<-r(d)P, 


while the most general value of y is 

{?(d)[a;-f]+m(d)[>-i/]+«(fl)[<-r]}i 

f a + V Wb->>]-r'W[<- r] ’ 
where m, n satisfy the relation 


25. Show that wave-functions of type yf{6) may be derived from 
solutions of Laplace’s equati'^n of this type by means of the results given 
on pp. Ill, 114. Hence show that there are wave-functions of type 
y f (d) which are not particular cases of a more general wave-function of 
type y f(aj where a, are defined by equations (280). 
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